function out=findcrit(func,vars,point,param)
   
  gradfunc=jacobian(func,vars);
hessfunc=jacobian(gradfunc,vars);
H=subs(hessfunc,vars,point);
v=subs(gradfunc,vars,point);
out=point-(param*v*inv(H))/norm(v*inv(H));