function dx=vanderpol(x,t)
dx=zeros(x);
dx(1) = x(2);
dx(2) = ((1-x(1)^2)*x(2)-x(1))/1.0d-6;
endfunction
