Re: serious NDSolve bug?
- To: mathgroup at smc.vnet.net
- Subject: [mg44211] Re: serious NDSolve bug?
- From: Paul Abbott <paul at physics.uwa.edu.au>
- Date: Tue, 28 Oct 2003 05:53:18 -0500 (EST)
- Organization: The University of Western Australia
- References: <bmod0f$hof$1@smc.vnet.net> <bmqptd$ri6$1@smc.vnet.net>
- Sender: owner-wri-mathgroup at wolfram.com
In article <bmqptd$ri6$1 at smc.vnet.net>,
Jens-Peer Kuska <kuska at informatik.uni-leipzig.de> wrote:
> and whe you write
>
> f[kampas_?NumericQ]:=(u[\[Xi]] /. First[sprendinys1 =
> NDSolve[{Derivative[2][\[Psi]][\[Xi]] == (\[Xi]^2 - energija)*
> \[Psi][\[Xi]], Derivative[1][u][\[Xi]] == \[Psi][\[Xi]]^2,
> \[Psi][\[Xi]min] == 0,
> Derivative[1][\[Psi]][\[Xi]min] == kampas, u[\[Xi]min] == 0},
> {\[Psi][\[Xi]], u[\[Xi]]},
> {\[Xi], \[Xi]min, \[Xi]max}]] /. \[Xi] -> \[Xi]max)
>
> and call
>
> teisingasKampas =
> FindRoot[f[kampas]==1, {kampas, -0.04, 0.04}]
>
> all works as expected *and* your code is better readable.
>
> BTW: Can you avoid postings with Mathematica's special characters ?
Why? I thought one of the major features of Mathematica is that we can
(and should) use special characters, and that Notebooks are portable
across all platforms because everything has an underlying ASCII
representation. Of course, the above ASCII is hard to read but when you
paste it into a Notebook it becomes perfectly readable.
Of course, it would be much if we could simply post Notebooks to
MathGroup ...
Cheers,
Paul
--
Paul Abbott Phone: +61 8 9380 2734
School of Physics, M013 Fax: +61 8 9380 1014
The University of Western Australia (CRICOS Provider No 00126G)
35 Stirling Highway
Crawley WA 6009 mailto:paul at physics.uwa.edu.au
AUSTRALIA http://physics.uwa.edu.au/~paul