In this article, we solve in closed form a system of nonlinear differential equations modelling the elastica in space of a thin, flexible, straight rod, loaded by a constant thrust at its free end. Common linearizations of strength of materials are of course not applicable any way, because we analyze great deformations, even if not so large to go off the linear elasticity range. By passing to cylindrical coordinates rho, theta, z, we earn a more tractable differential system evaluating rho as elliptic function of polar anomaly theta and also providing z through elliptic integrals of I and III kind. Deformed rod's centerline is then completely described under both tensile or compressive load. Finally, the planar case comes out as a degeneracy, where the Bernoulli lemniscatic integral appears.