IF(ABS(B-BO).GT.BMIN) GO TO 500 



GO TO 900 

 C 



C CORRECTION TERM NOT FINE ENOUGH, START NEW ITERATION 

 C 



500 ITERAT=ITERAT+1 



AO=A 



BO=B 

 POLD=PNEW 



IF(ILIST.EQ.1)WRITE(6,520)ITERAT,NBIS,A0,B0,P0LD 

 520 FORMATC ' ,I6,10X,I4,6X,3(4X,F10.4)) 



NBIS=0 



GO TO 150 

 900 ITERAT=ITERAT+1 



IF(ILIST.EQ.1)WRITE(6,920)ITERAT,NBIS,A0,B0,P0LD 

 920 FORMATC ' ,I6,10X,I4,6X,3(4X,F10.4)) 



RETURN 



END 

 C 

 C 



C FUNCTIONS TO CALCULATE POWER LAW FUNCTION AND PARTIAL 

 C DERIVATIVES WITH RESPECT TO A AND B 



C 



FUNCTION F1(X3,A3,B3) 

 C CALCULATE PARTIAL OF A*X**B WITH RESPECT TO A 

 F1=X3**B3 



A3=A3 



RETURN 



END 

 C 



FUNCTION F2(X4,A4,B4) 

 C CALCULATE PARTIAL OF A*X**B WITH RESPECT TO B 



F2=(X4**B4)*AL0G(X4) 



A4=A4 



RETURN 



END 

 C 



FUNCTION F3(X5,A5,B5) 

 C CALCULATE POWER LAW A*X**B 



F3=A5*X5**B5 



RETURN 



END 



134 



