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

