2楼: Originally posted by Yan_Jordan at 2017-06-23 16:57:50
首先 程序不全，没有f函数， 建议贴出错误信息。

从以上代码看，楼主对x y两个数组的定义（没有维度信息）很有问题，需要提前初始化好数组，维度应该是N+1（如果N不定，可以用动态数组）

我尝试猜测了下楼主的 ...

PROGRAM test
        IMPLICIT NONE ! 强制所有变量都有初始化定义

        real :: h,yc
        INTEGER :: N, i
        real, dimension(:), ALLOCATABLE :: x, y !定义动态数组
        h=0.1
        N=floor(2/H)
       
        ALLOCATE(x(N+1), y(N+1)) !分配动态数组内存大小
        x(1)=0
        y(1)=-1  !迭代的初值
       
        do i=0,N
                x(i+1)=x(i)+h
                y(i+1)=y(i)+h*f(x(i),y(i))
                yc=y(i)+h*f(x(i+1),y(i+1))
                y(i+1)=(y(i+1)+yc)/2
                print*,y(i+1)
        enddo

        DEALLOCATE(x,y) !释放动态数组内存
        
   CONTAINS
        
        real FUNCTION f(x,y)
                IMPLICIT NONE
                real, INTENT(in) :: x, y
                f=sqrt(x**2+y**2)
        ENDFUNCTION

END PROGRAM test

2楼: Originally posted by Yan_Jordan at 2017-06-23 16:57:50
首先 程序不全，没有f函数， 建议贴出错误信息。

从以上代码看，楼主对x y两个数组的定义（没有维度信息）很有问题，需要提前初始化好数组，维度应该是N+1（如果N不定，可以用动态数组）

我尝试猜测了下楼主的 ...

PROGRAM test IMPLICIT NONE ! 强制所有变量都有初始化定义 real :: h,yc INTEGER :: N, i real, dimension(:), ALLOCATABLE :: x, y !定义动态数组 h=0.1 N=floor(2/H) ALLOCATE(x(N+1), y(N+1)) !分配动态数组内存大小 x(1)=0 y(1)=-1 !迭代的初值 do i=0,N x(i+1)=x(i)+h y(i+1)=y(i)+h*f(x(i),y(i)) yc=y(i)+h*f(x(i+1),y(i+1)) y(i+1)=(y(i+1)+yc)/2 print*,y(i+1) enddo DEALLOCATE(x,y) !释放动态数组内存 CONTAINS real FUNCTION f(x,y) IMPLICIT NONE real, INTENT(in) :: x, y f=sqrt(x**2+y**2) ENDFUNCTION END PROGRAM test