function out=findmax(func,vars,point,param)
   
  gradfunc=jacobian(func,vars);
v=subs(gradfunc,vars,point);
out=point+(param*v)/norm(v);