本例求解DDE方程y' = ddex1de(t,y,Z)使用dde23，然后绘制出解决方案。

用dde23．

Sol = dde23(@ddex1de， [1 0.2]， @ddex1hist， [0 5]);

在区间内的100点处评估解[0 5]．

X = linspace(0,5);Y = deval(sol,x);

画出解。

情节(x, y)

这个例子解决了系统的问题Y ' = vdp1(t, Y)使用数值，然后绘制出解的第一个组成部分。

用数值．

Sol = ode45(@vdp1， [0 20]， [2 0]);

在区间内的100个点处评估解的第一个分量20 [0]．

X = linspace(0,20,100);Y = deval(sol,x,1);

画出解。

情节(x, y)

求解简单ODEY ' = t^2在初始条件下Y0 = 0在这段时间 $$［0，3.］$$使用ode23．

Sol = ode23(@(t,y) t^2， [0 3]， 0);

从7个方面评估解决方案。解决方案结构索尔包含一个插值函数德瓦尔用于在这些点上产生连续溶液。指定第二个输出参数德瓦尔还要返回插值函数在指定点处的导数。

X = linspace(0,3,7);[y,yp] = deval(sol,x)

y =1×70 0.0417 0.3333 1.1250 2.6667 5.2083 9.0000

yp =1×70 0.2500 1.0000 2.2500 4.0000 6.2500 9.0000