写一个求解常微分方程的matlab代码

下面是一个简单的 Matlab 代码，可以用于求解常微分方程：

function y = ode_solver(f, tspan, y0)
% f: 函数句柄，表示 y' = f(t, y)
% tspan: 时间跨度，如 [0, 10]
% y0: 初始条件，如 1
% y: 返回求解的结果

opts = odeset('RelTol',1e-6,'AbsTol',1e-8);  % 设置求解选项

[t, y] = ode45(f, tspan, y0, opts);  % 使用 ode45 求解常微分方程

plot(t, y);  % 绘制解曲线
xlabel('t');
ylabel('y');

使用示例：

% 定义常微分方程
f = @(t, y) -y + t^2 + 1;

% 设置时间跨度和初始条件
tspan = [0, 2];
y0 = 0;

% 求解常微分方程
y = ode_solver(f, tspan, y0);

上述代码中，`ode45` 是 Matlab 自带的一个求解常微分方程的函数，可以根据需要选择其他的求解函数，如 `ode23`, `ode113` 等。`opts` 是求解选项，可以根据实际需要进行设置。