Multivariable Calculus: Partial Differentiation, Jacobians, Extrema, and Leibniz Rule

Functions of Several Variables and Partial Differentiation

  • Functions of One vs. Two Variables:

    • A function of one variable is written as f(x)f(x), for example: f(x)=3x+2f(x) = 3x + 2.
    • A function of two variables is written as f(x,y)f(x,y), for example: f(x,y)=x2+3ayf(x,y) = x^2 + 3ay.
  • Formal Definitions of Partial Derivatives:

    • The partial derivative of f(x,y)f(x,y) with respect to xx is defined as the limit:     fx=limΔx0f(x+Δx,y)f(x,y)Δx\frac{\partial f}{\partial x} = \lim_{\Delta x \rightarrow 0} \frac{f(x + \Delta x, y) - f(x,y)}{\Delta x}
    • The partial derivative of f(x,y)f(x,y) with respect to yy is defined as the limit:     fy=limΔy0f(x,y+Δy)f(x,y)Δy\frac{\partial f}{\partial y} = \lim_{\Delta y \rightarrow 0} \frac{f(x, y + \Delta y) - f(x,y)}{\Delta y}
  • Example 1: First and Second Order Partial Derivatives:

    • Given the function:     z=x3+y33axyz = x^3 + y^3 - 3axy
    • First-order partial derivative with respect to xx:     zx=3x23ay\frac{\partial z}{\partial x} = 3x^2 - 3ay
    • First-order partial derivative with respect to yy:     zy=3y23ax\frac{\partial z}{\partial y} = 3y^2 - 3ax
    • Second-order partial derivative with respect to xx:     2zx2=6x\frac{\partial^2 z}{\partial x^2} = 6x
    • Second-order partial derivative with respect to yy:     2zy2=6y\frac{\partial^2 z}{\partial y^2} = 6y
    • Mixed partial derivatives:     2zxy=3a\frac{\partial^2 z}{\partial x \partial y} = -3a2zyx=3a\frac{\partial^2 z}{\partial y \partial x} = -3a
    • Equivalence of mixed partial derivatives:     2zxy=2zyx\frac{\partial^2 z}{\partial x \partial y} = \frac{\partial^2 z}{\partial y \partial x}
  • Example 2: Verification of Inverse Tangent Partial Derivatives:

    • Given:     u=tan1(x2+y2xy)u = \tan^{-1}\left(\frac{x^2+y^2}{x-y}\right)
    • Differentiating partially with respect to xx:     ux=11+(x2+y2xy)2×(xy)(2x)(x2+y2)(1)(xy)2=2x22xyx2y2(xy)2+(x2+y2)2\frac{\partial u}{\partial x} = \frac{1}{1 + \left(\frac{x^2+y^2}{x-y}\right)^2} \times \frac{(x-y)(2x) - (x^2+y^2)(1)}{(x-y)^2} = \frac{2x^2 - 2xy - x^2 - y^2}{(x-y)^2 + (x^2+y^2)^2}
    • Differentiating partially with respect to yy:     uy=11+(x2+y2xy)2×(xy)(2y)(x2+y2)(1)(xy)2=2xy2y2+x2+y2(xy)2+(x2+y2)2\frac{\partial u}{\partial y} = \frac{1}{1 + \left(\frac{x^2+y^2}{x-y}\right)^2} \times \frac{(x-y)(2y) - (x^2+y^2)(-1)}{(x-y)^2} = \frac{2xy - 2y^2 + x^2 + y^2}{(x-y)^2 + (x^2+y^2)^2}
  • Example 3: Three-Variable Partial Derivative Identity:

    • Given:     u=(x2+y2+z2)1/2u = (x^2 + y^2 + z^2)^{-1/2}
    • Prove that:     xux+yuy+zuz=ux \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y} + z \frac{\partial u}{\partial z} = -u
    • Proof Steps:
    • Differentiating uu partially with respect to xx:       ux=12(x2+y2+z2)3/2(2x)=x(x2+y2+z2)3/2\frac{\partial u}{\partial x} = -\frac{1}{2}(x^2 + y^2 + z^2)^{-3/2} (2x) = -x(x^2 + y^2 + z^2)^{-3/2}
    • Multiplying by xx:       xux=x2(x2+y2+z2)3/2x \frac{\partial u}{\partial x} = -x^2 (x^2 + y^2 + z^2)^{-3/2}
    • Similarly, for yy and zz:       yuy=y2(x2+y2+z2)3/2y \frac{\partial u}{\partial y} = -y^2 (x^2 + y^2 + z^2)^{-3/2}zuz=z2(x2+y2+z2)3/2z \frac{\partial u}{\partial z} = -z^2 (x^2 + y^2 + z^2)^{-3/2}
    • Adding the three expressions:       xux+yuy+zuz=(x2+y2+z2)(x2+y2+z2)3/2x \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y} + z \frac{\partial u}{\partial z} = -(x^2 + y^2 + z^2)(x^2 + y^2 + z^2)^{-3/2}=(x2+y2+z2)1/2=u= -(x^2 + y^2 + z^2)^{-1/2} = -u
  • Practice Problems on Partial Differentiation:

    • Problem 1: If u=x2(yz)+y2(zx)+z2(xy)u = x^2(y-z) + y^2(z-x) + z^2(x-y), show that:     ux+uy+uz=0\frac{\partial u}{\partial x} + \frac{\partial u}{\partial y} + \frac{\partial u}{\partial z} = 0
    • Problem 2: If V=(12xy+y2)1/2V = (1 - 2xy + y^2)^{-1/2}, prove that:     xVxyVy=y2V3x \frac{\partial V}{\partial x} - y \frac{\partial V}{\partial y} = y^2 V^3
    • Problem 3: If u=log(x3+y3+z33xyz)u = \log(x^3 + y^3 + z^3 - 3xyz), show that:
    • (i) ux+uy+uz=3x+y+z\frac{\partial u}{\partial x} + \frac{\partial u}{\partial y} + \frac{\partial u}{\partial z} = \frac{3}{x+y+z}
      • Proof:         ux=3x23yzx3+y3+z33xyz\frac{\partial u}{\partial x} = \frac{3x^2 - 3yz}{x^3 + y^3 + z^3 - 3xyz}uy=3y23xzx3+y3+z33xyz\frac{\partial u}{\partial y} = \frac{3y^2 - 3xz}{x^3 + y^3 + z^3 - 3xyz}uz=3z23xyx3+y3+z33xyz\frac{\partial u}{\partial z} = \frac{3z^2 - 3xy}{x^3 + y^3 + z^3 - 3xyz}ux+uy+uz=3(x2+y2+z2xyyzzx)(x+y+z)(x2+y2+z2xyyzzx)=3x+y+z\frac{\partial u}{\partial x} + \frac{\partial u}{\partial y} + \frac{\partial u}{\partial z} = \frac{3(x^2 + y^2 + z^2 - xy - yz - zx)}{(x+y+z)(x^2 + y^2 + z^2 - xy - yz - zx)} = \frac{3}{x+y+z}
    • (ii) (x+y+z)2u=9(x+y+z)2\left(\frac{\partial}{\partial x} + \frac{\partial}{\partial y} + \frac{\partial}{\partial z}\right)^2 u = \frac{-9}{(x+y+z)^2}
      • Proof:         (x+y+z)2u=(x+y+z)(3x+y+z)\left(\frac{\partial}{\partial x} + \frac{\partial}{\partial y} + \frac{\partial}{\partial z}\right)^2 u = \left(\frac{\partial}{\partial x} + \frac{\partial}{\partial y} + \frac{\partial}{\partial z}\right) \left(\frac{3}{x+y+z}\right)=x(3x+y+z)+y(3x+y+z)+z(3x+y+z)= \frac{\partial}{\partial x}\left(\frac{3}{x+y+z}\right) + \frac{\partial}{\partial y}\left(\frac{3}{x+y+z}\right) + \frac{\partial}{\partial z}\left(\frac{3}{x+y+z}\right)=3(x+y+z)2+3(x+y+z)2+3(x+y+z)2=9(x+y+z)2= \frac{-3}{(x+y+z)^2} + \frac{-3}{(x+y+z)^2} + \frac{-3}{(x+y+z)^2} = \frac{-9}{(x+y+z)^2}
  • Example 4: Heat/Diffusion Differential Equation:

    • Given:     θ=tner24t\theta = t^n e^{-\frac{r^2}{4t}}
    • Find the value of nn which satisfies:     1r2r(r2θr)=θt\frac{1}{r^2} \frac{\partial}{\partial r}\left(r^2 \frac{\partial \theta}{\partial r}\right) = \frac{\partial \theta}{\partial t}
    • Solution Steps:
    • Compute first derivative with respect to rr:       θr=tner24t(2r4t)=r2tn1er24t\frac{\partial \theta}{\partial r} = t^n e^{-\frac{r^2}{4t}} \left(-\frac{2r}{4t}\right) = -\frac{r}{2} t^{n-1} e^{-\frac{r^2}{4t}}
    • Multiply by r2r^2:       r2θr=r32tn1er24tr^2 \frac{\partial \theta}{\partial r} = -\frac{r^3}{2} t^{n-1} e^{-\frac{r^2}{4t}}
    • Compute derivative of r2θrr^2 \frac{\partial \theta}{\partial r} with respect to rr:       r(r2θr)=3r22tn1er24tr32tn1er24t(2r4t)\frac{\partial}{\partial r}\left(r^2 \frac{\partial \theta}{\partial r}\right) = -\frac{3r^2}{2} t^{n-1} e^{-\frac{r^2}{4t}} - \frac{r^3}{2} t^{n-1} e^{-\frac{r^2}{4t}} \left(-\frac{2r}{4t}\right)=(3r22tn1+r44tn2)er24t= \left(-\frac{3r^2}{2} t^{n-1} + \frac{r^4}{4} t^{n-2}\right) e^{-\frac{r^2}{4t}}
    • Divide by r2r^2:       1r2r(r2θr)=(32tn1+r24tn2)er24t\frac{1}{r^2} \frac{\partial}{\partial r}\left(r^2 \frac{\partial \theta}{\partial r}\right) = \left(-\frac{3}{2} t^{n-1} + \frac{r^2}{4} t^{n-2}\right) e^{-\frac{r^2}{4t}}
    • Compute derivative of θ\theta with respect to tt:       θt=ntn1er24t+tner24t(r24t2)=(ntn1+r24tn2)er24t\frac{\partial \theta}{\partial t} = n t^{n-1} e^{-\frac{r^2}{4t}} + t^n e^{-\frac{r^2}{4t}} \left(\frac{r^2}{4t^2}\right) = \left(n t^{n-1} + \frac{r^2}{4} t^{n-2}\right) e^{-\frac{r^2}{4t}}
    • Equating the two expressions:       ntn1=32tn1    n=32n t^{n-1} = -\frac{3}{2} t^{n-1} \implies n = -\frac{3}{2}
  • Example 5: One-Dimensional Wave Equation Solution:

    • If θ=ϕ(y+ax)+ψ(yax)\theta = \phi(y + ax) + \psi(y - ax), show that:     2θx2=a22θy2\frac{\partial^2 \theta}{\partial x^2} = a^2 \frac{\partial^2 \theta}{\partial y^2}
    • Proof:     θx=aϕ(y+ax)aψ(yax)\frac{\partial \theta}{\partial x} = a \phi'(y + ax) - a \psi'(y - ax)2θx2=a2ϕ(y+ax)+a2ψ(yax)=a2[ϕ(y+ax)+ψ(yax)]\frac{\partial^2 \theta}{\partial x^2} = a^2 \phi''(y + ax) + a^2 \psi''(y - ax) = a^2 [\phi''(y + ax) + \psi''(y - ax)]θy=ϕ(y+ax)+ψ(yax)\frac{\partial \theta}{\partial y} = \phi'(y + ax) + \psi'(y - ax)2θy2=ϕ(y+ax)+ψ(yax)\frac{\partial^2 \theta}{\partial y^2} = \phi''(y + ax) + \psi''(y - ax)2θx2=a22θy2\frac{\partial^2 \theta}{\partial x^2} = a^2 \frac{\partial^2 \theta}{\partial y^2}

