disp('   stress_intensity.m   ver 1.0 ')
disp(' ')
disp(' by Tom Irvine  Email: tomirvine@aol.com ')
disp(' ')
disp(' Stress Corrosion Cracking (SCC)')
disp(' ')
disp(' Stress required to propagate a partial through-the-thickness')
disp(' crack-like discontinuity ')
disp(' ')
disp(' Reference:  ME File: 044-1118 Rev -')
disp(' ')
%
K=input(' Enter stress intensity K IC  (ksi sqrt(in))  ')
a=input(' Enter depth  a (in)  ')
t=input(' Enter thickness  t (in)  ')
%
arg=(pi*a)/(2.*t);
%
den=sqrt(pi*a)*sqrt(sec(arg));
stress=K/den;
%
out5 = sprintf('\n Stress = %10.3g ksi\n',stress);
disp(out5);