% Find the longitudinal young's modulus of a composite material 
datafile='c:/usr/stressstrain.txt'
fh=fopen(datafile,'r')
nval=fgets(fh)
n=eval(nval)

for i=1:1:n
    a=fgets(fh);
    b=sscanf(a,'%f  %f');
% Dividing by 100 as strain is giving as %
    strain(i)=b(1)/100;
% Multiplying stress by 10^6 as stress is given in MPa
    stress(i)=b(2)*10^6;
end
    
numer=0
for i=1:1:n
    numer=numer+stress(i)*strain(i);
end

denom=0
for i=1:1:n
    denom=denom+strain(i)^2;
end
  
Ym=numer/denom

fprintf('Youngs modulus in GPa %f',Ym/10^9)