inital:=proc(neqn,t,u0)
#
# Function inital sets the initial condition vector
# of the linear PDE
#
# Type variables
  local xl, xu, pi, x, i:
#
# Problem parameters
  xl:=0.0:
  xu:=1.0:
#
# Initial condition
  pi:=evalf(Pi):
  for i from 1 to neqn do
    x:=xl+(i-1)/(neqn-1)*(xu-xl):
    u0[i]:=sin(pi*x):
  end do:
#
# End of inital
  end: