2                   method: 1-upw, 2-leapfrog, 3-l-w, 4-c-n
 600                   ntmax (max no of time steps)
 6.0                   tmax: max time for integration
 -2.0                  xmin
  4.0                  xmax
   .02                   alpha: diffusion coefficient 
  1.0                  u (velocity)
  0.1                  c = u*dt/dx
  0.02                 dt
  0                    cordt (dt from: 0-c or 1-dt directly
  0.0                  delta (param for mass operator)
  0.0                  que (param for upwind operator)
 10                    nout: output every nout steps
  1.0                  f0: normalization of function (t=f0*fnew)
   3                   init: 1-trunc sine, 2-sine, 3-step function
  1.0                  lambda: half wavelength for initial condition
   2                   parbound: 1-periodic, 2-dirichlet, 3-neumann
  1.0                  fb(1): boundary value or derivative at xmin 
  0.0                  fb(2): boundary value or derivative at xmax 
 101                   maxex: no of terms in exact solution<=ntex