function m=prof(xs,ys)
for k=length(ys):-1:1
    m=[m;2-((xs-1).^2+4*(ys(k)-1)^2+2*xs*ys(k))];
end
