x0 = 5;
xp0 = 0;
c = 1 ; 
K = 2;
sim('vanderpol');
close all;
plot(t,x);
comet(x,xp)