2-1.f90

用fortran编写的抛物线方程的数值算法

    DIMENSION U(100),UO(100)
	    PRINT*,'请输入空间的步长DX'
	READ (*,*) DX
		PRINT*,'请输入时间的步长DT'
	READ (*,*) DT
		PRINT*,'请输入空间序号最大值IMAX'
	READ (*,*) IMAX

    IMAX1=IMAX-1
	DO  I=1,IMAX
	U(I)=0.
  	UO(I)=0.
	ENDDO
    T=0.
	N=0
  

    100   	T=T+DT
            N=N+1
	     IF (T.LT..1) THEN
		   U(IMAX)=T*100
		 ELSE
	       U(IMAX)=10.
	     ENDIF
        
            DO  I=2,IMAX1
           	U(I)=UO(I)+DT*1./DX**2*(UO(I+1)+UO(I-1)-2.*UO(I))
	        ENDDO
             DO I=1,IMAX
  	         UO(I)=U(I)
             ENDDO
 

	
	WRITE (6,5) (UO(I),I=1,IMAX)
	  5 FORMAT(10 F10.4)
	  IF (N.LT.10) GO TO  100
	STOP
	END