From: W. Trevor King Date: Wed, 13 Feb 2013 16:30:42 +0000 (-0500) Subject: K should be 10, not -5 X-Git-Url: http://git.tremily.us/?a=commitdiff_plain;h=ba14c6205256e048f54c20619f3ad9fca1b8b21c;p=assignment-template.git 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 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