Euler's Theorem on Homogeneous Functions

  • Definition of a Homogeneous Function:

    • A function f(x,y)f(x,y) is said to be homogeneous of degree nn in xx and yy if it can be expressed in the form:     z=xnf(yx)z = x^n f\left(\frac{y}{x}\right)
    • Example: For z=x3y3x+yz = \frac{x^3 - y^3}{x + y}, factoring x3x^3 from numerator and xx from denominator:     z=x3(1y3x3)x(1+yx)=x2f(yx)z = \frac{x^3 \left(1 - \frac{y^3}{x^3}\right)}{x \left(1 + \frac{y}{x}\right)} = x^2 f\left(\frac{y}{x}\right)     Therefore, zz is a homogeneous function of degree n=2n = 2
  • Euler's Theorem Statement:

    • If z=f(x,y)z = f(x,y) is a homogeneous function of xx and yy of order nn, then:     xzx+yzy=nzx \frac{\partial z}{\partial x} + y \frac{\partial z}{\partial y} = n z     for all (x,y)(x,y) in the domain of the function.
  • Complete Proof of Euler's Theorem:

    • Since zz is a homogeneous function of xx and yy of order nn:     z=xnf(yx)z = x^n f\left(\frac{y}{x}\right)
    • Differentiating partially with respect to xx:     zx=nxn1f(yx)+xnf(yx)(yx2)=nxn1f(yx)yxn2f(yx)\frac{\partial z}{\partial x} = n x^{n-1} f\left(\frac{y}{x}\right) + x^n f'\left(\frac{y}{x}\right) \left(-\frac{y}{x^2}\right) = n x^{n-1} f\left(\frac{y}{x}\right) - y x^{n-2} f'\left(\frac{y}{x}\right)
    • Multiplying both sides by xx:     xzx=nxnf(yx)yxn1f(yx)x \frac{\partial z}{\partial x} = n x^n f\left(\frac{y}{x}\right) - y x^{n-1} f'\left(\frac{y}{x}\right)xzx=nzyxn1f(yx)x \frac{\partial z}{\partial x} = n z - y x^{n-1} f'\left(\frac{y}{x}\right)
    • Differentiating zz partially with respect to yy:     zy=xnf(yx)(1x)=xn1f(yx)\frac{\partial z}{\partial y} = x^n f'\left(\frac{y}{x}\right) \left(\frac{1}{x}\right) = x^{n-1} f'\left(\frac{y}{x}\right)
    • Multiplying both sides by yy:     yzy=yxn1f(yx)y \frac{\partial z}{\partial y} = y x^{n-1} f'\left(\frac{y}{x}\right)
    • Adding the equations for xzxx \frac{\partial z}{\partial x} and yzyy \frac{\partial z}{\partial y}:     xzx+yzy=nzyxn1f(yx)+yxn1f(yx)=nzx \frac{\partial z}{\partial x} + y \frac{\partial z}{\partial y} = n z - y x^{n-1} f'\left(\frac{y}{x}\right) + y x^{n-1} f'\left(\frac{y}{x}\right) = n z
  • Second-Order Extension of Euler's Theorem:

    • Statement: If uu is a homogeneous function of degree nn, then:     x22ux2+2xy2uxy+y22uy2=n(n1)ux^2 \frac{\partial^2 u}{\partial x^2} + 2xy \frac{\partial^2 u}{\partial x \partial y} + y^2 \frac{\partial^2 u}{\partial y^2} = n(n-1) u
    • Proof:
    • From Euler's theorem:       xux+yuy=nux \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y} = n u
    • Differentiating partially with respect to xx:       ux+x2ux2+y2uxy=nux\frac{\partial u}{\partial x} + x \frac{\partial^2 u}{\partial x^2} + y \frac{\partial^2 u}{\partial x \partial y} = n \frac{\partial u}{\partial x}x2ux2+y2uxy=(n1)uxx \frac{\partial^2 u}{\partial x^2} + y \frac{\partial^2 u}{\partial x \partial y} = (n-1) \frac{\partial u}{\partial x}
    • Differentiating partially with respect to yy:       x2uyx+uy+y2uy2=nuyx \frac{\partial^2 u}{\partial y \partial x} + \frac{\partial u}{\partial y} + y \frac{\partial^2 u}{\partial y^2} = n \frac{\partial u}{\partial y}x2uxy+y2uy2=(n1)uyx \frac{\partial^2 u}{\partial x \partial y} + y \frac{\partial^2 u}{\partial y^2} = (n-1) \frac{\partial u}{\partial y}
    • Multiplying the first derivative equation by xx and the second by yy, and adding:       x22ux2+2xy2uxy+y22uy2=(n1)(xux+yuy)=n(n1)ux^2 \frac{\partial^2 u}{\partial x^2} + 2xy \frac{\partial^2 u}{\partial x \partial y} + y^2 \frac{\partial^2 u}{\partial y^2} = (n-1)\left(x \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y}\right) = n(n-1)u
  • Applications & Worked Examples of Euler's Theorem:

    • Example 1: If u=sin1(x2+y2x+y)u = \sin^{-1}\left(\frac{x^2+y^2}{x+y}\right), show that xux+yuy=tan(u)x \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y} = \tan(u).
    • Solution:       Let z=sin(u)=x2+y2x+y=x2(1+y2x2)x(1+yx)=x1f(yx)z = \sin(u) = \frac{x^2+y^2}{x+y} = \frac{x^2\left(1 + \frac{y^2}{x^2}\right)}{x\left(1 + \frac{y}{x}\right)} = x^1 f\left(\frac{y}{x}\right).       zz is a homogeneous function of degree n=1n = 1.       By Euler's Theorem:       xzx+yzy=1×zx \frac{\partial z}{\partial x} + y \frac{\partial z}{\partial y} = 1 \times z       Since z=sin(u)z = \sin(u), zx=cos(u)ux\frac{\partial z}{\partial x} = \cos(u) \frac{\partial u}{\partial x} and zy=cos(u)uy\frac{\partial z}{\partial y} = \cos(u) \frac{\partial u}{\partial y}.       xcos(u)ux+ycos(u)uy=sin(u)x \cos(u) \frac{\partial u}{\partial x} + y \cos(u) \frac{\partial u}{\partial y} = \sin(u)xux+yuy=sin(u)cos(u)=tan(u)x \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y} = \frac{\sin(u)}{\cos(u)} = \tan(u)
    • Example 2: If u=log(x2+y2x+y)u = \log\left(\frac{x^2+y^2}{x+y}\right), show that xux+yuy=1x \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y} = 1
    • Solution:       Let z=eu=x2+y2x+yz = e^u = \frac{x^2+y^2}{x+y}, which is homogeneous of degree n=1n = 1.       By Euler's Theorem:       xzx+yzy=1×zx \frac{\partial z}{\partial x} + y \frac{\partial z}{\partial y} = 1 \times zxeuux+yeuuy=eux e^u \frac{\partial u}{\partial x} + y e^u \frac{\partial u}{\partial y} = e^uxux+yuy=1x \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y} = 1
  • Practice Exercises on Euler's Theorem:

    • Problem 1: If u=log(x5+y5+z5x2+y2+z2)u = \log\left(\frac{x^5 + y^5 + z^5}{x^2 + y^2 + z^2}\right), show that:     xux+yuy+zuz=3x \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y} + z \frac{\partial u}{\partial z} = 3
    • Problem 2: If u=sin1(x3+y3+z3ax+by+cz)u = \sin^{-1}\left(\frac{x^3 + y^3 + z^3}{ax + by + cz}\right), prove that:     xux+yuy+zuz=2tan(u)x \frac{\partial u}{\partial x} + y \frac{\partial u}{\partial y} + z \frac{\partial u}{\partial z} = 2 \tan(u)
    • Problem 3: If u=f(yx)+g(xy)u = f\left(\frac{y}{x}\right) + g\left(\frac{x}{y}\right), prove that:     x22ux2+2xy2uxy+y22uy2=0x^2 \frac{\partial^2 u}{\partial x^2} + 2xy \frac{\partial^2 u}{\partial x \partial y} + y^2 \frac{\partial^2 u}{\partial y^2} = 0
    • Problem 4: If u=tan1(yx)u = \tan^{-1}\left(\frac{y}{x}\right), prove that:     x22ux2+2xy2uxy+y22uy2=sin(2u)sin2(u)x^2 \frac{\partial^2 u}{\partial x^2} + 2xy \frac{\partial^2 u}{\partial x \partial y} + y^2 \frac{\partial^2 u}{\partial y^2} = -\sin(2u) \sin^2(u)

