function y=f(x)
a=1.42; b=0.3; 
s=size(x); k=s(2); y=0*x;
for j=1:k ,
 y(:,j)=[a-x(1,j)^2+b*x(2,j);x(1,j)]; end

