fix

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);
 				}
 			}