www.pudn.com > Fortran.zip > JJSIM0.FOR, change:1994-06-10,size:672b

 
	EXTERNAL FS 
	DIMENSION X(2),XX(2,3),F(3),XT(2),XF(2) 
	DOUBLE PRECISION X,Z,XX,F,FS,XT,XF 
	N=2 
	M=3 
	D=1.0 
	U=1.6 
	V=0.4 
	EPS=1.0E-30 
	CALL JJSIM(N,M,D,U,V,X,Z,EPS,K,FS,XX,F,XT,XF) 
	WRITE(*,*) 
	WRITE(*,10) K 
10	FORMAT(1X,'K=',I4) 
	WRITE(*,*) 
	WRITE(*,20) 
20	FORMAT(7X,'X(1)',11X,'X(2)',11X,'F') 
	DO 30 I=1,M 
30	WRITE(*,40) XX(1,I),XX(2,I),F(I) 
40	FORMAT(1X,3D15.6) 
	WRITE(*,*) 
	WRITE(*,50) (I,X(I),I=1,N) 
50	FORMAT(1X,'X(',I2,' )=',D15.6) 
	WRITE(*,60) Z 
60	FORMAT(1X,'Z=',D15.6) 
	WRITE(*,*) 
	END 
 
	FUNCTION FS(N,X) 
	DIMENSION X(N) 
	DOUBLE PRECISION X,S,FS 
	S=X(2)-X(1)*X(1) 
	S=100.0*S*S 
	FS=S+(1.0-X(1))*(1.0-X(1)) 
	RETURN 
	END