Differentialgleichungssysteme < Sonstiges < Hochschule < Informatik < Vorhilfe
|
Status: |
(Frage) überfällig | Datum: | 17:16 Mi 05.03.2008 | Autor: | Lufti |
Aufgabe | Schreiben Sie ein Fortran-Programm, das ein nichtlineares Differentialgleichungssystem (Anfangswertproblem) loesen kann.
Folgende Methoden sollen verwendet werden:
forward Euler
backward Euler
Heun
Runge-Kutta 4. Ordnung |
Hallo,
ich habe damit angefangen, ein Programm zu schreiben, das eine einzelne Differentialgleichung loesen kann. So sieht z.B der Teil fuer forward Euler aus:
!EULER_FORWARD
x=x0
y=y0
i=0
WRITE (*,*)'x=',x,'y=',y
DO WHILE (i<n)
y=euler_forward(h,x,y)
x=x+h
WRITE (*,*)'x=',x,'y=',y
i=i+1
END DO
FUNCTION euler_forward (h,x,y)
IMPLICIT NONE
REAL, INTENT(IN)::h,x,y
REAL, external::f
REAL:: euler_forward
euler_forward=y+h*f(x,y)
RETURN
END FUNCTION euler_forward
FUNCTION f (x,y)
IMPLICIT NONE
REAL, INTENT(IN):: x,y
REAL::f
f=(y**(-2))*(cos(x)+1)
RETURN
END FUNCTION f
Ich muesste das doch relativ einfach umschreiben koennen, oder nicht? Ich wollte y als Feld deklarieren und dann eine zweite Schleife laufen lassen, die mir dann an jeder x-Stelle die Werte fuer y1, y2 usw ausrechnet. Da bin ich dann aber nicht weitergekommen, weil ich ja irgendwie mehrere Funktionen brauche und ich weiss nicht, wie meine Verfahren diese aufrufen koennen.
Kann mir jemand helfen?
Ich habe diese Frage in keinem Forum auf anderen Internetseiten gestellt.
|
|
|
|
Status: |
(Mitteilung) Reaktion unnötig | Datum: | 17:20 Sa 08.03.2008 | Autor: | matux |
$MATUXTEXT(ueberfaellige_frage)
|
|
|
|