数值分析 四阶龙格库塔法

【数值分析 四阶龙格库塔法】四阶常微分龙格-库塔法求解常微分方程初值问题 。四阶龙格 库塔法的局部截断误差阶为O(h5) , 欧拉法的绝对稳定实域为2 MATLAB四阶龙格库塔法解微分方程;Functiondytest(t,clc%清除内存中的变量%数值yinline( x * exp(x) 2 * x1 );y(1)%四阶龙格库塔法y0 clear,CLC%清除内存中的变量%数值solves yinline( x * exp(x) 2 * x1 );y(1)%四阶龙格库塔法y0 , y0);y1(end,1)%y(1) functiondyfun(x,y)dyzeros(3,1)的值;dy(1)y(2);dy(2)y(3);dy(3)y(3) y(2)y(1) 2 * x3;运行结果ans3.7183ans3.718 。
132];[x1,y1]ode45(@fun,[0#include#include/*n表示几个部分,n 1表示他输出的个数*/intrungkutta (doubley0,doublea,doubleb,intn,double * x,double * y,intstyle,double (* function) (double,double)) {doubleh
k2,k3,k4;inti//x(double *)malloc((n 1)* sizeof(double));//y(double *)malloc((n 1)* sizeof(double));x[0]a;y[0]y0;switch(style){ case 2:for(i0; 。

    推荐阅读