www.pudn.com > Fortran.zip > KFOUR0.FOR, change:1995-03-31,size:381b

 
 
	DIMENSION F(101), A(51), B(51) 
	DOUBLE PRECISION F,A,B,C,H 
	WRITE(*,*) 
	WRITE(*,10) 
10	FORMAT(3X,'N',12X,'A(N)',16X,'B(N)') 
	H=6.283185306/101.0 
	DO 20 I=1,101 
	  C=H*(I-0.5) 
	  F(I)=C*C 
20	CONTINUE 
	N1=51 
	N2=101 
	CALL KFOUR(F,N1,N2,A,B) 
	DO 30 I=1,N1 
	  K=I-1 
	  WRITE(*,40) K,A(I),B(I) 
30	CONTINUE 
	WRITE(*,*) 
40	FORMAT(1X,I3,5X,D15.6,4X,D15.6) 
	END