We introduce a new high accuracy second-order finite-difference time-domain (FDTD) algorithm based on non-standard finite differences, and use it to investigate light propagation in subwavelength structures such as photonic crystals.