怎么用四阶龙格库塔方法求解二阶微分方程

这个我用ode45的方法求出来了

但是不知道怎么用四阶龙格库塔方法求解，求指教

参考：Matlab中龙格-库塔(Runge-Kutta)方法原理及实现https://blog.csdn.net/xiaokun19870825/article/details/78763739?utm_medium=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-1.control&dist_request_id=&depth_1-utm_source=distribute.pc_relevant.none-task-blog-BlogCommendFromMachineLearnPai2-1.control