summaryrefslogtreecommitdiffstats
diff options
context:
space:
mode:
-rw-r--r--tests/EquationSystem.cpp10
1 files changed, 7 insertions, 3 deletions
diff --git a/tests/EquationSystem.cpp b/tests/EquationSystem.cpp
index 0a929aa..6919506 100644
--- a/tests/EquationSystem.cpp
+++ b/tests/EquationSystem.cpp
@@ -25,7 +25,11 @@
typedef vZ::EquationSystem<2> Y;
-// y'' = y (y == C*exp(t))
+// y'' = y (y == C*exp(t) + D*exp(-t))
+//
+// Split as:
+// y' = v
+// v' = y
Y
f(double t, Y y)
{
@@ -40,7 +44,7 @@ main()
{
Y y;
y[0] = 1.0;
- y[1] = 1.0;
+ y[1] = 0.0;
vZ::GenericDP45Integrator<Y> integrator(f);
integrator.tol(1e-6)
.y(y)
@@ -50,7 +54,7 @@ main()
integrator.integrate(2.0);
double actual = integrator.y()[0];
- double expected = std::exp(2.0);
+ double expected = std::cosh(2.0);
std::cout << std::setprecision(10)
<< "Numerical: " << actual << std::endl