TITLE 'Mg-diffusion'     { the problem identification }
COORDINATES cartesian1  { coordinate system, 1D,2D,3D, etc }
VARIABLES        { system variables }
  C              { Mg concentration }
! SELECT         { method controls }
 DEFINITIONS    { parameter definitions }
 D=4.60
 C0=1E19
!INITIAL VALUES
 
EQUATIONS        { PDE's, one for each variable }
  dxx(D*C)=dt(C)
! CONSTRAINTS    { Integral constraints }
BOUNDARIES       { The domain definition }
  REGION 1       { For each material region }
    START(0)  POINT VALUE(C)=C0  { Walk the domain boundary }
    LINE TO (80)
 TIME 0 TO 7.5  by 0.01 { if time dependent }
MONITORS         { show progress }
PLOTS    { save result displays }
  FOR t=7.5
  ELEVATION(C) FROM(0) TO (80)  LINLOG  AS  'CONCENTRATION DISTRIBUTION'
  TABLE(C)
END