How should I calculate the stationary value of ## S[y] ##?

  • #1
Math100
797
221
Homework Statement
Use the method of Lagrange multipliers to find the function ## y(x) ## that makes the functional ## S[y]=\int_{1}^{2}x^2y^2dx ## stationary subject to the two constraints ## \int_{1}^{2}xydx=1 ## and ## \int_{1}^{2}x^2ydx=2 ##. Calculate the stationary value of ## S[y] ##.
Relevant Equations
None.
Consider the functional ## S[y]=\int_{1}^{2}x^2y^2dx ## stationary subject to the two constraints ## \int_{1}^{2}xydx=1 ## and ## \int_{1}^{2}x^2ydx=2 ##.
Then the auxiliary functional is ## \overline{S}[y]=\int_{1}^{2}(x^2y^2+\lambda_{1}xy+\lambda_{2}x^2y)dx, y(1)=y(2)=0 ## where ## \lambda_{1} ## and ## \lambda_{2} ## are the Lagrange multipliers.
By definition, the Euler-Lagrange equation is ## \frac{d}{dx}(\frac{\partial F}{\partial y'})-\frac{\partial F}{\partial y}=0, y(a)=A, y(b)=B ## for the functional ## S[y]=\int_{a}^{b}F(x, y, y')dx, y(a)=A, y(b)=B ##.
Let ## F(x, y, y')=x^2y^2+\lambda_{1}xy+\lambda_{2}x^2y ##.
This gives ## \frac{\partial F}{\partial y}=2x^2y+\lambda_{1}x+\lambda_{2}x^2 ##.
Thus, the Euler-Lagrange equation is ## 2x^2y+\lambda_{1}x+\lambda_{2}x^2=0\implies 2x^2y=-\lambda_{1}x-\lambda_{2}x^2\implies y=-\frac{\lambda_{1}}{2x}-\frac{\lambda_{2}}{2} ##.
The first constraint ## \int_{1}^{2}xydx=1 ## gives ## 1=\int_{1}^{2}x[-\frac{\lambda_{1}}{2x}-\frac{\lambda_{2}}{2}]dx\implies 1=-\frac{\lambda_{1}}{2}-\frac{3\lambda_{2}}{4} ## and the second constraint ## \int_{1}^{2}x^2ydx=2 ## gives ## 2=\int_{1}^{2}x^2[-\frac{\lambda_{1}}{2x}-\frac{\lambda_{2}}{2}]dx\implies 2=-\frac{3\lambda_{1}}{4}-\frac{7\lambda_{2}}{6} ##, so ## \lambda_{1}=16 ## and ## \lambda_{2}=-12 ##.
Hence ## y(x)=-\frac{\lambda_{1}}{2x}-\frac{\lambda_{2}}{2}\implies y(x)=-\frac{8}{x}+6 ##.
Observe that ## S[y]=\int_{1}^{2}x^2y^2dx\implies S[y]=\int_{1}^{2}x^2(-\frac{8}{x}+6)^2dx\implies S[y]=\int_{1}^{2}x^2(\frac{64}{x^2}-\frac{96}{x}+36)dx\implies S[y]=4 ##.
Therefore, the function ## y(x) ## that makes the functional ## S[y]=\int_{1}^{2}x^2y^2dx ## stationary subject to the two constraints ## \int_{1}^{2}xydx=1 ## and ## \int_{1}^{2}x^2ydx=2 ## is ## y(x)=-\frac{8}{x}+6 ## and the stationary value of ## S[y] ## is ## 4 ##.

I just want to know if the work and answer shown above is correct or not. Please check/verify/confirm to see if this is correct or not.
 
Last edited:
Physics news on Phys.org
  • #2
This is the correct method; I hve not checked your arithmetic.

Note that [itex]y'[/itex] does not appear in [itex]F[/itex], so you do not end up with an ODE and you do not have boundary conditions.
 
  • Like
Likes Math100
  • #3
pasmith said:
This is the correct method; I hve not checked your arithmetic.

Note that [itex]y'[/itex] does not appear in [itex]F[/itex], so you do not end up with an ODE and you do not have boundary conditions.
Thank you so much for confirming/verifying!
 

FAQ: How should I calculate the stationary value of ## S[y] ##?

What is the functional ## S[y] ##?

The functional ## S[y] ## is an expression that maps a function ## y(x) ## to a real number, typically represented as an integral involving ## y(x) ## and its derivatives. It is often written in the form ## S[y] = \int_{a}^{b} L(x, y, y') \, dx ##, where ## L ## is a given function called the Lagrangian.

What does it mean for a functional to have a stationary value?

A functional has a stationary value if its first variation vanishes for a given function ## y(x) ##. This means that small perturbations or variations in ## y(x) ## do not change the value of the functional to the first order. In mathematical terms, the condition for a stationary value is given by ## \delta S[y] = 0 ##.

What is the Euler-Lagrange equation?

The Euler-Lagrange equation is a differential equation derived from the condition that the functional ## S[y] ## has a stationary value. It is given by ## \frac{\partial L}{\partial y} - \frac{d}{dx} \left( \frac{\partial L}{\partial y'} \right) = 0 ##. Solving this equation provides the function ## y(x) ## that makes the functional stationary.

How do I apply boundary conditions when calculating the stationary value?

When calculating the stationary value of a functional, boundary conditions must be specified to ensure a well-defined problem. These conditions are typically given as ## y(a) = y_a ## and ## y(b) = y_b ##. They are used in the integration by parts process to eliminate boundary terms, ensuring that the variations in ## y(x) ## at the endpoints do not affect the stationary condition.

Can you provide an example of calculating the stationary value of a functional?

Sure! Consider the functional ## S[y] = \int_{0}^{1} \left( (y')^2 + y^2 \right) \, dx ##. The Lagrangian is ## L = (y')^2 + y^2 ##. The Euler-Lagrange equation is ## \frac{\partial L}{\partial y} - \frac{d}{dx} \left( \frac{\partial L}{\partial y'} \right) = 0 ##, which simplifies to ## y'' - y = 0 ##. Solving this differential equation with boundary conditions ## y(0) = 0 ## and ## y(1) = 0 ##, we find that the solution is ## y(x) = 0 ##. Thus, the stationary value is achieved when ## y(x) ##

Back
Top