限制plot窗格为0<y<100

要限制plot窗格为0<y<100，可以使用`ylim`函数，将y轴的取值范围限制在0到100之间。修改后的代码示例如下：

% 定义常数和步长
a = 0;
b = 3000;
h = 0.01;
N = (b-a)/h;

% 定义初始条件
u = zeros(N+1, 1);
v = zeros(N+1, 1);
u(1) = 0;
v(1) = 1;
f = @(x, u, v) [u(2); 1000*(1-u(1)^2)*u(2)-u(1)];

% 使用龙格库塔法进行数值求解
for i = 1:N
    k1 = h * f(a+(i-1)*h, [u(i);v(i)]);
    k2 = h * f(a+(i-1)*h+h/2, [u(i)+k1(1)/2;v(i)+k1(2)/2]);
    k3 = h * f(a+(i-1)*h+h/2, [u(i)+k2(1)/2;v(i)+k2(2)/2]);
    k4 = h * f(a+i*h, [u(i)+k3(1);v(i)+k3(2)]);
    u(i+1) = u(i) + 1/6*(k1(1)+2*k2(1)+2*k3(1)+k4(1));
    v(i+1) = v(i) + 1/6*(k1(2)+2*k2(2)+2*k3(2)+k4(2));
end

% 绘制结果图形，限制y轴范围为0到100
x = a:h:b;
plot(x, u);
title('龙格库塔法解微分方程u''=y'', v''=1000(1-u^2)v-u');
xlabel('x');
ylabel('u');
ylim([0, 100]);

这段代码将使用龙格库塔法对微分方程进行数值求解，并绘制出结果图形，同时限制了y轴范围在0到100之间。