谢谢！没系统学过fortran，不知到interface还能这么用……真的谢谢啦

[Fortran] 查看源码 复制源码

program ode
use rksuite_90
implicit none
……
external::f
dimension(2)::f
……
call range_integrate_r1(comm,f,t_want,t_got,y_got,yderiv_got,flag)
stop
end program ode

function f(t,y)
use rksuite_90
implicit none
real(kind=wp)::t
real(kind=wp)::y(2)
real(kind=wp)::f(2)     !声名返回值类型
f(1)=y(2)
f(2)=-y(1)
end function