bux fix for single segment pendula

main
Benjamin Kraft 1 year ago
parent c9a281fbb3
commit 719c710618
  1. 9
      src/Pendulum.cpp

@ -71,11 +71,10 @@ void Pendulum::update(double h, double g) {
if (i > 0){ if (i > 0){
V[i - 1] = (p1 - X[i - 1]) / h; V[i - 1] = (p1 - X[i - 1]) / h;
X[i - 1] = p1; X[i - 1] = p1;
}
if (i == X.size() - 1){ if (i == X.size() - 1){
V[i] = (p2 - X[i]) / h; V[i] = (p2 - X[i]) / h;
X[i] = p2; X[i] = p2;
}
} }
p1 = p2; p1 = p2;

Loading…
Cancel
Save