functionalDerivative

Find the functional derivative of the functional $$S[y]={\int }_b^{a}y(x)\sin (y(x))dx$$ with respect to the function $$y$$, where the integrand is $$f[y(x)]=y(x)\sin (y(x))$$.

Declare y(x) as a symbolic function and define f as the integrand of $$S$$. Use f and y as the parameters of functionalDerivative.

syms y(x)
f = y*sin(y);
G = functionalDerivative(f,y)

G(x) = $$\sin \left(y\left(x\right)\right)+\cos \left(y\left(x\right)\right)\hspace{0.17em}y\left(x\right)$$

Find the functional derivative of the functional $$S[u,v]={\displaystyle {\int }_b^a(u^2(x)\frac{dv(x)}{dx}+v(x)\frac{d^{2}u(x)}{d{x}^2})dx}$$ with respect to the functions $$u$$ and $$v$$, where the integrand is $$f[u(x),v(x),u^{\prime \prime }(x),v^{\prime }(x)]=u^2\frac{dv}{dx}+v\frac{d^{2}u}{d{x}^2}$$.

Declare u(x) and v(x) as symbolic functions, and define f as the integrand of $$S$$.

syms u(x) v(x)
f = u^2*diff(v,x) + v*diff(u,x,x);

Specify a vector of symbolic functions [u v] as the second input argument in functionalDerivative.

G = functionalDerivative(f,[u v])

G(x) = 
$$\left(\begin{array}{c}\frac{{\partial }^2}{\partial x^2}\mathrm{ }v\left(x\right)+2\hspace{0.17em}u\left(x\right)\hspace{0.17em}\frac{\partial }{\partial x}\mathrm{ }v\left(x\right)\\ \frac{{\partial }^2}{\partial x^2}\mathrm{ }u\left(x\right)-2\hspace{0.17em}u\left(x\right)\hspace{0.17em}\frac{\partial }{\partial x}\mathrm{ }u\left(x\right)\end{array}\right)$$

functionalDerivative returns a vector of symbolic functions containing the functional derivatives of the integrand f with respect to u and v, respectively.

Find the Euler–Lagrange equation of a mass m that is connected to a spring with spring constant k.

Define the kinetic energy T, potential energy V, and Lagrangian L of the system. The Lagrangian is the difference between the kinetic and potential energy.

syms m k x(t)
T = 1/2*m*diff(x,t)^2;
V = 1/2*k*x^2;
L = T - V

L(t) = 
$$\frac{m\hspace{0.17em}{\left(\frac{\partial }{\partial t}\mathrm{ }x\left(t\right)\right)}^2}{2}-\frac{k\hspace{0.17em}{x\left(t\right)}^2}{2}$$

In Lagrangian mechanics, the action functional of the system is equal to the integral of the Lagrangian over time, or $$S[x]={\int }_{t_1}^{t_2}L[t,x(t),\underset{}{\overset{\dot{}}{x}}(t)]dt$$. The Euler–Lagrange equation describes the motion of the system for which $$S[x(t)]$$ is stationary.

Find the Euler–Lagrange equation by taking the functional derivative of the integrand L and setting it equal to 0.

eqn = functionalDerivative(L,x) == 0

eqn(t) = 
$$-m\hspace{0.17em}\frac{{\partial }^2}{\partial t^2}\mathrm{ }x\left(t\right)-k\hspace{0.17em}x\left(t\right)=0$$

eqn is the differential equation that describes mass-spring oscillation.

Solve eqn using dsolve. Assume the mass m and spring constant k are positive. Set the initial conditions for the oscillation amplitude as $$x(0)=10$$ and the initial velocity of the mass as $$\underset{}{\overset{\dot{}}{x}}(0)=0$$.

assume(m,'positive')
assume(k,'positive')
Dx(t) = diff(x(t),t);
xSol = dsolve(eqn,[x(0) == 10, Dx(0) == 0])

xSol = 
$$10\hspace{0.17em}\cos \left(\frac{\sqrt{k}\hspace{0.17em}t}{\sqrt{m}}\right)$$

Clear assumptions for further calculations.

assume([k m],'clear')

The Brachistochrone problem is to find the quickest path of descent of a particle under gravity without friction. The motion is confined to a vertical plane. The time for a body to move along a curve $$y(x)$$ from point $$a$$ to $$b$$ under gravity $$g$$ is given by

Find the quickest path by minimizing the change in $$t$$ with respect to small variations in the path $$y$$. The condition for a minimum is $$\frac{\delta t}{\delta y}(x)=0$$.

Compute the functional derivative to obtain the differential equation that describes the Brachistochrone problem. Use simplify to simplify the equation to its expected form.

syms g y(x)
assume(g,'positive')
f = sqrt((1 + diff(y)^2)/(2*g*y));
eqn = functionalDerivative(f,y) == 0;
eqn = simplify(eqn)

eqn(x) = 
$$2\hspace{0.17em}y\left(x\right)\hspace{0.17em}\frac{{\partial }^2}{\partial x^2}\mathrm{ }y\left(x\right)+{\left(\frac{\partial }{\partial x}\mathrm{ }y\left(x\right)\right)}^2=-1$$

This equation is the standard differential equation of the Brachistochrone problem. To find the solutions of the differential equation, use dsolve. Specify the 'Implicit' option to true to return implicit solutions, which have the form $$F(y(x))=g(x)$$.

sols = dsolve(eqn,'Implicit',true)

