# You may run this example by doing: # # ode < soliton.ode | graph -T X -C # # or alternatively, to create a real-time plot, # # ode < soliton.ode | graph -T X -C -x 5 15 5 -y -3 3 # This example numerically solves an initial value problem arising in # John Selker's soliton thesis: # # x'' = (2 / t) * x' + B * x # B = -1 # x(5) = 0 # x'(5) = 1 x' = vx vx' = 2 / t * vx + B * x B = -1 x = 0 vx = 1 print t, x step 5, 15