Indefinite integrals involving the Jacobi Zeta and Heuman Lambda functions

Jacobian elliptic functions are used to obtain formulas for deriving indefinite integrals for the Jacobi Zeta function and Heuman's Lambda function. Only sample results are presented, mostly obtained from powers of the twelve Glaisher elliptic functions. However, this sample includes all such integrals in the literature, together with many new integrals. The method used is based on the differential equations obeyed by these functions when the independent variable is the argument u of elliptic function theory. The same method was used recently, in a companion paper, to derive similar integrals for the three canonical incomplete elliptic integrals.