Hallo Frank,
Meine Frage ist also, ob man mit NDSolve eine DGL mit
Dirichlet-Randbedingung lösen kann. Scheinbar kann das NDSolve alleine
nicht.
Dieser Schein trügt schlussendlich,
In[1]:= sol = NDSolve[{D[y[x] D[y[x], x], x] == 0, y[0] == 100, y[1] ==
10},
y, {x, 0, 1},
Method -> {"Shooting",
"StartingInitialConditions" -> {y[0] == 1/1000, y'[0] == -1/100}}]
In[3]:= Plot[sol[[1, 1, 2]][x], {x, 0, 1}]
funktioniert mit Bemerkungen. Das hatte ich schon vorige Woche probiert,
dann dem Wolfram Support einen bug report gesendet, der dort nicht
reproduzierbar war, jetzt geht's hier in einem frischen Notebook auch,
hallelujah.
Gruss
Udo.
|