K should be 10, not -5

The homework calls for an energy barrier of five units.

diff --git a/one_gaussian_bump.c b/one_gaussian_bump.c
index 1124022144d069ec733bcdb6a8286e7c97b9ccdf..081ab204581053531660116c4e6efa268489d222 100644 (file)
--- a/one_gaussian_bump.c
+++ b/one_gaussian_bump.c
@@ -323,7 +323,7 @@ int main(int argc, char **argv) {
        system.dx_dt_fn = &single_gaussian_dx_dt;
        system.energy_fn = &single_gaussian_energy;
        system.dx_dt_state = (void *)&gaussian;
-       gaussian.k = -5;
+       gaussian.k = 10;
        gaussian.b = 1;
        state.x[2] = -sqrt(2*0.5/system.m);  /* v = -sqrt(2E/m),  with E = 0.5 */
 

diff --git a/phase.gp b/phase.gp
index 7cc02aa459822a572e275cf3324a762393104946..c6861bc12fccad90c7c2c5b330481244c956daf3 100644 (file)
--- a/phase.gp
+++ b/phase.gp
@@ -4,8 +4,7 @@ set title 'Phase space'
 set xlabel 'x'
 set ylabel 'v_x'
 unset key
-set yrange [-4.5:0.5]
-k = -5
+k = 10
 b = 1
 V(r) = 0.5 * k * b**2 * exp(-r**2/b**2)
 set style data lines