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.