2-3.f90

来自「用fortran编写的双曲线方程的数值算法」· F90 代码 · 共 49 行

  !SOLUTION OF HYPERBOLIC EQN
  !EXPLICIT SCHEME
    DIMENSION F(100), FO(100)
	 
		PRINT*,'请输入时间的步长DT'
	READ (*,*) DT
	PRINT*,'请输入空间的步长DX'
	READ (*,*) DX
       PRINT*,'请输入C'
	READ (*,*) C

		PRINT*,'请输入空间序号最大值IMAX'
	READ (*,*) IMAX


    DO  I=1, IMAX
    FO(I)=0.
	ENDDO
   
    READ (*,*) (FO(I),I=1,4)

 
       DO I=1, IMAX
          F(I)=FO(I)
       ENDDO
  
    N=0
	T=0.
   DO WHILE(N<100)
       T=T+DT
       N=N+1
       DO  I=2, IMAX
       F(I)=FO(I)-C*DT/DX*(FO(I)-FO(I-1))
       ENDDO
     
	   DO  I=1, IMAX
       FO(I)=F(I)
   	   ENDDO
 	 
	   WRITE (6,11) (FO(I),I=1,IMAX)
	   
   	ENDDO
	
	WRITE (6,11) (FO(I),I=1,IMAX)
	   11 FORMAT (10F10.4)
	STOP
	END