MathGroup Archive 1996

[Date Index] [Thread Index] [Author Index]

Search the Archive

Re: help on critical points


>From: Jose Luis Lugo Goytia <lugo at hplara.iquimica.unam.mx>
>Subject: [mg4037] help on critical points

>Could you please give me some advise?

Jose: one small change is needed to your code.

   gradiente[f_,var_]:=Outer[D,{f},var];
   ceros[f_,vars_]:=
      vars /.
         Solve[
             Numerator[Together[f]] == Table[0,{i,1,Length[vars]}],
             vars
         ];
(*
Change
   PuntosCriticos[f_,vars_]:=
      ceros[gradiente[f @@ vars,vars][[1]],vars] //N;
   to
*)

   PuntosCriticos[f_,vars_]:= ceros[gradiente[f,vars][[1]],vars] //N;

(*test*)
   PuntosCriticos[x^2 y - y^2 -2 x y,{x,y}]

      {{1., -0.5}, {0, 0}, {2., 0}}

Allan Hayes
hay at haystack.demon.co.uk

==== [MESSAGE SEPARATOR] ====


  • Prev by Date: Re: help on critical points
  • Next by Date: Please Help
  • Previous by thread: Re: help on critical points
  • Next by thread: Please Help