DMUG-Archiv 2010

Frühere   Chronologischer Index   Spätere
Vorherige   Thematischer Index   Nächste

Re: NDSolve für Dirichlet-Randbedingung

Hallo Frank

ich habe folgendes Problem mit Mathematica 7.01.0 nicht lösen können:

NDSolve[{0==D[x[r]x'[r],r], x[0]==100, x[1]==10}, x[r], {r,0,1}]

Diese Aufgabe ist äquivalent zu

In[68]:= Clear[y, x];
NDSolve[{y[x] y'[x] == -4950, y[0] == 100}, y[x], {x, 0, 1}]
Out[69]= {{y[x] -> InterpolatingFunction[{{0., 1.}}, <>][x]}}

falls man die Konstante -4950 kennt. Zu dem Zweck tippt man <snip>

Noch simpler ist

In[1]:= DSolve[{y[x] y'[x] == c, y[0] == p}, y[x], x]
Out[1]= {{y[x] -> -Sqrt[p^2 + 2 c x]}, {y[x] -> Sqrt[p^2 + 2 c x]}}

In[2]:= (* Suche das c so, dass y[1] == q wird *)
Solve[-Sqrt[p^2 + 2 c] == q, c]
Out[2]= {{c -> 1/2 (-p^2 + q^2)}}

In[3]:= Solve[Sqrt[p^2 + 2 c] == q, c]
Out[3]= {{c -> 1/2 (-p^2 + q^2)}}

++Gruss
Udo.



Verweise:
Frühere   Chronologischer Index   Spätere
Vorherige   Thematischer Index   Nächste

DMUG DMUG-Archiv, http://www.mathematica.ch/archiv.html