public void TestCase()
{
RK4Solver solver = new RK4Solver(ODEFunction, 0, Vector <double> .Build.DenseOfArray(new double[] { 0 }));
//solver.Init(0, Vector<double>.Build.DenseOfArray(new double[] { 0.0 }));
solver.Step(1);
Console.WriteLine(solver.State [0]);
solver.Step(1);
Console.WriteLine(solver.State [0]);
solver.Step(1);
Console.WriteLine(solver.State [0]);
}
Omega; // rotational velocity
public InflowModel()
{
solver = new RK4Solver(ODEFunction, 0, Vector <double> .Build.DenseOfArray(new double[] { 0 }));
mu = lambdabar = vbar = dCt = Ct0 = Omega = 0;
}