5        me parameter controlling choice of method (~> follow.table)
     2        ipr controls output, 1~>td solution
    11        jmax, no gridpoints along x
    10        maxex. no terms in exact sol
   500        nmax, max no of time steps
   0.01        alph thermal diffusivity alpha
   1.0        s=alpha*dt/dx**2
  12.0        tmax, max time
   4.5        tst initial time
   0.00       gam parameter gamma used in (7.26)
Various schemes implemented in program DIFIM:
  ME Description:
     1 fdm-2nd order ( 0 , 1  ,  0  )       bet = 0.5 + gam
     2 fem-2nd order (1/6, 2/3, 1/6 )	    bet = 0.5 + gam
     3 fdm-4th order ( 0 , 1  ,  0  )	    bet = 0.5 + gami - 1/(12s)
     4 fem-4th order (1/6, 2/3, 1/6 )	    bet = 0.5 + gaml + 1/(12s)
     5 composite     (1/12,5/6, 1/12)	    bet = 0.5 + gam
Other Variables:
  bet parameter beta used in (7.26) and table 7.4
  emx mass operator, M_x
  elx difference operator, L_xx
  ad, bd, cd ~> aj, bj and cj used in (7.31)