inital:=proc(neqn,t,u0)
#
# Function inital sets the initial condition vector
# of the nonlinear PDE
#
# Type variables
  local ui, i:
#
# Problem parameters
  ui:=298.0:
#
# Initial condition
  for i from 1 to neqn do
    u0[i]:=ui:
  end do:
#
# End of inital
  end: