MathGroup Archive 2008

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

Search the Archive

Re: Re: Two FindRoot questions

  • To: mathgroup at smc.vnet.net
  • Subject: [mg89888] Re: Re: [mg89872] Two FindRoot questions
  • From: Bob Hanlon <hanlonr at cox.net>
  • Date: Tue, 24 Jun 2008 03:21:16 -0400 (EDT)
  • Reply-to: hanlonr at cox.net

Or you could put the SetPrecision inside the function definition

Clear[f];

f[k_ /; k <= 0] := 
  x /. FindRoot[x^2 + SetPrecision[k, 20], {x, 0, 10}, 
    WorkingPrecision -> 20];

{f[-10], f[10]}

{3.1622776601683793320,f(10)}

Plot[f[k], {k, -10, 10}]


Bob Hanlon

---- Bob Hanlon <hanlonr at cox.net> wrote: 
> Clear[f];
> 
> f[k_ /; k <= 0] := 
>   x /. FindRoot[x^2 + k, {x, 0, 10}, WorkingPrecision -> 20];
> 
> {f[-10], f[10]}
> 
> {3.1622776601683793320,f(10)}
> 
> Plot[f[SetPrecision[k, 20]], {k, -10, 10}]
> 
> 
> Bob Hanlon
> 
> ---- Aaron Fude <aaronfude at gmail.com> wrote: 
> > Hi,
> > 
> > These are not FindRoot questions, per se...
> > 
> > Here's a simple example which I want to ask three questions about:
> > 
> > f[k_] := x /. FindRoot[x^2 + k, {x, 0, 10}];
> > f[-10]
> > Plot[f[k], {k, -10, 10}]
> > 
> > First, I want the plot to only show where there exists a root.
> > Is the right solution to make f[] return Null?
> > How do I make f[] return Null? (Is there a way to "catch" the
> > warnings?)
> > 
> > Finally, I need to solve my equations to 20 digits. How do I do that?
> > I've read about Accuracy and Precision but it didn't help.
> > 
> > Thanks!
> > 



  • Prev by Date: Re: Export a matrix in C form
  • Next by Date: Re: Solving a Sum
  • Previous by thread: Re: Two FindRoot questions
  • Next by thread: Re: Two FindRoot questions