function z=pr08_fce(p)
x=p(1);
y=p(2);
z=-(x.^2+x+y+y.^2./(1+x.^2)).*exp(-(x.^2+y.^2)./25);