Merci Tom! Das hab ich gesucht !
Sorry das ich erst jetzt Antworte und bewerte.... ich mach das in meiner Freizeit und komme nicht immer direkt dazu zu antworten.
Ich habe jetzt nochmal mit den Rückgabewerten rumgespielt und habe jetzt mit einem Array gearbeitet. Der Grund hierfür ist, dass ich mit einem System aus ODE's arbeite....d.h.
die Funktionen sind von einander abhängig in dem Rechenvorgang.
das sieht in etwa so aus:
y1 = (a*y1)/y2
y2 = (b*y1)/y2
das Beispiel dazu siehst du hier:
{
public class Rechnen
{
public double[] Funktion(double a, double b)
{
double[] Calci;
Calci = new double[2];
Calci[0] = a + Calci[1];
Calci[1] = a + Calci[0];
return Calci;
}
}
public class Program
{
static void Main(string[] args)
{
Rechnen Test = new Rechnen();
double[] nummern;
nummern = Test.Funktion(3, 5);
Console.WriteLine("{0,3:F2}{1,10:F2}",nummern[0],nummern[1]);
Console.ReadKey();
}
}
}
Funktioniert wunderbar. Das muss ich jetzt in diesen Block einarbeiten:
{
// Definition des Runge Kutta Rechenverfahrens
class RungeKutta_Methode
{
public delegate double GetDerivative(double x, double y);
GetDerivative getDerivative;
double x0, y0, h;
double d0;
public RungeKutta_Methode(GetDerivative getDerivative, double x0,
double y0, double h, double d0)
{
this.getDerivative = getDerivative;
this.x0 = x0;
this.y0 = y0;
this.h = h;
this.d0 = d0;
}
public int Step(out double x1, out double y1, out double d1)
{
double k1 = d0;
double k2 = getDerivative(x0 + h / 2, y0 + h / 2 * k1);
double k3 = getDerivative(x0 + h / 2, y0 + h / 2 * k2);
double k4 = getDerivative(x0 + h, y0 + h * k3);
y1 = y0 + h / 6 * (k1 + 2 * k2 + 2 * k3 + k4);
x1 = x0 + h;
d1 = getDerivative(x1, y1);
x0 = x1;
y0 = y1;
d0 = d1;
return 0;
}
}
// Beschreibung des Rechenverfahrens und die Berechnungsvorschrift für den ersten Runge Kutta Schritt (Randbedindung)
class Funktion
{
public double[] funktion(double x, double y)
{
double[] fun;
fun = new double[1];
//double u = q / c;
fun[1] = y + Math.Exp(x);
return fun;
}
}
class Program
{
static void Main(string[] args)
{
// Programmtest
Funktion Test = new Funktion();
double x0 = 0;
double x_max = 0.2;
double y0 = 1;
double h = 0.05;
double[] i0;
i0 = Test.funktion(x0, y0);
RungeKutta_Methode integratorRK4D = new RungeKutta_Methode(Test.funktion, x0, y0, h, i0[1]);
Console.WriteLine("TEST Funktion");
Console.WriteLine("x y y'");
Console.WriteLine("-------------------");
Console.WriteLine("{0,3:F2}{1,10:F6}", x0, y0);
double x, y, i;
do
{
integratorRK4D.Step(out x, out y, out i);
Console.WriteLine("{0,3:F2}{1,10:F6}", x, y);
} while (x < x_max);
Console.WriteLine();
Console.ReadKey();
}
}
}
In der Klasse Funktion hab ich jetzt erstmal nur eine Gleichung [1] geschrieben.
Das Problem liegt bei
RungeKutta_Methode integratorRK4D = new RungeKutta_Methode(Test.funktion, x0, y0, h, i0[1]);
Hier sagt er mir, dass "Test.funktion" den falschen Rückgabetyp hat ?!
Ich stehe etwas auf dem Schlauch... :(
