A=[ln37.225 1 ln6 ln30 ln4000

ln122.15 1 ln6 ln30 ln4500

ln160.475 1 ln6 ln30 ln5000

ln150.725 1 ln6 ln30 ln5500

ln184.375 1 ln6 ln30 ln6000

ln166.175 1 ln6 ln44 ln4000

ln135.225 1 ln6 ln44 ln4500

ln52.075 1 ln6 ln44 ln5000

ln90.575 1 ln6 ln44 ln5500

ln127.75 1 ln6 ln44 ln6000

ln129.35 1 ln6 ln66 ln4000

ln132.525 1 ln6 ln66 ln4500

ln145.125 1 ln6 ln66 ln5000

ln145.475 1 ln6 ln66 ln5500

ln141.075 1 ln6 ln66 ln6000

ln130.625 1 ln6 ln92 ln4000

ln152.875 1 ln6 ln92 ln4500

ln157.1 1 ln6 ln92 ln5000

ln130.65 1 ln6 ln92 ln5500

ln161.025 1 ln6 ln92 ln6000

ln152.55 1 ln6 ln120 ln4000

ln114.85 1 ln6 ln120 ln4500

ln158.4 1 ln6 ln120 ln5000

ln132.7 1 ln6 ln120 ln5500

ln120.325 1 ln6 ln120 ln6000

ln36.2 1 ln4 ln30 ln4000

ln55.175 1 ln4 ln44 ln4000

ln70.6 1 ln4 ln66 ln4000

ln32.53 1 ln4 ln92 ln4000

ln79.2 1 ln4 ln120 ln4000

76.775 1 ln4 ln30 ln6000

ln103.325 1 ln4 ln44 ln6000

ln99.05   ln4 ln66 ln6000

ln99.7 1 ln4 ln92 ln6000

y=A(:,1),x=A(:,2:5)

B=REGRESS(y, x, alpha)

e=B(1),f=B(2),g=B(3),h=B(4)

[b,bint,r,rint,stats]=regress(y, x)

b,bint,stats