clear
close

load u1_d;
size(u1_d)
n1=u1_d(1,1)+1
r1=u1_d(1,2)
n2=u1_d(2,1)
if n2~=0
   n2=n2+1
end
r2=u1_d(2,2)

subplot(1,2,1)
for j=1:n1
  r(j)=u1_d(j+2,1);
  u(j)=u1_d(j+2,2);
end
rmax=max(r);
umax=max(u);
plot(r,u)
title('u1 (Zoom Up)')
set(gca, 'box','on')
axis([0 rmax 0 umax*1.1])
xlabel('r (0 - R1)')
ylabel('integrand')

subplot(1,2,2)
for j=1:n1+n2
  r(j)=u1_d(j+2,1);
  u(j)=u1_d(j+2,2);
end
rmax=max(r);
umax=max(u);
sum1=u1_d(n1+n2+3,1)
sum2=u1_d(n1+n2+3,2)
Zs=u1_d(n1+n2+4,1)
sum=u1_d(n1+n2+4,2)
plot(r,u)
axis([0 rmax 0 umax*1.1])
text(rmax*0.45,umax*1.01,['Zs =',num2str(Zs)])
text(rmax*0.45,umax*0.96,['R1 =',num2str(r1)])
text(rmax*0.45,umax*0.91,['R2 =',num2str(r2)])
text(rmax*0.45,umax*0.85,['sum1=',num2str(sum1)])
text(rmax*0.45,umax*0.80,['sum2=',num2str(sum2)])
text(rmax*0.45,umax*0.75,['Total=',num2str(sum)])
title('u1+u2 (Total)')
set(gca, 'box','on')
xlabel('r (0 - R2)')
ylabel('integrand')

print circular2.ps



