matlab常微分方程的差分方程如何代码实现

数舌哆猢筢学的很多计算可以利用matlab来做。对于常微分方程 : u'(t)=f(t,u) u(0)=u0。如何利用matlab计算函数u在一些节点上的数值？以向后euler法为例：u[n+1]=u[n]+h*f(t[n+1],u[n+1])

工具/原料

matlab

分析

1、向后euler法由于等式前后都有u[n+1],所以，不能向向前euler法一样，只需要做出循环：u[n+1]=u[n]+h*f(t[n],u[n])即可。

2、迭代部分：temp是迭代后的值，用于比较差是否小于误差精度。