Total Differentials, Chain Rule, and Composite Functions

  • Definitions & Concepts:

    • Let z=f(x,y)z = f(x,y) be a function of variables xx and yy, where xx and yy are themselves functions of another variable tt (x=h(t)x = h(t), y=g(t)y = g(t)). Then zz is called a composite function of single variable tt.
    • Theorem on Total Differentials: If z=f(x,y)z = f(x,y) possesses continuous first-order partial derivatives, the total differential of zz is given by:     dz=fxdx+fydydz = \frac{\partial f}{\partial x} dx + \frac{\partial f}{\partial y} dy
    • Total Derivative (Chain Rule for 1 Variable tt):     dzdt=fxdxdt+fydydt\frac{dz}{dt} = \frac{\partial f}{\partial x} \frac{dx}{dt} + \frac{\partial f}{\partial y} \frac{dy}{dt}
  • Example 1: Total Differential Calculation:

    • Find total differential of f(x,y)=yex+yf(x,y) = y e^{x+y}.     df=fxdx+fydydf = \frac{\partial f}{\partial x} dx + \frac{\partial f}{\partial y} dyfx=yex+y\frac{\partial f}{\partial x} = y e^{x+y}fy=yex+y+ex+y=(y+1)ex+y\frac{\partial f}{\partial y} = y e^{x+y} + e^{x+y} = (y+1) e^{x+y}df=yex+ydx+(y+1)ex+ydydf = y e^{x+y} dx + (y+1) e^{x+y} dy
  • Example 2: Total Derivative Calculation:

    • Find total derivative of f(x,y)=x2+3xyf(x,y) = x^2 + 3xy with respect to uu, given x=sin(u)x = \sin(u), y=cos(u)y = \cos(u).     dfdu=fxdxdu+fydydu\frac{df}{du} = \frac{\partial f}{\partial x} \frac{dx}{du} + \frac{\partial f}{\partial y} \frac{dy}{du}fx=2x+3y,fy=3x\frac{\partial f}{\partial x} = 2x + 3y, \quad \frac{\partial f}{\partial y} = 3xdxdu=cos(u),dydu=sin(u)\frac{dx}{du} = \cos(u), \quad \frac{dy}{du} = -\sin(u)dfdu=(2x+3y)cos(u)+(3x)(sin(u))=(2sin(u)+3cos(u))cos(u)3sin(u)sin(u)\frac{df}{du} = (2x + 3y) \cos(u) + (3x) (-\sin(u)) = (2 \sin(u) + 3 \cos(u)) \cos(u) - 3 \sin(u) \sin(u)
  • Change of Variables (Chain Rule for Two Independent Variables):

    • Let u=f(x,y)u = f(x,y) where x=ϕ(s,t)x = \phi(s,t) and y=ψ(s,t)y = \psi(s,t). Then:     us=uxxs+uyys\frac{\partial u}{\partial s} = \frac{\partial u}{\partial x} \frac{\partial x}{\partial s} + \frac{\partial u}{\partial y} \frac{\partial y}{\partial s}ut=uxxt+uyyt\frac{\partial u}{\partial t} = \frac{\partial u}{\partial x} \frac{\partial x}{\partial t} + \frac{\partial u}{\partial y} \frac{\partial y}{\partial t}
  • Example 3: Partial Derivatives in Polar Coordinates:

    • Given z=x3xy+y3z = x^3 - xy + y^3 where x=rcos(θ)x = r \cos(\theta) and y=rsin(θ)y = r \sin(\theta), find zr\frac{\partial z}{\partial r} and zθ\frac{\partial z}{\partial \theta}.
    • Solution:     zx=3x2y,zy=x+3y2\frac{\partial z}{\partial x} = 3x^2 - y, \quad \frac{\partial z}{\partial y} = -x + 3y^2xr=cos(θ),yr=sin(θ)\frac{\partial x}{\partial r} = \cos(\theta), \quad \frac{\partial y}{\partial r} = \sin(\theta)zr=(3x2y)cos(θ)+(x+3y2)sin(θ)\frac{\partial z}{\partial r} = (3x^2 - y) \cos(\theta) + (-x + 3y^2) \sin(\theta)     Substitute x=rcos(θ)x = r \cos(\theta) and y=rsin(θ)y = r \sin(\theta):     =[3(rcos(θ))2rsin(θ)]cos(θ)+[rcos(θ)+3(rsin(θ))2]sin(θ)= [3(r \cos(\theta))^2 - r \sin(\theta)] \cos(\theta) + [-r \cos(\theta) + 3(r \sin(\theta))^2] \sin(\theta)=3r2(cos3(θ)+sin3(θ))2rsin(θ)cos(θ)= 3r^2 (\cos^3(\theta) + \sin^3(\theta)) - 2r \sin(\theta) \cos(\theta)zθ=zx(rsin(θ))+zy(rcos(θ))\frac{\partial z}{\partial \theta} = \frac{\partial z}{\partial x} (-r \sin(\theta)) + \frac{\partial z}{\partial y} (r \cos(\theta))=(3x2y)(rsin(θ))+(x+3y2)(rcos(θ))= (3x^2 - y) (-r \sin(\theta)) + (-x + 3y^2) (r \cos(\theta))=r[3r2sin2(θ)cos(θ)3r2cos2(θ)sin(θ)+rcos(θ)sin(θ)rsin(θ)cos(θ)]= r [3r^2 \sin^2(\theta) \cos(\theta) - 3r^2 \cos^2(\theta) \sin(\theta) + r \cos(\theta) \sin(\theta) - r \sin(\theta) \cos(\theta)]=r2(sin(θ)cos(θ))[3rsin(θ)cos(θ)+1]= r^2 (\sin(\theta) - \cos(\theta)) [3r \sin(\theta) \cos(\theta) + 1]
  • Problem 4: Proof for Composite Exponential Function:

    • If z=eax+byf(axby)z = e^{ax+by} f(ax-by), show that:     bzx+azy=2abzb \frac{\partial z}{\partial x} + a \frac{\partial z}{\partial y} = 2ab z
    • Proof:     Let u=ax+byu = ax + by and v=axbyv = ax - by. Then z=euf(v)z = e^u f(v).     zx=zuux+zvvx=[euf(v)]a+[euf(v)]a=aeu[f(v)+f(v)]\frac{\partial z}{\partial x} = \frac{\partial z}{\partial u} \frac{\partial u}{\partial x} + \frac{\partial z}{\partial v} \frac{\partial v}{\partial x} = [e^u f(v)] a + [e^u f'(v)] a = a e^u [f(v) + f'(v)]zy=zuuy+zvvy=[euf(v)]b+[euf(v)](b)=beu[f(v)f(v)]\frac{\partial z}{\partial y} = \frac{\partial z}{\partial u} \frac{\partial u}{\partial y} + \frac{\partial z}{\partial v} \frac{\partial v}{\partial y} = [e^u f(v)] b + [e^u f'(v)] (-b) = b e^u [f(v) - f'(v)]bzx+azy=abeu[f(v)+f(v)]+abeu[f(v)f(v)]=2abeuf(v)=2abzb \frac{\partial z}{\partial x} + a \frac{\partial z}{\partial y} = ab e^u [f(v) + f'(v)] + ab e^u [f(v) - f'(v)] = 2ab e^u f(v) = 2ab z
  • Problem 5: Proof for Cyclic Variable Function:

    • If u=f(yz,zx,xy)u = f(y-z, z-x, x-y), prove that:     ux+uy+uz=0\frac{\partial u}{\partial x} + \frac{\partial u}{\partial y} + \frac{\partial u}{\partial z} = 0
    • Proof:     Let X=yzX = y-z, Y=zxY = z-x, Z=xyZ = x-y.     ux=fX(0)+fY(1)+fZ(1)=fY+fZ\frac{\partial u}{\partial x} = \frac{\partial f}{\partial X}(0) + \frac{\partial f}{\partial Y}(-1) + \frac{\partial f}{\partial Z}(1) = -\frac{\partial f}{\partial Y} + \frac{\partial f}{\partial Z}uy=fX(1)+fY(0)+fZ(1)=fXfZ\frac{\partial u}{\partial y} = \frac{\partial f}{\partial X}(1) + \frac{\partial f}{\partial Y}(0) + \frac{\partial f}{\partial Z}(-1) = \frac{\partial f}{\partial X} - \frac{\partial f}{\partial Z}uz=fX(1)+fY(1)+fZ(0)=fX+fY\frac{\partial u}{\partial z} = \frac{\partial f}{\partial X}(-1) + \frac{\partial f}{\partial Y}(1) + \frac{\partial f}{\partial Z}(0) = -\frac{\partial f}{\partial X} + \frac{\partial f}{\partial Y}     Adding all three equations yields zero.

