Given the initial shape of the longitudinal axis of a curved bar is given by
y₀(x) = a*sin(Pi*x/L)

where a is a constant and L is the length of the bar

and proceeding to calculate the deflections due to deformation from the diff. eq.

EI (dy₁/dx) = -P(y₀+y₁)

the ode for deflections is obtained

dy₁/dx + k*y₁ = -k*a* sin(Pi*x/L)

where k is a constant = sqrt (P/EI)

I am trying to develop a solution to a problem which consists of taking a given vector field (force) and converting that field from one domain to another while maintaining an equivalent system of forces. The real world problem occurs in posttensioning design where steel tendons are place in concrete members and then a tensile force is applied to the tendon, which is then restrained at the end points, creating a compressive vector field in the concrete.

I am trying to solve a one-dimensional heat equation problem with a piecewise continuous input of temperatures.  Dr.

I would like to use Maple to provide an independent check on a finite element model solution (ANSYS) of a thermal transient analysis of a reinforced concrete shell.  The following information is known about the shell:

Recently Dr. Israel responded to my request for help in extending the EllipticF function past the limit of Pi/2 for the amplitude (see topic titled Elliptic Integrals). After reviewing A&S Chapter 17, I have tried to duplicate the results using the JacobiAM function in Maple. The help page for this function indicates that there is no limit on the amplitude. The attached worksheet evaluates the form suggested by A&S, Eq. 17.4.3 and the JacobiAM function. It is interesting to note the only when the argument given to the EllipticF function is equal to the remainder of Pi/2 - beta that the two expressions are equal. I would think that the JacobiAM form is a more compact representation of EllipticF for amplitudes greater than Pi/2. Are the two functions equivalent as used in the worksheet?