function f = foscilw(y,t)
ome=1.;
f(1)=y(2);
if (y(1)>0)
 f(2)=-ome*ome*y(1);
else
 f(2)=-100*ome*ome*y(1);
endif
end