二阶泰勒展开法求解微分方程

最新推荐文章于 2024-05-16 10:30:00 发布

老李今天学习了吗

最新推荐文章于 2024-05-16 10:30:00 发布

阅读量4.1k

点赞数

分类专栏：偏微分方程数值解文章标签：算法

本文链接：https://blog.csdn.net/weixin_29732003/article/details/105185497

版权

偏微分方程数值解专栏收录该内容

7 篇文章 4 订阅

订阅专栏

日期：3030-3-29

二阶泰勒展开法求解微分方程

在这里插入图片描述

代码示例

%% 1
h = 0.5; t = 2:h:3;
N = length(t);
y = zeros(1,N); y(1) = 1;
%求fi
fi = @(x,y,t)(1+(t-y).^2-h*((t-y).^3));
%求y(3)
for i =1:N-1
    y(i+1) = y(i)+h*fi(y(i),t(i),h);
end
y(N)
%% 2
h = 0.25; t = 0:h:1;
N = length(t);
y = zeros(1,N); y(1) = 1;
%求fi
fi = @(x,y,t)(cos(2*t)+sin(3*t)+(h/2)*(-2*sin(2*t)+3*cos(3*t)));
%求y(3)
for i =1:N-1
    y(i+1) = y(i)+h*fi(y(i),t(i),h);
end
y(N)