1 [x,y] = meshgrid(linspace(-3,3,201)); 2 3 %Scattering from a cylindrical cavity 4 [du, dv] = cylindrical_cavity(x,y,0,4*pi,1,1,1); 5 z = sqrt(du.^2+dv.^2); 6 mask = 1.0*(x.^2+y.^2 >= 1); 7 subplot(1,4,1) 8 surf(x,y,z./mask) 9 shading interp 10 11 axis([-3 3 -3 3]) 12 view(0,90) 13 14 %Scattering from a cylindrical inclusion 15 [du, dv,du_p, dv_p] = cylindrical_inclusion(x,y,x,y,0,4*pi,1,1,1,1/2,1/10,1); 16 17 z = sqrt(du.^2+dv.^2); 18 z_p = sqrt(du_p.^2+dv_p.^2); 19 mask = 1.0*(x.^2+y.^2 >= 1); 20 mask_p = 1.0*(x.^2+y.^2 < 1); 21 subplot(1,4,2) 22 surf(x,y,z./mask) 23 hold on 24 surf(x,y,z_p./mask_p) 25 shading interp 26 27 axis([-3 3 -3 3]) 28 view(0,90) 29 30 31 %surface waves on a convex boundary 32 c = 6.712844035888430/6; 33 B = 0.259859421131935i; 34 [du0,dv0] = surface_wave_convex(x,y,1,1,c,B); 35 subplot(1,4,3) 36 z = sqrt(du0.^2+dv0.^2); 37 mask = 1.0*(x.^2+y.^2 <= 1); 38 surf(x,y,z./mask) 39 shading interp 40 41 axis([-3 3 -3 3]) 42 view(0,90) 43 44 %surface waves on a concave boundary 45 subplot(1,4,4) 46 c = (4.000739349461448 + 0.616417207729038i)/6; 47 B = -2.741262468068052 - 16.740023608389567i; 48 [du, dv] = surface_waves_concave(x,y,0,6,c,B); 49 z = sqrt(du.^2+dv.^2); 50 mask = 1.0*(x.^2+y.^2 >= 1); 51 surf(x,y,z./mask) 52 shading interp 53 54 axis([-3 3 -3 3]) 55 view(0,90) 56 57 pause 58 59 60 61 62 63 64