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