Direct and inverse interpolation for Jacobian elliptic functions, zeta function of Jacobi and complete elliptic integrals of the second kind (Q1352792)

The problems mentioned in this paper for the ellipsoidal shape of the internal boundary of the cavity or inclusion suggests the use of the ellipsoidal coordinates. The corresponding coordinate transformation, which relates the Cartesian coordinates \((x,y,z)\) to the ellipsoidal coordinates \(a_i\), \(i= 1,2,3\) is most conveniently introduced by using Jacobian elliptic functions. Computing the value of the Jacobian elliptic functions, given the argument \(u\) and the parameter \(m\), is a problem, whose solution can be found either tabulated in tables of elliptic functions or by use of existing software, such as Mathematica, etc. The inverse problem, finding the argument, given the Jacobian elliptic function and the parameter \(m\), is a problem whose solution is found only in tables of elliptic functions. Standard polynomial inverse procedures fail, due to ill conditioning of the system of linear equations of the unknown. This paper contains a numerical procedure for inverse interpolation which gives good results in the computation of the argument of the Jacobian elliptic function given the Jacobian elliptic function and the parameter. Also, a direct interpolation is described which gives the zeta function of Jacobi and the complete elliptic integral of the second kind given the argument and the parameter. These new interpolation procedures are important in problems involving cavities or inclusions of ellipsoidal shape encountered in the mechanical design of bearing, filters and composite materials. They are also important in the modeling of porosity of bones. This porosity may lead to osteoporosis, a disease which affects bone material density in human with bad consequences.