Re: NDSolve für Dirichlet-Randbedingung

["Udo und Susanne Krause" <su.krause@XXXXXXX.ch>]

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.