求助四阶龙格库塔

fcode

咋说呢，我也没有头绪，发现了好多个问题。
因为不知道你的意图，也不知道怎么改。

例如：
273行开始
DO j = 1, M
        DO !//此处死循环
           IF(SX<X(1).OR.SX>X(N)) THEN
死循环了


再例如。
主程序中，p0(i)=i*1.d0
但p0并没有传入rungekutta龙格库塔子程序中。该子程序中的 p0 没有定义，也没有值。

fcode

你先把死循环和p0没有传入的问题修改了吧。

帮助你理解需求并完全纠正逻辑问题，需要太多的精力，我不一定能帮助你。

你可以把上述问题修改后再看看，有问题再继续帖新的代码求助。

另外，断点单步调试会对你有很大的帮助： http://debug.w.fcode.cn  和 http://v.fcode.cn/video-debug.html