In the paper, a novel algorithm employing pseudo-spectral approach is developed for the PKN model of hydrofracturing. The respective solvers compute both the solution and its temporal derivative. In comparison with conventional solvers, they demonstrate significant cost effectiveness in terms of balance between the accuracy of computations and densities of the temporal and spatial meshes. Various fluid flow regimes are considered.