sub dif.for

来自「非线性回归问题SQP解法」· FOR 代码 · 共 11 行

	SUBROUTINE DIF(N,ME,MI,M,NAMAX,X,OGRA,CGRA,Y,CY1,CY2,NF)
	DOUBLE PRECISION X(N),Y(N),OGRA(N),CGRA(NAMAX,M),FY1,FY2,CY1(N),
     *		CY2(N),DELTA,STEP,DSTEP
	DATA DELTA/1.0D-5/
	DO 3 J=1,N
		DO 1 I=1,N
			Y(I)=X(I)
			CGRA(J,K)=(CY1(K)-CY2(K))/DSTEP
1		CONTINUE
3	CONTINUE
	END