function sol=myfun(sol)
x=sol(:,end-1);
z =  3*(1-x(1)).^2.*exp(-(x(1).^2) - (x(2)+1).^2) ... 
     - 10*(x(1)/5 - x(1).^3 - x(2).^5).*exp(-x(1).^2-x(2).^2)- 1/3*exp(-(x(1)+1).^2 - x(2).^2);
 sol(:,end)=z;