// laxwen1.hpp		1-D Lax Wendroff (single step) method//		u_t + f_x = rhs//				where rhs is optional, assumed zero otherwise#ifndef _LAX_WENDROFF1_HPP_#define _LAX_WENDROFF1_HPP_#include <dynamics.hpp>#ifndef SPACE_DERIV#define SPACE_DERIVtypedef float (*SpaceDeriv)(Vector&,int);#endifclass LaxWendroff1 : public Dynamics{	private:          float r;	  Vector u_new;	  SpaceDeriv f;	  SpaceDeriv rhs;	public:	  LaxWendroff1(SpaceDeriv f, float dx_in, float dt_in, int n_in,		      SpaceDeriv user_rhs = NULL);	 ~LaxWendroff1() { }	  void step(const int steps = 1);};#endif