function [nx,xmax,tmax,dt,alpha,mu,dx]=default

nx=600;
xmax=20*pi;
dx=xmax/nx;
tmax=20;
dt=0.05;
alpha=1;
mu=0.5;

