用Matlab编程：试用4、5阶龙格-库塔法求下列微分方程组的数值解，并画出解的曲线图 dx/dt=yz dy/dt=-xz dz/dt=-0.5xy x（0）=0 y（0）=1 z（0）=1 写出完整的运行代码

要使用Matlab编程解决给定的微分方程组并绘制解的曲线图，我们首先需要定义微分方程组的函数，然后使用Matlab内置的数值解法函数来求解。这里，我们将使用4阶和5阶龙格-库塔法，分别对应于ode45和ode45函数（Matlab中ode45是4阶和5阶的自适应龙格-库塔法）。

以下是完整的Matlab代码，用于求解给定的微分方程组并绘制结果曲线图：

function main
    % 初始条件
    x0 = [0, 1, 1]; % [x(0), y(0), z(0)]
    tspan = [0 10]; % 定义时间范围

    % 使用ode45求解微分方程组
    [t, y] = ode45(@odefun, tspan, x0);

    % 绘制结果曲线图
    figure;
    subplot(3, 1, 1);
    plot(t, y(:,1));
    title('x(t)');
    xlabel('Time (t)');
    ylabel('x');

    subplot(3, 1, 2);
    plot(t, y(:,2));
    title('y(t)');
    xlabel('Time (t)');
    ylabel('y');

    subplot(3, 1, 3);
    plot(t, y(:,3));
    title('z(t)');
    xlabel('Time (t)');
    ylabel('z');
end

% 微分方程组的函数
function dydt = odefun(t, y)
    dydt = zeros(3, 1);
    dydt(1) = y(2) * y(3); % dx/dt = yz
    dydt(2) = -y(1) * y(3); % dy/dt = -xz
    dydt(3) = -0.5 * y(1) * y(2); % dz/dt = -0.5xy
end

在这段代码中，我们首先定义了一个名为main的函数，它设置了初始条件和时间范围。然后，我们调用ode45函数来求解微分方程组，并将结果存储在t和y中。ode45函数的第一个参数是一个函数句柄@odefun，它指向定义微分方程组的函数odefun。odefun函数接受时间t和状态向量y作为输入，并返回状态的变化率dydt。

最后，我们使用plot函数绘制了三个子图，分别表示x(t)、y(t)和z(t)随时间的变化。

请注意，这段代码是一个示例，实际使用时需要在Matlab环境中运行。1