function [xl,xu]=par
%
% Function par sets the parameters for the parabolic 
% PDE problem
%
  xl=0.0;
  xu=1.0;


