
In this work we present a theoretical model for differentiable programming.
We construct an algebraic language that encapsulates formal semantics of
differentiable programs by way of Operational Calculus. The algebraic nature of
Operational Calculus can alter the properties of the programs that are
expressed within the language and transform them into their solutions.
In our model programs are elements of programming spaces and viewed as maps
from the virtual memory space to itself. Virtual memory space is an algebra of
programs, an algebraic data structure one can calculate with. We define the
operator of differentiation ($\partial$) on programming spaces and, using its
powers, implement the general shift operator and the operator of program
composition. We provide the formula for the expansion of a differentiable
program into an infinite tensor series in terms of the powers of $\partial$. We
express the operator of program composition in terms of the generalized shift
operator and $\partial$, which implements a differentiable composition in the
language. Such operators serve as abstractions over the tensor series algebra,
as main actors in our language.
We demonstrate our models usefulness in differentiable programming by using
it to analyse iterators, deriving fractional iterations and their iterating
velocities, and explicitly solve the special case of ReduceSum.

In this article we use algebrogeometric tools to describe the structure of a
superintegrable system. We study degenerate Neumann system with potential
matrix that has some eigenvalues of multiplicity greater than one. We show that
the degenerate Neumann system is superintegrable if and only if its spectral
curve is reducible and that its flow can be linearized on the generalized
Jacobians of the spectral curve. We also show that the generalized Jacobians of
the hypereliptic component of the spectral curves are models for the minimal
invariant tori of the flow. Moreover the spectral invariants generate local
actions that span the invariant tori, while the moment maps for the rotational
symmetries provide additional first integrals of the system. Using our results
we reproduce already known facts that the degenerate Neumann system is
superintegrable, if its potencial matrix has eigenvalues of multiplicity
greater or equal than $3$.

In this paper we study the Neumann system, which describes the harmonic
oscillator (of arbitrary dimension) constrained to the sphere. In particular we
will consider the confluent case where two eigenvalues of the potential
coincide, which implies that the system has S^{1} symmetry. We will prove
complete algebraic integrability of confluent Neumann system and show that its
flow can be linearized on the generalized Jacobian torus of some singular
algebraic curve. The symplectic reduction of S^{1} action will be described and
we will show that the general Rosochatius system is a symplectic quotient of
the confluent Neumann system, where all the eigenvalues of the potential are
double. This will give a new mechanical interpretation of the Rosochatius
system.