视频中的Matlab代码为：

%%
clear; clc;
dsolve('Du = 1+u^2', 't')
syms u(t)
dsolve(diff(u,1) == 1+u^2)

%%
clear;clc;
[t,y] = ode23(@(t,y) 2*t, [0 5], 0);
plot(t,y, '-o')