Message Boards Message Boards

0
|
5135 Views
|
5 Replies
|
1 Total Likes
View groups...
Share
Share this post:

FindRoot not recognising replaced symbol in nested function

Archie Watts-Farmer, Carnot Ltd
Posted 4 years ago

Hi All,

I have the following example FindRoot function as part of a Nest function.

Clear[del]
del[x_] := Module[{y, dP, res, OP},
  res := {y} /. FindRoot[y^2 + 2 y - x, {y, 0}];
  OP := res[[1]];
  dP := OP + 1;
  Nest[(# + dP/#)/2 &, 1.5, 2]
  ]
del[name] /. name -> 3

Mathematica doesn't recognise the "name" symbol and returns the following error:

"FindRoot::nlnum: The function value {0. -1. name} is not a list of numbers with dimensions {1} at {y$21486} = {0.}.

ReplaceAll::reps: {FindRoot[y$21486^2+2 y$21486-name,{y$21486,0}]} is neither a list of replacement rules nor a valid dispatch table, and so cannot be used for replacing."

If you set Nest to apply zero times or replace it with a function that doesn't repeat the code works fine. It also works if you replace FindRoot with a simple assignment for symbol "OP". However, it doesn't work if you use Solve instead of FindRoot. Can anyone shine any light on this?

Many thanks for any help.

Attachments:
Replace.nb
POSTED BY: Archie Watts-Farmer
5 Replies
Sort By:
Rohit Namjoshi
Rohit Namjoshi
Posted 4 years ago

Crossposted here.
POSTED BY: Rohit Namjoshi

Hi,

there is a clue in the error message.

{FindRoot[y214862+2y21486-name,{y$21486,0}]}

FindRoot attempts to find a root numerically. As you see the error message says that it fails to do so, because there is the variable y (with the $21486, because of the localisation of Module), and the variable: name. That cannot work. Solve uses other methods and in particular symbolic solutions which makes it work in this case (it doesn't care that name is substituted later.

This is basically a combination of (i) when are the variables substituted by numbers and (ii) do you use symbolics or numerics, I would think.

Best wishes, Marco
POSTED BY: Marco Thiel

Many thanks for the comprehensive explanation Marco.
POSTED BY: Archie Watts-Farmer

Many thanks Marco, much appreciated. Interestingly, if you use solve with the expression equated to zero it also works (in my original post I'd tried Solve without equating the expression to zero). 

Clear[del]
del[x_] := Module[{y, dP, res, OP},
  res := {y} /. Solve[y^2 + 2 y - x == 0, y];
  OP := res[[1]];
  dP := OP + 1;
  Nest[(# + dP/#)/2 &, 1.5, 2]
  ]
del[name] /. name -> 3

Can't figure out why this works and FindRoot doesn't.
POSTED BY: Archie Watts-Farmer

Hi,

this works:

Clear[del]
del[x_] := 
 Module[{y, dP, res, OP}, 
  res := {y} /. FindRoot[y^2 + 2 y - x, {y, 0}];
  OP := res[[1]];
  dP := OP + 1;
  Nest[(# + dP/#)/2 &, 1.5, 2]]
del[name /. name -> 3]

Of course,

del[3]

is more concise.

Cheers,

Marco
POSTED BY: Marco Thiel
Reply to this discussion
Community posts can be styled and formatted using the Markdown syntax.
Tag limit exceeded
Note: Only the first five people you tag will receive an email notification; the other tagged names will appear as links to their profiles.
Reply Preview
Attachments
Remove
or Discard

Group Abstract Group Abstract