Re: help on critical points
- To: mathgroup at smc.vnet.net
- Subject: [mg4095] Re: [mg4037] help on critical points
- From: Allan Hayes <hay at haystack.demon.co.uk>
- Date: Tue, 4 Jun 1996 02:17:27 -0400
- Sender: owner-wri-mathgroup at wolfram.com
>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] ====