function [u]=inital(n,t)
%
% Function inital sets the initial condition vector
% for the nonlinear PDE
%
% Problem parameters
  [alpha,k,sigma,L,ui,ua]=par1;
%
% Inital condition
  for i=1:n
    u(i)=ui;
  end
%
% End of inital
