NeuralPDR: neural differential equations as surrogate models for photodissociation regions