hdfp.for

优化算法,是FORTRUN语言写。都经过了运行。

	SUBROUTINE HDFP(H,N,DELG,DELX,IE)
	DIMENSION H(20,20),DELX(20),DELG(20),DHG(20),DGH(20)
	DXDG=0.0
	DGHDG=0.0
	DO 20 I=1,N
	DHG(I)=0.0
	DGH(I)=0.0
	DO 10 J=1,N
	DHG(I)=DHG(I)-H(I,J)*DELG(J)
10	DGH(I)=DGH(I)+DELG(J)*H(J,I)
	DXDG=DXDG+DELX(I)*DELG(I)
20	DGHDG=DGHDG+DGH(I)*DELG(I)
	IE=-1
	IF(ABS(DXDG).LT.0.1E-60) GO TO 40
	IF(DGHDG.LT.0.1E-60) GO TO 40
	IE=0
	DO 30 I=1,N
	DO 30 J=1,N
30	H(I,J)=H(I,J)+DELX(I)*DELX(J)/DXDG+DHG(I)*DGH(J)/DGHDG
40	RETURN
	END