A physics-informed neural network puts the differential equation INTO the loss. Write the equation as a residual, for u'' + ω²u = 0 the residual is f(t) = u''(t) + ω²u(t), and drive f to zero at randomly sampled 'collocation' points, adding terms for the boundary or initial conditions. Because there is no data term inside the domain, the network is supervised entirely by physics; it interpolates the true solution between the conditions. The technical key, and the reason this was new: the loss needs the network's OWN derivatives with respect to its input, u'(t) and u''(t), used as outputs. ferromotion computes them exactly by propagating a second-order jet (u, u', u'') through the network on the same autodiff tape that computes the parameter gradient, so the whole thing trains with exact gradients, no finite differences. Two payoffs beyond elegance: a trained PINN evaluates a PDE solution in microseconds where a classical solver takes minutes, and, the version that matters for robotics, you can mix a physics-residual loss where you have no data with a data loss where you do, so the model extrapolates safely into the unmeasured regions.