function [nx,xmax,tmax,dt,dx,D]=default

nx=10;
xmax=10;
dx=xmax/nx;
tmax=20;
dt=0.02;
D=-.4;

