function yh=shfun(v0)
  vx=v0(1); vy=v0(2);
  y10=-60+13*vx+2*vy-vx^2;
  y14=-272+7*vx+31*vy-vy^2;
  yh=[y10;y14];
end
