function [u]=inital(n,t)
%
% Function inital sets the initial condition vector
% of the linear PDE problem
%
% Problem parameters
  xl=0.0;
  xu=1.0;
%
% Initial condition
  for i=1:n
    x=xl+(i-1)/(n-1)*(xu-xl);
    u(i)=sin(pi*x);
  end

