An Algorithm for Symbolic Computing of Singular Limits of Dynamical Systems