function out = w10oneway_ode(t,u) % % w10oneway_ode.m -- djm -- 06 mar 06 % - one-way wave, method of lines % % - - - - - - - - - - - - - - - - - - - - - - - - - - - - - - global c dx out = -(c/2/dx)*([u(2:end) ; u(1)] - [u(end) ; u(1:end-1)]);