Transformation of Coordinates and Differential Operators

  • Transformation of Laplace Equation into Polar Coordinates:
    • Transform the equation:     2ux2+2uy2=0\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = 0     into polar coordinates (r,θ)(r, \theta).
    • Coordinate Relations:     x=rcos(θ),y=rsin(θ)x = r \cos(\theta), \quad y = r \sin(\theta)r2=x2+y2    r=x2+y2r^2 = x^2 + y^2 \implies r = \sqrt{x^2+y^2}tan(θ)=yx    θ=tan1(yx)\tan(\theta) = \frac{y}{x} \implies \theta = \tan^{-1}\left(\frac{y}{x}\right)
    • Derivative Relations:     rx=xr=cos(θ),ry=yr=sin(θ)\frac{\partial r}{\partial x} = \frac{x}{r} = \cos(\theta), \quad \frac{\partial r}{\partial y} = \frac{y}{r} = \sin(\theta)θx=yr2=sin(θ)r,θy=xr2=cos(θ)r\frac{\partial \theta}{\partial x} = -\frac{y}{r^2} = -\frac{\sin(\theta)}{r}, \quad \frac{\partial \theta}{\partial y} = \frac{x}{r^2} = \frac{\cos(\theta)}{r}
    • First Derivatives:     ux=cos(θ)ursin(θ)ruθ\frac{\partial u}{\partial x} = \cos(\theta) \frac{\partial u}{\partial r} - \frac{\sin(\theta)}{r} \frac{\partial u}{\partial \theta}uy=sin(θ)ur+cos(θ)ruθ\frac{\partial u}{\partial y} = \sin(\theta) \frac{\partial u}{\partial r} + \frac{\cos(\theta)}{r} \frac{\partial u}{\partial \theta}
    • Operator Equivalences:     x=cos(θ)rsin(θ)rθ\frac{\partial}{\partial x} = \cos(\theta) \frac{\partial}{\partial r} - \frac{\sin(\theta)}{r} \frac{\partial}{\partial \theta}y=sin(θ)r+cos(θ)rθ\frac{\partial}{\partial y} = \sin(\theta) \frac{\partial}{\partial r} + \frac{\cos(\theta)}{r} \frac{\partial}{\partial \theta}
    • Second Derivatives:     2ux2=(cos(θ)rsin(θ)rθ)(cos(θ)ursin(θ)ruθ)\frac{\partial^2 u}{\partial x^2} = \left(\cos(\theta) \frac{\partial}{\partial r} - \frac{\sin(\theta)}{r} \frac{\partial}{\partial \theta}\right) \left(\cos(\theta) \frac{\partial u}{\partial r} - \frac{\sin(\theta)}{r} \frac{\partial u}{\partial \theta}\right)2uy2=(sin(θ)r+cos(θ)rθ)(sin(θ)ur+cos(θ)ruθ)\frac{\partial^2 u}{\partial y^2} = \left(\sin(\theta) \frac{\partial}{\partial r} + \frac{\cos(\theta)}{r} \frac{\partial}{\partial \theta}\right) \left(\sin(\theta) \frac{\partial u}{\partial r} + \frac{\cos(\theta)}{r} \frac{\partial u}{\partial \theta}\right)
    • Summing the two expressions yields the polar form of the Laplacian:     2ux2+2uy2=2ur2+1rur+1r22uθ2=0\frac{\partial^2 u}{\partial x^2} + \frac{\partial^2 u}{\partial y^2} = \frac{\partial^2 u}{\partial r^2} + \frac{1}{r} \frac{\partial u}{\partial r} + \frac{1}{r^2} \frac{\partial^2 u}{\partial \theta^2} = 0

