500                     ntmax (max no of time steps)
  6000.0                  tmax: max time for integration
     0.0                  xmin
     1.0                  xmax
    .10e-4                alpha: diffusion coefficient 
   0.5                    s = alpha*dt/dx/dx
    8                     nout: output every nout steps
   100.0                  f0: normalization of Temperature (t=f0*fnew)
     0.0                  finit: parameter for initial condition
     0.5                  fb0(1): boundary value at xmin (initially)
     0.5                  fb0(2): boundary value at xmax (initially)
     1.0                  fb1(1): boundary value at xmin 
     1.0                  fb1(2): boundary value at xmax 
   20                     maxex: no of terms in exact solution