diff --git a/src/Solver.cpp b/src/Solver.cpp index fbbbec4..9f7beb5 100644 --- a/src/Solver.cpp +++ b/src/Solver.cpp @@ -66,7 +66,7 @@ void Solver::SolveExplicit(System& sys, double tstop) const { double tx = (node->r()->T() - 2 * node->T() + node->l()->T()) / pow(sys.step(), 2); double ty = (node->u()->T() - 2 * node->T() + node->d()->T()) / pow(sys.step(), 2); - double t = delta * (tx + ty) * sys.a()/pow(sys.step(), 2) + node->T(); + double t = delta * (tx + ty) * sys.a() + node->T(); node->SetT(t); } }