RE: FindRoot problem

• To: mathgroup at smc.vnet.net
• Subject: [mg40038] RE: [mg40032] FindRoot problem
• From: "David Park" <djmp at earthlink.net>
• Date: Mon, 17 Mar 2003 03:33:09 -0500 (EST)
• Sender: owner-wri-mathgroup at wolfram.com

```Robert,

Try using Ted Ersek's RootSearch package from MathSource. It even gets both
roots for you.

Needs["Enhancements`RootSearch`"]

f[x_ /; x >= 0 && x <= 1.54] :=
Module[{}, Sec[x] - 0.0018167 (Sec[x] - 1) - 0.002875 (Sec[x] - 1)^2 -
0.000808 (Sec[x] - 1)^3 ];

RootSearch[f[x] == 10, {x, 0, 1.54}]
{{x -> 1.4811}, {x -> 1.53501}}

David Park

From: Pigeon, Robert [mailto:Robert.Pigeon at drdc-rddc.gc.ca]
To: mathgroup at smc.vnet.net

Good day,
I have a problem with FindRoot that I do not understand, so I am
asking for help!  I have a function defined as a Module.  I use FindRoot to
find a root when the function is given a value.  Everything works fine
except when I put a condition to my function, when I do I get an error
message.  Here what I try to do:

f[x_]:= Module[{}, Sec[x] - 0.0018167 (Sec[x]-1) - 0.002875 (Sec[x]-1)^2 -
0.000808 (Sec[x]-1)^3 ];
f::usage="explanation of f";

If I do
FindRoot[f[x]==10,{x,1.4}]
I get
{x->1.53501}
All that is correct.

But now if I do
f[x_ /; x>= 0 && x<= 1.54]:= same as above
and I try to find the root the same way, I get an error message
FindRoot[f[x]==10,{x,1.4}]
I get
FindRoot::frjc: Could not symbolically find the Jacobian of {f[x]-10.  Try
giving two starting values for each variable.

I tried putting my conditions inside the Module, after the Module.  I tried
to use Which, If.  But, as soon as I put a condition, I get the same error.

Any idea ???

Thanks,

Robert

```

• Prev by Date: Re: FindRoot problem
• Next by Date: Re: Re: How to apply one plot range to another plot?
• Previous by thread: RE: FindRoot problem
• Next by thread: Re: FindRoot problem