The partial differential Equation (REF ) is solved numerically to
obtain the time dependent behaviour of an accretion disc using the standard
finite differencing scheme with the courant condition imposed as described
in [1]}.
At an outer
boundary radius, the disc values are fixed to the steady state ones
while for the inner one we maintain the viscous stress free
condition.