Differentiation of Implicit Functions

Implicit function differentiation formulas and example

  • Definitions:

    • Implicit Function: A relation of the form f(x,y)=cf(x,y) = c in which xx and yy cannot be separated out explicitly.
    • Example: The relation x3+y33axy=0x^3 + y^3 - 3axy = 0 is an implicit function.
    • Explicit Function: A relation in which one variable can be isolated completely in terms of the other.
    • Example: The relation y2=x2+2xy^2 = x^2 + 2x is an explicit function.
  • Formulas for Differentiation of Implicit Functions:

    • Let f(x,y)=cf(x,y) = c be an implicit function. Then:
    • (i) First derivative formula:       dydx=fxfy,fy0\frac{dy}{dx} = -\frac{f_x}{f_y}, \quad f_y \neq 0
    • (ii) Second derivative formula:       d2ydx2=fxx(fy)22fxfyfxy+fyy(fx)2(fy)3,fy0provided fxy=fyx\frac{d^2 y}{dx^2} = -\frac{f_{xx}(f_y)^2 - 2 f_x f_y f_{xy} + f_{yy}(f_x)^2}{(f_y)^3}, \quad f_y \neq 0 \quad \text{provided } f_{xy} = f_{yx}
  • Example 1: Exponential Variable Power:

    • Find dydx\frac{dy}{dx} if xy+yx=cx^y + y^x = c
    • Solution:     Let f(x,y)=xy+yxc=0f(x,y) = x^y + y^x - c = 0     Partial derivative with respect to xx:     fx=yxy1+yxlog(y)f_x = y x^{y-1} + y^x \log(y)     Partial derivative with respect to yy:     fy=xylog(x)+xyx1f_y = x^y \log(x) + x y^{x-1}     Applying formula dydx=fxfy\frac{dy}{dx} = -\frac{f_x}{f_y}:     dydx=yxy1+yxlog(y)xylog(x)+xyx1\frac{dy}{dx} = -\frac{y x^{y-1} + y^x \log(y)}{x^y \log(x) + x y^{x-1}}
  • Example 2: Folium of Descartes Derivative:

    • If x3+y33axy=0x^3 + y^3 - 3axy = 0, find dydx\frac{dy}{dx} at x=a,y=ax = a, y = a:     fx=3x23ay,fy=3y23axf_x = 3x^2 - 3ay, \quad f_y = 3y^2 - 3axdydx=3x23ay3y23ax=x2ayy2ax\frac{dy}{dx} = -\frac{3x^2 - 3ay}{3y^2 - 3ax} = -\frac{x^2 - ay}{y^2 - ax}

Exact and Inexact Differentials

  • Definitions:

    • Exact Differential: Differentials of the form:     M(x,y)dx+N(x,y)dy=d(f(x,y))M(x,y) dx + N(x,y) dy = d(f(x,y))     which can be integrated directly, are called exact differentials.
    • Inexact Differential: Differentials that cannot be written directly as an exact differential of a single function without an integrating factor.
  • Condition for Exactness:

    • A differential Mdx+NdyM dx + N dy is exact if and only if:     My=Nx\frac{\partial M}{\partial y} = \frac{\partial N}{\partial x}
  • Examples:

    • Example of Exact Differential: xdy+ydx=d(xy)    f(x,y)=xy+cx dy + y dx = d(xy) \implies f(x,y) = xy + c.     Here M=yM = y and N=xN = x, so My=1\frac{\partial M}{\partial y} = 1 and Nx=1\frac{\partial N}{\partial x} = 1.
    • Example of Inexact Differential: Show that xdy+3ydxx dy + 3y dx is inexact.     Here M=3yM = 3y and N=xN = x. My=3\frac{\partial M}{\partial y} = 3 whereas Nx=1\frac{\partial N}{\partial x} = 1. Since 313 \neq 1, the differential is inexact.

