三阶ODE四阶龙格-库塔法求解问题：洛伦兹方程MATLAB代码修正请求

问题详情 #

• 求解对象：洛伦兹方程
• 参数配置：σ=10，b=8/3，r=28
• 初始条件：x=y=z=5
• 积分区间：$t \in [0,20]$
• 求解方法：步长h=0.03125的四阶龙格-库塔（RK4）法
• 预期输出：得到x（对应变量y(1)）随t变化的混沌曲线
• 当前异常：计算结果收敛至0（yl(end,1)=-6.49e-48），生成错误曲线

涉及代码文件 #

需检查修正以下三个MATLAB代码文件：

• LorenzPlot.m
• rk4sys.m
• lorenz.m

请求修正上述代码，以得到符合预期的正确求解结果。

内容的提问来源于stack exchange，提问作者Turgut Ataseven