%
disp(' ');
disp(' arclength_sine.m  ver 1.1  May 29, 2007 ');
disp(' by Tom Irvine  Email: tomirvine@aol.com ');
%
disp(' This program integrates the arc length of a half-sine function ');
%
disp(' ');
A=input(' Enter amplitude ');
L=input(' Enter length ');
x=L; % integration limit
%
disp(' ');
Q = quad(@arc_function,0,x,[],[],A,L);
%
out4 = sprintf('     Length = %g ',L);
disp(out4);
%
out5 = sprintf(' Arc Length = %g ',Q);
disp(out5);
%
delta=Q-L;
%
out6 = sprintf('\n delta = %g ', delta);
disp(out6);
%
strain = (delta/L)*1.0e+06;
%
out7 = sprintf(' Strain = %g micro unit/unit',strain);
disp(out7);