Jacobians (Functional Determinants)

  • Definition of Jacobian for Two Variables:

    • If uu and vv are functions of two independent variables xx and yy, then the Jacobian of u,vu, v with respect to x,yx, y is defined as the determinant:     J=(u,v)(x,y)=J(u,vx,y)=uxuyvxvyJ = \frac{\partial(u,v)}{\partial(x,y)} = J\left(\frac{u,v}{x,y}\right) = \begin{vmatrix} \frac{\partial u}{\partial x} & \frac{\partial u}{\partial y} \\ \frac{\partial v}{\partial x} & \frac{\partial v}{\partial y} \end{vmatrix}
  • Definition of Jacobian for Three Variables:

    • Similarly, the Jacobian of u,v,wu, v, w with respect to x,y,zx, y, z is:     (u,v,w)(x,y,z)=uxuyuzvxvyvzwxwywz\frac{\partial(u,v,w)}{\partial(x,y,z)} = \begin{vmatrix} \frac{\partial u}{\partial x} & \frac{\partial u}{\partial y} & \frac{\partial u}{\partial z} \\ \frac{\partial v}{\partial x} & \frac{\partial v}{\partial y} & \frac{\partial v}{\partial z} \\ \frac{\partial w}{\partial x} & \frac{\partial w}{\partial y} & \frac{\partial w}{\partial z} \end{vmatrix}
  • Property 1: Reciprocal Property of Jacobians:

    • Statement: If J=(u,v)(x,y)J = \frac{\partial(u,v)}{\partial(x,y)} and J=(x,y)(u,v)J' = \frac{\partial(x,y)}{\partial(u,v)}, then:     JJ=1J J' = 1
    • Proof:     Let u=f(x,y)u = f(x,y) and v=g(x,y)v = g(x,y). Solving for xx and yy gives x=ϕ(u,v)x = \phi(u,v) and y=ψ(u,v)y = \psi(u,v).     Total differentials:     du=uxdx+uydydu = \frac{\partial u}{\partial x} dx + \frac{\partial u}{\partial y} dydv=vxdx+vydydv = \frac{\partial v}{\partial x} dx + \frac{\partial v}{\partial y} dy     Substituting dx=xudu+xvdvdx = \frac{\partial x}{\partial u} du + \frac{\partial x}{\partial v} dv and dy=yudu+yvdvdy = \frac{\partial y}{\partial u} du + \frac{\partial y}{\partial v} dv yields identity matrix product of partial derivatives:     (uxuyvxvy)(xuxvyuyv)=(1001)\begin{pmatrix} \frac{\partial u}{\partial x} & \frac{\partial u}{\partial y} \\ \frac{\partial v}{\partial x} & \frac{\partial v}{\partial y} \end{pmatrix} \begin{pmatrix} \frac{\partial x}{\partial u} & \frac{\partial x}{\partial v} \\ \frac{\partial y}{\partial u} & \frac{\partial y}{\partial v} \end{pmatrix} = \begin{pmatrix} 1 & 0 \\ 0 & 1 \end{pmatrix}     Taking determinants on both sides gives JJ=1J J' = 1
  • Property 2: Chain Rule for Jacobians:

    • If u,vu, v are functions of r,sr, s, and r,sr, s are functions of x,yx, y, then:     (u,v)(x,y)=(u,v)(r,s)×(r,s)(x,y)\frac{\partial(u,v)}{\partial(x,y)} = \frac{\partial(u,v)}{\partial(r,s)} \times \frac{\partial(r,s)}{\partial(x,y)}
  • Jacobian of Implicit Functions:

    • If u1,u2u_1, u_2 are functions of x,yx, y defined implicitly by two equations:     f1(u1,u2,x,y)=0f_1(u_1, u_2, x, y) = 0f2(u1,u2,x,y)=0f_2(u_1, u_2, x, y) = 0     Then:     (u1,u2)(x,y)=(1)2(f1,f2)(x,y)(f1,f2)(u1,u2)\frac{\partial(u_1, u_2)}{\partial(x,y)} = (-1)^2 \frac{\frac{\partial(f_1, f_2)}{\partial(x,y)}}{\frac{\partial(f_1, f_2)}{\partial(u_1, u_2)}}
    • Similarly, if u1,u2,u3u_1, u_2, u_3 are functions of x,y,zx, y, z defined implicitly by three equations f1=0,f2=0,f3=0f_1 = 0, f_2 = 0, f_3 = 0:     (u1,u2,u3)(x,y,z)=(1)3(f1,f2,f3)(x,y,z)(f1,f2,f3)(u1,u2,u3)\frac{\partial(u_1, u_2, u_3)}{\partial(x,y,z)} = (-1)^3 \frac{\frac{\partial(f_1, f_2, f_3)}{\partial(x,y,z)}}{\frac{\partial(f_1, f_2, f_3)}{\partial(u_1, u_2, u_3)}}
  • Worked Examples on Jacobians:

    • Example 1: Given u=x2+xyu = x^2 + xy and v=xyv = xy, find (u,v)(x,y)\frac{\partial(u,v)}{\partial(x,y)}.     ux=2x+y,uy=x\frac{\partial u}{\partial x} = 2x + y, \quad \frac{\partial u}{\partial y} = xvx=y,vy=x\frac{\partial v}{\partial x} = y, \quad \frac{\partial v}{\partial y} = x(u,v)(x,y)=2x+yxyx=x(2x+y)xy=2x2+xyxy=2x2\frac{\partial(u,v)}{\partial(x,y)} = \begin{vmatrix} 2x+y & x \\ y & x \end{vmatrix} = x(2x+y) - xy = 2x^2 + xy - xy = 2x^2
    • Example 2: Verification of JJ=1J J' = 1 for Polar Coordinates:     x=rcos(θ),y=rsin(θ)x = r \cos(\theta), \quad y = r \sin(\theta)J=(x,y)(r,θ)=cos(θ)rsin(θ)sin(θ)rcos(θ)=rcos2(θ)+rsin2(θ)=rJ = \frac{\partial(x,y)}{\partial(r,\theta)} = \begin{vmatrix} \cos(\theta) & -r \sin(\theta) \\ \sin(\theta) & r \cos(\theta) \end{vmatrix} = r \cos^2(\theta) + r \sin^2(\theta) = r     For inverse functions r=x2+y2r = \sqrt{x^2+y^2} and θ=tan1(yx)\theta = \tan^{-1}\left(\frac{y}{x}\right):     J=(r,θ)(x,y)=xryryr2xr2=x2r3+y2r3=r2r3=1rJ' = \frac{\partial(r,\theta)}{\partial(x,y)} = \begin{vmatrix} \frac{x}{r} & \frac{y}{r} \\ -\frac{y}{r^2} & \frac{x}{r^2} \end{vmatrix} = \frac{x^2}{r^3} + \frac{y^2}{r^3} = \frac{r^2}{r^3} = \frac{1}{r}JJ=r×1r=1J J' = r \times \frac{1}{r} = 1
    • Example 3: Implicit Jacobian System:     If x2+y2+u2v2=0x^2 + y^2 + u^2 - v^2 = 0 and uv+xy=0uv + xy = 0, prove that:     (u,v)(x,y)=x2y2u2+v2\frac{\partial(u,v)}{\partial(x,y)} = \frac{x^2 - y^2}{u^2 + v^2}
    • Solution:       Let f1=x2+y2+u2v2=0f_1 = x^2 + y^2 + u^2 - v^2 = 0 and f2=uv+xy=0f_2 = uv + xy = 0.       (f1,f2)(x,y)=2x2yyx=2x22y2=2(x2y2)\frac{\partial(f_1, f_2)}{\partial(x,y)} = \begin{vmatrix} 2x & 2y \\ y & x \end{vmatrix} = 2x^2 - 2y^2 = 2(x^2 - y^2)(f1,f2)(u,v)=2u2vvu=2u2+2v2=2(u2+v2)\frac{\partial(f_1, f_2)}{\partial(u,v)} = \begin{vmatrix} 2u & -2v \\ v & u \end{vmatrix} = 2u^2 + 2v^2 = 2(u^2 + v^2)(u,v)(x,y)=(1)22(x2y2)2(u2+v2)=x2y2u2+v2\frac{\partial(u,v)}{\partial(x,y)} = (-1)^2 \frac{2(x^2 - y^2)}{2(u^2 + v^2)} = \frac{x^2 - y^2}{u^2 + v^2}
    • Example 4: Three-Variable Symmetric Determinant:     If u=xyzu = xyz, v=x2+y2+z2v = x^2 + y^2 + z^2, w=x+y+zw = x + y + z, find (u,v,w)(x,y,z)\frac{\partial(u,v,w)}{\partial(x,y,z)}.
    • Solution:       (u,v,w)(x,y,z)=yzzxxy2x2y2z111=2yzzxxyxyz111\frac{\partial(u,v,w)}{\partial(x,y,z)} = \begin{vmatrix} yz & zx & xy \\ 2x & 2y & 2z \\ 1 & 1 & 1 \end{vmatrix} = 2 \begin{vmatrix} yz & zx & xy \\ x & y & z \\ 1 & 1 & 1 \end{vmatrix}       Applying column operations C2C2C1C_2 \rightarrow C_2 - C_1 and C3C3C1C_3 \rightarrow C_3 - C_1:       =2(xy)(yz)(zx)= 2 (x-y)(y-z)(z-x)

