a=60
h=15
f=@(x) 2*sqrt(1+4*h^2*x.^2/a^4)
quadl(f,0,a)
disp('press any key to continue')
pause
clc
eps=8.85e-12;
z=0.05
R=0.06
sig=300e-6;
field=sig*z/4/eps*quadl(@(r) 2*r./(z.^2+r.^2).^(3/2),0,R)
disp('press any key to continue')
pause
clc
x=quadl(@(s) -13.6*cos(53*s).*sin(8*(0.119-s)), 0, 0.119)