An Algorithm for Qualitative Simulation of Gene Regulatory Networks with Steep Sigmoidal Response Functions