用matlab自行编写脚本，使用Runge-Kutta法求解微分方程，并分析误差。 完成函数y = RungeKutta(x,h)，函数参数含义同上，例如调用函数RungeKutta(3,0.25)即可计算出h=0.25参数下的y(3)。

function y = RungeKutta(x,h)
% x为要求解的自变量值，h为步长
% 定义微分方程dy/dx=f(x,y)，这里取f(x,y)=x+y
f = @(x,y) x + y;
% 初始值y0=1
y0 = 1;
% 计算迭代次数
n = round((x - 0) / h);
% 初始化y
y = y0;
% 迭代求解
for i = 1:n
k1 = f(0 + (i - 1) * h, y);
k2 = f(0 + (i - 1/2) * h, y + h/2 * k1);
k3 = f(0 + (i - 1/2) * h, y + h/2 * k2);
k4 = f(0 + i * h, y + h * k3);
y = y + h/6 * (k1 + 2*k2 + 2*k3 + k4);
end
% 计算精确解并输出误差
y_exact = exp(x) - x - 1;
error = abs(y_exact - y);
fprintf('当x=%.2f时，使用h=%.2f的Runge-Kutta法求解得到y=%.6f，\n',x,h,y);
fprintf('精确解为y=%.6f，误差为%.6f。\n',y_exact,error);
end

运行示例：

>> RungeKutta(3,0.25)
当x=3.00时，使用h=0.25的Runge-Kutta法求解得到y=17.342759，
精确解为y=17.367879，误差为0.025120。

相关问题

如何在MATLAB中利用4阶Runge-Kutta方法求解给定的常微分方程并展示完整的源码和操作步骤？

在MATLAB中实现4阶Runge-Kutta方法来求解常微分方程，关键在于准确地将算法的每一步转换为代码。以下是一个具体示例和步骤，帮助你掌握整个实现过程。

参考资源链接：[MATLAB实现4阶Runge-Kutta法求解ODE](https://wenku.csdn.net/doc/4uatihzf5e?spm=1055.2569.3001.10343)

首先，确保你有《MATLAB实现4阶Runge-Kutta法求解ODE》这一资源，它包含了详细的代码和解释，将帮助你理解算法的每个细节。

步骤1：定义微分方程
假设我们要求解的微分方程是 dy/dx = f(x, y)，例如 y' = -2xy，初始条件 y(0) = 1。我们首先定义函数 f(x, y)。

function dydx = f(x, y)
    dydx = -2 * x * y;
end

步骤2：编写4阶Runge-Kutta方法的实现代码
创建一个名为 `runge_kutta_4.m` 的文件，并编写以下代码：

function [x, y] = runge_kutta_4(f, x0, y0, x_end, n)
    % f - 微分方程函数句柄
    % x0 - 初始x值
    % y0 - 初始y值
    % x_end - 结束x值
    % n - 步数
    h = (x_end - x0) / n; % 步长
    x = x0:h:x_end; % x值数组
    y = zeros(1, length(x)); % 初始化y值数组
    y(1) = y0; % 初始条件
    for i = 1:n
        k1 = h * f(x(i), y(i));
        k2 = h * f(x(i) + 0.5 * h, y(i) + 0.5 * k1);
        k3 = h * f(x(i) + 0.5 * h, y(i) + 0.5 * k2);
        k4 = h * f(x(i) + h, y(i) + k3);
        y(i+1) = y(i) + (k1 + 2 * k2 + 2 * k3 + k4) / 6;
    end
end

步骤3：调用函数并绘制结果
在MATLAB命令窗口或脚本文件中调用 `runge_kutta_4` 函数：

[x, y] = runge_kutta_4(@f, 0, 1, 2, 100);
plot(x, y);
xlabel('x');
ylabel('y');
title('Runge-Kutta Method');

执行上述代码后，你将得到在区间 [0,2] 上使用100个步长计算的数值解，并将结果绘制成图表。

至此，你就完成了使用MATLAB的4阶Runge-Kutta方法来求解常微分方程的全过程。如果你希望深入理解算法背后的数学原理，或想进一步探索其它数值分析方法，可以深入阅读《MATLAB实现4阶Runge-Kutta法求解ODE》这一资料，它提供了全面的理论背景和实践指导。

参考资源链接：[MATLAB实现4阶Runge-Kutta法求解ODE](https://wenku.csdn.net/doc/4uatihzf5e?spm=1055.2569.3001.10343)