We have developed an finite-difference time-domain program that can analyze photonic devices with gain and/or dispersion. As an example, a two-dimensional photonic-crystal laser is simulated. The simulation can show the relaxation oscillation behavior at extremely high current injection.