//*********************************************************************************

// The procedure that sets final values.  It fills the array u with the final values
// starting with u(xMin).  There are nu values in all, with x(k+1) = x(k) + dx.

int FinalValues( double* u, double xMin, double dx, int nu );

//*********************************************************************************

// The procedure that takes a time step.  uOld holds the values
// at the old time.  uNew gets the values at the new time.
// It uses nu values of uOld to compute nu-2 values of uNew, with
// one value of uNew missing from each end.

int TimeStep( double* uNew, double* uOld, double dx, double dt, double D, int nu );