sols = 
$$\begin{array}{l}\left(\begin{array}{c}y\left(x\right)=C_5+x\hspace{0.17em}\mathrm{i}\\ y\left(x\right)=C_4-x\hspace{0.17em}\mathrm{i}\\ {\sigma }_1=C_3-x\\ \frac{C_1+y\left(x\right)}{y\left(x\right)}=0\\ {\sigma }_1=C_2+x\end{array}\right)\\ \\ \mathrm{where}\\ \\ \mathrm{ }{\sigma }_1=C_1\hspace{0.17em}\mathrm{atan}\left(\sqrt{-\frac{C_1}{y\left(x\right)}-1}\right)-y\left(x\right)\hspace{0.17em}\sqrt{-\frac{C_1}{y\left(x\right)}-1}\end{array}$$

The symbolic solver dsolve returns general solutions in the complex space. Symbolic Math Toolbox™ does not accept the assumption that the symbolic function $$y(x)$$ is real.

Depending on the boundary conditions, there are two real-space solutions to the Brachistochrone problem. One of the two solutions below describes a cycloid curve in real space.

solCycloid1 = sols(3)

solCycloid1 = 
$$C_1\hspace{0.17em}\mathrm{atan}\left(\sqrt{-\frac{C_1}{y\left(x\right)}-1}\right)-y\left(x\right)\hspace{0.17em}\sqrt{-\frac{C_1}{y\left(x\right)}-1}=C_3-x$$

solCycloid2 = sols(5)

solCycloid2 = 
$$C_1\hspace{0.17em}\mathrm{atan}\left(\sqrt{-\frac{C_1}{y\left(x\right)}-1}\right)-y\left(x\right)\hspace{0.17em}\sqrt{-\frac{C_1}{y\left(x\right)}-1}=C_2+x$$

Another solution in real space is a horizontal straight line, where $$y$$ is a constant.

solStraight = simplify(sols(4))

solStraight = $$C_1+y\left(x\right)=0$$

To illustrate the cycloid solution, consider an example with boundary conditions $$y(0)=5$$ and $$y(4)=1$$. In this case, the equation that can satisfy the given boundary conditions is solCycloid1. Substitute the two boundary conditions into solCycloid1.

eq1 = subs(solCycloid1,[x y(x)],[0 5]);
eq2 = subs(solCycloid1,[x y(x)],[4 1]);

The two equations, eq1 and eq2, have two unknown coefficients, $$C_1$$ and $$C_3$$. Use vpasolve to find numerical solutions for the coefficients. Substitute these solutions into solCycloid1.

coeffs = vpasolve([eq1 eq2]);
eqCycloid = subs(solCycloid1,{'C1','C3'},{coeffs.C1,coeffs.C3})

eqCycloid = 
$$-6.4199192418473511250705556729108\hspace{0.17em}\mathrm{atan}\left(\sqrt{\frac{6.4199192418473511250705556729108}{y\left(x\right)}-1}\right)-y\left(x\right)\hspace{0.17em}\sqrt{\frac{6.4199192418473511250705556729108}{y\left(x\right)}-1}=-x-5.8078336827583088482183433150164$$

The implicit equation eqCycloid describes the cycloid solution of the Brachistochrone problem in terms of $$x$$ and $$y(x)$$.

You can then use fimplicit to plot eqCycloid. Since fimplicit only accepts implicit symbolic equations that contain symbolic variables $$x$$ and $$y$$, convert the symbolic function $$y(x)$$ to a symbolic variable $$y$$. Use mapSymType to convert $$y(x)$$ to $$y$$. Plot the cycloid solution within the boundary conditions $$0<x<4$$ and $$1<y<5$$.

funToVar = @(obj) sym('y');
eqPlot = mapSymType(eqCycloid,'symfun',funToVar);
fimplicit(eqPlot,[0 4 1 5])

For a function $$u(x,y)$$ that describes a surface in 3-D space, the surface area can be determined by the functional

where $$u_x$$ and $$u_y$$ are the partial derivatives of $$u$$ with respect to $$x$$ and $$y$$.

Find the functional derivative of the integrand f with respect to u.

syms u(x,y)
f = sqrt(1 + diff(u,x)^2 + diff(u,y)^2);
G = functionalDerivative(f,u)

G(x, y) = 
$$\begin{array}{l}-\frac{{\left(\frac{\partial }{\partial y}\mathrm{ }u\left(x,y\right)\right)}^2\hspace{0.17em}\frac{{\partial }^2}{\partial x^2}\mathrm{ }u\left(x,y\right)+\frac{{\partial }^2}{\partial x^2}\mathrm{ }u\left(x,y\right)+{{\sigma }_1}^2\hspace{0.17em}\frac{{\partial }^2}{\partial y^2}\mathrm{ }u\left(x,y\right)-2\hspace{0.17em}\frac{\partial }{\partial y}\mathrm{ }{\sigma }_1\hspace{0.17em}\frac{\partial }{\partial y}\mathrm{ }u\left(x,y\right)\hspace{0.17em}{\sigma }_1+\frac{{\partial }^2}{\partial y^2}\mathrm{ }u\left(x,y\right)}{{\left({{\sigma }_1}^2+{\left(\frac{\partial }{\partial y}\mathrm{ }u\left(x,y\right)\right)}^2+1\right)}^{3/2}}\\ \\ \mathrm{where}\\ \\ \mathrm{ }{\sigma }_1=\frac{\partial }{\partial x}\mathrm{ }u\left(x,y\right)\end{array}$$

The result is the equation G that describes the minimal surface of a 3-D surface defined by u(x,y). The solutions to this equation describe minimal surfaces in 3-D space, such as soap bubbles.