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