Maxima and Minima for Functions of Two Variables

  • Definitions:

    • Local Maximum: A function f(x,y)f(x,y) has a maximum value at x=a,y=bx = a, y = b if f(a,b)>f(a+h,b+k)f(a,b) > f(a+h, b+k) for small positive or negative values of hh and kk
    • Local Minimum: A function f(x,y)f(x,y) has a minimum value at x=a,y=bx = a, y = b if f(a,b)<f(a+h,b+k)f(a,b) < f(a+h, b+k) for small positive or negative values of hh and kk
    • Critical / Stationary Points: Points (x1,y1),(x2,y2)(x_1, y_1), (x_2, y_2) where fx=0f_x = 0 and fy=0f_y = 0 simultaneously.
  • Working Method for Maxima and Minima:

    • Step I: Find fxf_x and fyf_y
    • Step II: Solve fx=0f_x = 0 and fy=0f_y = 0 simultaneously to find critical points (x1,y1),(x2,y2),(x_1, y_1), (x_2, y_2), \dots
    • Step III: Compute second partial derivatives:     A=fxx=2fx2,B=fxy=2fxy,C=fyy=2fy2A = f_{xx} = \frac{\partial^2 f}{\partial x^2}, \quad B = f_{xy} = \frac{\partial^2 f}{\partial x \partial y}, \quad C = f_{yy} = \frac{\partial^2 f}{\partial y^2}
    • Step IV: Evaluate ACB2AC - B^2 at each critical point:
    • (i) If ACB2>0AC - B^2 > 0 and A<0A < 0, then f(x,y)f(x,y) has a maximum value at that point, given by f(x1,y1)f(x_1, y_1).
    • (ii) If ACB2>0AC - B^2 > 0 and A>0A > 0, then f(x,y)f(x,y) has a minimum value at that point, given by f(x1,y1)f(x_1, y_1).
    • (iii) If ACB2<0AC - B^2 < 0, then f(x,y)f(x,y) has neither maximum nor minimum at that point. This point is called a saddle point.
    • (iv) If ACB2=0AC - B^2 = 0, the case is doubtful and requires further investigation of sign changes of f(a,b)f(a+h,b+k)f(a,b) - f(a+h, b+k).
  • Example 1: Exponential Function Extrema:

    • Show that f(x,y)=x3e(x2+y2)f(x,y) = x^3 e^{-(x^2+y^2)} has a maximum at (3/2,0)(\sqrt{3/2}, 0), a minimum at (3/2,0)(-\sqrt{3/2}, 0), and a stationary point at the origin whose nature is doubtful.
    • Solution:     fx=(3x22x4)e(x2+y2)=0f_x = (3x^2 - 2x^4) e^{-(x^2+y^2)} = 0fy=2x3ye(x2+y2)=0f_y = -2x^3 y e^{-(x^2+y^2)} = 0     Stationary points are (0,0)(0,0), (3/2,0)(\sqrt{3/2}, 0), and (3/2,0)(-\sqrt{3/2}, 0).
    • At (0,0)(0,0): A=0,B=0,C=0    ACB2=0A = 0, B = 0, C = 0 \implies AC - B^2 = 0 (Doubtful case).       Evaluating f(0,0)f(h,k)=h3e(h2+k2)f(0,0) - f(h,k) = -h^3 e^{-(h^2+k^2)} shows sign changes with sign of hh, confirming (0,0)(0,0) is not an extremum.
    • At (3/2,0)(\sqrt{3/2}, 0): A=63/2e3/2<0A = -6 \sqrt{3/2} e^{-3/2} < 0 and ACB2>0    AC - B^2 > 0 \implies Maximum.
    • At (3/2,0)(-\sqrt{3/2}, 0): A=63/2e3/2>0A = 6 \sqrt{3/2} e^{-3/2} > 0 and ACB2>0    AC - B^2 > 0 \implies Minimum.
  • Example 2: Cubic Polynomial Extrema:

    • Find all maxima and minima of f(x,y)=x3+y363(x+y)+12xyf(x,y) = x^3 + y^3 - 63(x+y) + 12xy
    • Solution Steps:
    • Step I & II: Set first derivatives to zero:       fx=3x263+12y=0    x2+4y21=0f_x = 3x^2 - 63 + 12y = 0 \implies x^2 + 4y - 21 = 0fy=3y263+12x=0    y2+4x21=0f_y = 3y^2 - 63 + 12x = 0 \implies y^2 + 4x - 21 = 0       Subtracting equations gives (xy)(x+y4)=0(x-y)(x+y-4) = 0       Solving system yields stationary points: (3,3)(3,3), (7,7)(-7,-7), (1,5)(1,5), and (5,1)(5,-1).
    • Step III & IV: Compute A=6xA = 6x, B=12B = 12, C=6yC = 6y:
      • At (3,3)(3,3): A=18>0A = 18 > 0, C=18C = 18, ACB2=324144=180>0    AC - B^2 = 324 - 144 = 180 > 0 \implies Minimum at (3,3)(3,3) with minimum value f(3,3)=216f(3,3) = -216
      • At (7,7)(-7,-7): A=42<0A = -42 < 0, C=42C = -42, ACB2=1764144=1620>0    AC - B^2 = 1764 - 144 = 1620 > 0 \implies Maximum at (7,7)(-7,-7) with maximum value f(7,7)=784f(-7,-7) = 784
      • At (1,5)(1,5) and (5,1)(5,-1): ACB2=36(5)(1)144=324<0    AC - B^2 = 36(5)(-1) - 144 = -324 < 0 \implies Saddle points.

