bvp4c mit Vektor in Funktion < Matlab < Mathe-Software < Mathe < Vorhilfe
|
Hallo,
bvp4c ist ja ein toller solver für Randwertprobleme und kann sogar unbekannte Parameter finden.
Mein Problem ist allerdings viel simpler:
Ich habe die DGL:
[mm] U''(r)=-4\pin(r), [/mm] wobei n(r) ein mir bekannter 1D-Vektor ist
Die Lösung hat die Form [mm] U(r)=\alpha\cdot [/mm] r
Außerdem kenne ich die Randbedingungen U(0)=0 und U(r_max)=c,
wobei c eine bekannte Konstante ist.
Ich dachte es müsste wie folgt gehen:
f = @(r,U) [ u(2); -4*pi.*r.*n];
g = @(ya,yb) [ya(1); yb(1)-c];
solinit=bvpinit(linspace(0,r_max,10),[0 0]);
sol = bvp4c(f,g,solinit);
Allerdings bekomme ich immer den Fehler:
Error using vertcat
Dimensions of matrices being concatenated are not consistent.
Ideen?
|
|
| |
|
Hallo.
Mit welcher Version von Matlab arbeitest Du?
Wie groß ist r bez. rmax? Oder habe ich das überlesen?
Beste Grüße.
|
|
|
| |
|
| Status: |
(Mitteilung) Reaktion unnötig  | | Datum: | 14:20 Sa 27.08.2016 | | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|
