public void Constructor_SetsDim()
{
//Create an ODE initial value problem
ODEInitialValueProblem ivp = new ODEInitialValueProblem(3, 0.2, 0);
Assert.AreEqual(3, ivp.GetDim());
}
public void Constructor_WillNotSetDimToZero()
{
//Create an ODE initial value problem
ODEInitialValueProblem ivp = new ODEInitialValueProblem(0, 0.1, 0);
Assert.AreNotEqual(0, ivp.GetDim());
}
//Returns the current state vector
public VectorND GetCurrentState()
{
if (ivp.GetDim() == 3)
{
return(new VectorND(GetCurrentState("x"),
GetCurrentState("y"),
GetCurrentState("z")));
}
else
{
return(new VectorND(GetCurrentState("x"),
GetCurrentState("x'"),
GetCurrentState("y"),
GetCurrentState("y'"),
GetCurrentState("z"),
GetCurrentState("z'")));
}
}