% surface plot of rosenbrock function
rosen = @(x,y)100*(y-x.^2).^2+(1-x).^2;
ezsurf(rosen)
ezsurf(rosen,[0,2,0,2])
rosen(1,1)