C First index is chip C Second index is test. dimension alpha(4,4),beta1(4,4),beta2(4,4) real nonlin dimension x(6),y(6) alpha(1,1)=1.1150 alpha(1,2)=1.1120 alpha(1,3)=1.1128 alpha(1,4)=1.1103 beta1(1,1)=-0.0338 beta1(1,2)=-0.0314 beta1(1,3)=-0.0367 beta1(1,4)=-0.0335 beta2(1,1)=+0.0057 beta2(1,2)=+0.0058 beta2(1,3)=+0.0074 beta2(1,4)=+0.0063 alpha(2,1)=1.0123 alpha(2,2)=1.0177 alpha(2,3)=1.0178 alpha(2,4)=1.0172 beta1(2,1)=-0.0123 beta1(2,2)=-0.0109 beta1(2,3)=-0.0150 beta1(2,4)=-0.0141 beta2(2,1)=0.0017 beta2(2,2)=0.0020 beta2(2,3)=0.0029 beta2(2,4)=0.0028 alpha(3,1)=1.0000 alpha(3,2)=1.0000 alpha(3,3)=1.0000 alpha(3,4)=1.0000 beta1(3,1)=-0.0099 beta1(3,2)=-0.0083 beta1(3,3)=-0.0124 beta1(3,4)=-0.0118 beta2(3,1)=+0.0014 beta2(3,2)=+0.0017 beta2(3,3)=+0.0027 beta2(3,4)=+0.0027 alpha(4,1)=1.0707 alpha(4,2)=1.0639 alpha(4,3)=1.0641 alpha(4,4)=1.0638 beta1(4,1)=-0.0163 beta1(4,2)=-0.0133 beta1(4,3)=-0.0178 beta1(4,4)=-0.0168 beta2(4,1)=+0.0027 beta2(4,2)=+0.0025 beta2(4,3)=+0.0037 beta2(4,4)=+0.0035 x(1)=100. x(2)=1000. x(3)=10000. x(4)=30000. x(5)=50000. x(6)=60000. do ichip=1,4 write(6,33) ichip 33 format(/'Chip:', i5) do itime=1,4 do i=1,6 xold=x(i) xf=xold/32000. nonlin=(1.+beta1(ichip,itime)*xf+beta2(ichip,itime)*xf**2) xnew=xold*alpha(ichip,itime)*nonlin y(i)=xnew Cwrite(6,1),itime,ichip,xold,nonlin,xnew/xold if (i.eq.6) write(6,*) itime, (y(6)/y(1))/(x(6)/x(1)) 1 format(' Test',i5,f10.2,f10.5,f10.5) end do end do end do stop end