y''-1000(1-y^2)*y'+y=0;y(0)=0,y'(0)=1

这是一个二阶常微分方程，可以通过变量代换将其化为一阶方程组的形式。令 $z_1=y$，$z_2=y'$，则有：

$$
\begin{cases}
z_1' = z_2 \\
z_2' = 1000(1-z_1^2)z_2-z_1
\end{cases}
$$

这是一个一阶微分方程组，可以使用常规的数值方法求解。比如，使用Matlab的ode45函数可以求解该微分方程组，代码如下：

% 定义微分方程组
function dydt = myODE(t,y)
    dydt = [y(2); 1000*(1-y(1)^2)*y(2)-y(1)];
end

% 指定初始条件
y0 = [0; 1];
tspan = [0 1];

% 调用函数求解微分方程组
[t,y] = ode45(@myODE,tspan,y0);

% 绘制解的图像
plot(t,y(:,1))
xlabel('t')
ylabel('y')

这段代码使用ode45函数求解微分方程组，并绘制出 $y$ 随时间变化的图像。运行后，可以得到该微分方程的解。