List Archives >
Maple User Group List Archive >
Archive by date >
This Month By Date >
This Month By Topic
[MUG] ODE with singularity
| [MUG] ODE with singularity |
|
Author: Jens-Uwe Herrmann
Posted: Mon, 18 Nov 2002 13:16:41 +0100
|
>> From: Jens-Uwe Herrmann "jens.herrmann"
hallo,
my problem is the following:
i have an ordinary diff. equ.
>ode256:=(1+3*alpha)/6*diff(diff(xi(eta),eta),eta)/(1+diff(xi(eta),eta)^2)^(3/2)+1+alpha-eta=0;
with alpha=[0,infinity], but i'm only interested in alpha=[0,3] anyway.
i.e. i say alpha:=0.51
i want to solve the ode with the following init. cond.
> ic256:=xi(1)=0,D(xi)(1)=-100000:
(-100000 is an aproximation for -infinity)
and i do it numeric. what comes out is fine and ok
> p:=dsolve({ode256,ic256},numeric,range=0..1);
Warning, cannot evaluate the solution further left of .11473470,
probably a singularity
p := proc(rkf45_x) ... end proc
now,there is a warning and i'm interested in this number. is there a way
to assign the singularity ( here eta=0.11473470 ) to some placeholder?
for example Sing:=0.11473470, without doing it by hand. something like
grep it.
thank you for listening...
jens-uwe
|
[View Complete Thread]
Previous by date: [MUG] Fourier Series, Kronecker-Delta ?, Andreas Wolf
Next by date: [MUG] Re: Too many levels of recursion, Maple User Group
Previous thread: [MUG] implicitdiff on generic polynoms, Charles James Leonardo Quarra Cappiello
Next thread: [MUG] Too many levels of recursion,
|
(NAG)