Lagrange's Method of Undetermined Multipliers

  • Method Description:

    • To find maximum or minimum values of f(x,y,z)f(x,y,z) subject to constraint equation g(x,y,z)=0g(x,y,z) = 0:
    • Define auxiliary function:     F(x,y,z,λ)=f(x,y,z)+λg(x,y,z)F(x,y,z,\lambda) = f(x,y,z) + \lambda g(x,y,z)     where λ\lambda is Lagrange's undetermined multiplier.
    • Necessary conditions for extrema are:     Fx=0,Fy=0,Fz=0,Fλ=0\frac{\partial F}{\partial x} = 0, \quad \frac{\partial F}{\partial y} = 0, \quad \frac{\partial F}{\partial z} = 0, \quad \frac{\partial F}{\partial \lambda} = 0
  • Problem 1: Distance-Squared Minimization under Linear Constraint:

    • Find minimum value of x2+y2+z2x^2 + y^2 + z^2 subject to ax + by + cz = p$.\n - Solution:\n    Auxiliary function: F = (x^2 + y^2 + z^2) + \lambda (ax + by + cz - p)\n    Partial derivatives:\n    \frac{\partial F}{\partial x} = 2x + \lambda a = 0 \implies x = -\frac{\lambda a}{2}\n    \frac{\partial F}{\partial y} = 2y + \lambda b = 0 \implies y = -\frac{\lambda b}{2}\n    \frac{\partial F}{\partial z} = 2z + \lambda c = 0 \implies z = -\frac{\lambda c}{2}\n    Substituting into constraint ax + by + cz = p:\n    a\left(-\frac{\lambda a}{2}\right) + b\left(-\frac{\lambda b}{2}\right) + c\left(-\frac{\lambda c}{2}\right) = p \implies -\frac{\lambda}{2}(a^2 + b^2 + c^2) = p\n    \lambda = -\frac{2p}{a^2 + b^2 + c^2}\n    Coordinates of minimum point:\n    x = \frac{ap}{a^2+b^2+c^2}, \quad y = \frac{bp}{a^2+b^2+c^2}, \quad z = \frac{cp}{a^2+b^2+c^2}\n    Minimum value:\n    x^2 + y^2 + z^2 = \frac{a^2 p^2 + b^2 p^2 + c^2 p^2}{(a^2+b^2+c^2)^2} = \frac{p^2}{a^2+b^2+c^2}\n\n- **Problem 2: Open Rectangular Box Material Minimization**:\n - A rectangular box open at top is to have volume 32\text{ m}^3. Find dimensions for least material.\n - Solution:\n    Let x, y, z be length, breadth, height.\n    Surface Area S = xy + 2yz + 2zx\n    Constraint Volume V = xyz = 32\n    Auxiliary function: F = xy + 2yz + 2zx + \lambda(xyz - 32)\n    \frac{\partial F}{\partial x} = y + 2z + \lambda yz = 0\n    \frac{\partial F}{\partial y} = x + 2z + \lambda xz = 0\n    \frac{\partial F}{\partial z} = 2x + 2y + \lambda xy = 0\n    Solving system yields x = y = 2z.\n    Substituting into xyz = 32 \implies (2z)(2z)(z) = 32 \implies 4z^3 = 32 \implies z^3 = 8 \implies z = 2\text{ m}.\n    Optimal dimensions: x = 4\text{ m}, y = 4\text{ m}, z = 2\text{ m}.\n\n- **Problem 3: Multiple Constraints Optimization**:\n - Find extreme values of f(x,y,z) = 2x + 3y + zsubjecttosubject tox^2 + y^2 = 5andandx + z = 1\n - Solution:\n    F = 2x + 3y + z + \lambda_1 (x^2 + y^2 - 5) + \lambda_2 (x + z - 1)\n    \frac{\partial F}{\partial z} = 1 + \lambda_2 = 0 \implies \lambda_2 = -1\n    \frac{\partial F}{\partial y} = 3 + 2\lambda_1 y = 0 \implies y = -\frac{3}{2\lambda_1}\n    \frac{\partial F}{\partial x} = 2 + 2\lambda_1 x + \lambda_2 = 0 \implies 1 + 2\lambda_1 x = 0 \implies x = -\frac{1}{2\lambda_1}\n    Substituting xandandyintointox^2 + y^2 = 5:\n    \left(-\frac{1}{2\lambda_1}\right)^2 + \left(-\frac{3}{2\lambda_1}\right)^2 = 5 \implies \frac{10}{4\lambda_1^2} = 5 \implies \lambda_1 = \pm \frac{1}{\sqrt{2}}\n - For \lambda_1 = \frac{1}{\sqrt{2}}::x = -\frac{1}{\sqrt{2}},,y = -\frac{3}{\sqrt{2}},,z = 1 + \frac{1}{\sqrt{2}},givingminimum, giving minimumf = 1 - 5\sqrt{2}.\n - For \lambda_1 = -\frac{1}{\sqrt{2}}::x = \frac{1}{\sqrt{2}},,y = \frac{3}{\sqrt{2}},,z = 1 - \frac{1}{\sqrt{2}},givingmaximum, giving maximumf = 1 + 5\sqrt{2}.\n\n- **Problem 4: Inscribed Box in Ellipsoid**:\n - Find maximum volume of rectangular parallelopiped inscribed in ellipsoid \frac{x^2}{a^2} + \frac{y^2}{b^2} + \frac{z^2}{c^2} = 1.\n - Vertices at (\pm x, \pm y, \pm z),volume, volumeV = (2x)(2y)(2z) = 8xyz\n - Auxiliary function F = 8xyz + \lambda\left(\frac{x^2}{a^2} + \frac{y^2}{b^2} + \frac{z^2}{c^2} - 1\right)\n - Differentiating and setting to zero yields \frac{x^2}{a^2} = \frac{y^2}{b^2} = \frac{z^2}{c^2} = \frac{1}{3}\n - Dimensions: x = \frac{a}{\sqrt{3}},,y = \frac{b}{\sqrt{3}},,z = \frac{c}{\sqrt{3}}\n - Maximum Volume V = \frac{8abc}{3\sqrt{3}}\n\n# Differentiation Under the Integral Sign (Leibniz Rule)\n\n![Differentiation under the integral sign and Leibniz rule example](https://assets.knowt.com/pdf-flow-prod/f357703a-3a80-4c71-81d8-7edee760f945-figures/1.jpg)\n\n- **General Leibniz Rule (Variable Limits of Integration)**:\n - If f(x,t)iscontinuousandcontinuouslydifferentiable(partialderivativesexistandarecontinuous)andintegrationlimitsis continuous and continuously differentiable (partial derivatives exist and are continuous) and integration limitsa(x)andandb(x)arecontinuouslydifferentiablefunctionsofare continuously differentiable functions ofx, then:\n    \frac{d}{dx} \int_{a(x)}^{b(x)} f(x,t) \, dt = f(x, b(x)) b'(x) - f(x, a(x)) a'(x) + \int_{a(x)}^{b(x)} \frac{\partial}{\partial x} f(x,t) \, dt\n\n- **Simpler Leibniz Integral Rule (Constant Limits)**:\n - In case where a(x) = aandandb(x) = b are constant functions, the formula reduces to:\n    \frac{d}{dx} \int_{a}^{b} f(x,t) \, dt = \int_{a}^{b} \frac{\partial}{\partial x} f(x,t) \, dt\n\n- **Example 1: Variable Limits Differentiation**:\n - Find derivative with respect to x of:\n    I(x) = \int_{x}^{x^2} \frac{\sin(xt)}{t} \, dt\n - Solution:\n    Here f(x,t) = \frac{\sin(xt)}{t},,a(x) = x \implies a'(x) = 1,,b(x) = x^2 \implies b'(x) = 2x\n    Applying Leibniz Rule:\n    \frac{dI}{dx} = f(x, x^2) (2x) - f(x, x) (1) + \int_{x}^{x^2} \frac{\partial}{\partial x} \left(\frac{\sin(xt)}{t}\right) dt\n    = \frac{\sin(x^3)}{x^2} (2x) - \frac{\sin(x^2)}{x} (1) + \int_{x}^{x^2} \cos(xt) \, dt\n    = \frac{2 \sin(x^3)}{x} - \frac{\sin(x^2)}{x} + \left[ \frac{\sin(xt)}{x} \right]{x}^{x^2}\n    = \frac{2 \sin(x^3)}{x} - \frac{\sin(x^2)}{x} + \frac{\sin(x^3)}{x} - \frac{\sin(x^2)}{x}\n    \frac{dI}{dx} = \frac{3 \sin(x^3) - 2 \sin(x^2)}{x}\n\n- **Example 2: Definite Integral Evaluation via Parameter Differentiation**:\n - Evaluate I(\alpha) = \int{0}^{1} \frac{x^\alpha - 1}{\log(x)} \, dx\n - Solution Steps:\n - At \alpha = 0,,I(0) = \int_{0}^{1} \frac{x^0 - 1}{\log(x)} \, dx = 0\n - Differentiating under integral sign with respect to parameter \alpha:\n      \frac{dI}{d\alpha} = \int_{0}^{1} \frac{\partial}{\partial \alpha} \left(\frac{x^\alpha - 1}{\log(x)}\right) dx = \int_{0}^{1} \frac{x^\alpha \log(x)}{\log(x)} \, dx = \int_{0}^{1} x^\alpha \, dx\n - Integrating with respect to x:\n      \frac{dI}{d\alpha} = \left[ \frac{x^{\alpha+1}}{\alpha + 1} \right]{0}^{1} = \frac{1}{\alpha + 1}\n - Integrating with respect to \alpha:\n      I(\alpha) = \log(\alpha + 1) + C\n - Since I(0) = 0 \implies C = 0, we deduce:\n      I(\alpha) = \log(\alpha + 1)\n\n- **Example 3: Logarithmic Rational Integral Evaluation**:\n - Evaluate F(\alpha) = \int{0}^{1} \frac{\log(1 + \alpha x)}{1 + x^2} \, dx\n - Differentiating under integral sign:\n    \frac{dF}{d\alpha} = \int_{0}^{1} \frac{x}{(1 + \alpha x)(1 + x^2)} \, dx\n - Using partial fractions and integrating yields:\n    F(\alpha) = \frac{1}{2} \log(1 + \alpha^2) \tan^{-1}(\alpha)\n - For \alpha = 1:\n    \int_{0}^{1} \frac{\log(1 + x)}{1 + x^2} \, dx = \frac{\pi}{8} \log(2)$$