From ba14c6205256e048f54c20619f3ad9fca1b8b21c Mon Sep 17 00:00:00 2001 From: "W. Trevor King" Date: Wed, 13 Feb 2013 11:30:42 -0500 Subject: [PATCH] K should be 10, not -5 The homework calls for an energy barrier of five units. --- one_gaussian_bump.c | 2 +- phase.gp | 3 +-- 2 files changed, 2 insertions(+), 3 deletions(-) diff --git a/one_gaussian_bump.c b/one_gaussian_bump.c index 1124022..081ab20 100644 --- 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 7cc02aa..c6861bc 100644 --- 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 -- 2.26.2