clc
close all
clear all
Dr = 0.3;
Dl = 0.2;
ar = 0.6;
al = -0.6;
r=1;l=-1;h=1;Delta = 0.05;

[panelsL,panelsR,panels0,panelsF] = semicircles_two(Dr,Dl,ar,al,l,r,h,Delta);

domainplot(panelsL,panelsR,panels0,panelsF)

print -djpeg wave_bem_example_domain

alpha = 0:0.01:5;
number_modes = 5;
R = zeros(size(alpha));
T = ones(size(alpha));
% we know the solution for alpha = 0 is R=0, T=1;

for j = 2:length(alpha)
    j
    [Phi,Phi_n,refl,trans] = fixed_body_twod(alpha(j),number_modes,panelsL,panelsR,panels0,panelsF);
    R(j) = refl(1);
    T(j) = trans(1);

end

set(gca,'FontSize',18,'LineWidth',2)
plot(alpha,abs(R),'-',alpha,abs(T),'--','LineWidth',2)
xlabel('\alpha','FontSize',18)
ylabel('|R|, |T|','FontSize',18)
axis([0,alpha(end),0,1])

print -djpeg wave_bem_example_RT