close all

%Plottar egenfuntionerna på området till de beräknade egenvärdena

R=0.01; %Radie
phi_0=(pi/4)*1.7;   %Vinkel phi
mnvalues=[1,1,2;1,2,1];  %Bestämda värden på m och n
k=[3.1137*10^5,7.0997*10^5,7.9156*10^5];    %Bestämda egenvärden

% Skapa ett mesh i polära- och kartesiskakoordinater
r = linspace(0, R, 500);
phi = linspace(0, phi_0, 500);
[r_mesh, phi_mesh] = meshgrid(r, phi);
[x, y] = pol2cart(phi_mesh, r_mesh);

% Egenfuntioner till motsvarande m och n
E = @(m, n, r_s, phi_s) sin(n .* pi .* phi_s ./ phi_0) * besselj(n * pi / phi_0, sqrt(k(m + 2*n - 2)) * r_s);

% Plotta egenfunktionerna till motsvarande egenvärden
for i = 1:3
 figure();
 surf(x, y, E(mnvalues(1, i), mnvalues(2, i), r_mesh, phi_mesh), 'EdgeColor', 'none');
 view([0 90]);
 grid off;
 axis off;
end