public void Param_c_is_zero_Test()
{
// r = 2, a_0 = a_1 = 1, y(t) = e^{0*t} = 1
var solver = new OdeSolverExp(0, 1, 1);
var pCoeffs = Range(0, 10).Select(solver.PCoeffs).ToArray();
Assert.That(pCoeffs, Is.EqualTo(Repeat(0d, 10).ToArray()).Within(Tolerance));
var yCoeffs = solver.CalcCoeffs(2);
Assert.That(yCoeffs, Is.EqualTo(new[] { 1d, 0 }).Within(Tolerance));
yCoeffs = solver.CalcCoeffs(6);
Assert.That(yCoeffs, Is.EqualTo(new[] { 1d, 0, 0, 0, 0, 0 }).Within(Tolerance));
}
public void PCoeffTwoTest()
{
var solver = new OdeSolverExp(0, 1, 1);
Assert.That(solver.PCoeffs(2), Is.EqualTo(0).Within(Tolerance));
}
