The math optional, made finite. Daily Practice

Quasilinear first-order PDEs (Lagrange’s method)

At a Glance

Why This Chapter Matters

Quasilinear first-order PDEs appear in 7 of the last 12 years, split evenly between 10- and 15-mark slots and rated the hardest PDE atom in Paper 2. The difficulty is almost entirely in finding the right multipliers for Lagrange’s auxiliary equations — once the two first integrals are found, the general solution is immediate. Three question types repeat: pure Lagrange (find two first integrals, write the general solution), integral surface through a curve (additionally eliminate a parameter to pin down a specific surface), and orthogonal trajectories (form the PDE from the family’s gradient and apply Lagrange). The multiplier patterns that recur — (x,y,z)(x,y,z), (y,x,0)(y,x,0), (3x2,3y2,0)(3x^2,3y^2,0), and “add all three ratios” — are worth memorising.

Minimum Theory

Lagrange’s quasilinear PDE. The PDE P(x,y,z)p+Q(x,y,z)q=R(x,y,z)P(x,y,z)p+Q(x,y,z)q=R(x,y,z) (where p=zxp=z_x, q=zyq=z_y) has general solution F(u1,u2)=0F(u_1,u_2)=0 where u1=c1u_1=c_1 and u2=c2u_2=c_2 are two independent first integrals of the Lagrange auxiliary (subsidiary) equations: dxP=dyQ=dzR.\frac{dx}{P}=\frac{dy}{Q}=\frac{dz}{R}.

Finding first integrals via multipliers. For a common ratio k=dx/P=dy/Q=dz/Rk=dx/P=dy/Q=dz/R, any multipliers (,m,n)(\ell,m,n) give k=(dx+mdy+ndz)/(P+mQ+nR)k=(\ell\,dx+m\,dy+n\,dz)/(\ell P+mQ+nR). If P+mQ+nR=0\ell P+mQ+nR=0 then dx+mdy+ndz=0\ell\,dx+m\,dy+n\,dz=0 along characteristics, giving an exact first integral when the numerator is integrable. Good multiplier choices: (x,y,z)(x,y,z) (looks for x2+y2+z2x^2+y^2+z^2 structure); (y,x,0)(y,x,0) (gives d(xy)d(xy)); (1,1,1)(1,1,1) (gives x+y+zx+y+z); (x,y,ϕ)(x,y,\phi) (gives x2+y2+ϕ2x^2+y^2+\phi^2).

Integral surface through a curve. Given parametric curve x=x0(t),  y=y0(t),  z=z0(t)x=x_0(t),\;y=y_0(t),\;z=z_0(t): evaluate c1=u1(x0,y0,z0)c_1=u_1(x_0,y_0,z_0) and c2=u2(x0,y0,z0)c_2=u_2(x_0,y_0,z_0) as functions of the parameter tt; eliminate tt to get the relation F(c1,c2)=0F(c_1,c_2)=0 that the integral surface satisfies.

Orthogonal trajectory surfaces. For a one-parameter family ϕ(x,y,z)=C\phi(x,y,z)=C, the normal direction at each point is ϕ\nabla\phi; orthogonal surfaces satisfy the PDE ϕxp+ϕyq=ϕz\phi_x p+\phi_y q=\phi_z (after eliminating CC). Solve by Lagrange’s method.

Question Archetypes

Three patterns cover every quasilinear PDE question in the corpus.

ArchetypeYou are seeing this when…
integral-surfaceLagrange PDE + a curve; find the specific integral surface through that curve
lagrange-linearLagrange PDE; find the general solution (two first integrals)
orthogonal-surfacesa one-parameter family of surfaces; find the family of orthogonal trajectory surfaces

integral-surface (4 question(s); 2018, 2019, 2020, 2024)

Recognition Cues

Solution Template

  1. Write the auxiliary equations dx/P=dy/Q=dz/Rdx/P=dy/Q=dz/R.
  2. Find two independent first integrals u1=c1u_1=c_1 and u2=c2u_2=c_2 using multipliers or direct pairing.
  3. Write the general solution F(u1,u2)=0F(u_1,u_2)=0 (or u2=G(u1)u_2=G(u_1) equivalently).
  4. On the given curve, compute c1(t)=u1(x0(t),y0(t),z0(t))c_1(t)=u_1(x_0(t),y_0(t),z_0(t)) and c2(t)=u2(x0(t),y0(t),z0(t))c_2(t)=u_2(x_0(t),y_0(t),z_0(t)).
  5. Eliminate the parameter tt between c1(t)c_1(t) and c2(t)c_2(t) to find the relation F(c1,c2)=0F(c_1,c_2)=0.
  6. Substitute back u1,u2u_1,u_2 to write the surface equation explicitly.

Worked Example(s)

2018 Paper 2, 2018-P2-Q6a (15 marks)

(y3x2x4)p+(2y4x3y)q=9z(x3y3)(y^3x-2x^4)p+(2y^4-x^3y)q=9z(x^3-y^3); integral surface through x=t,y=t2,z=1x=t,y=t^2,z=1.

Auxiliary equations: dxx(y32x3)=dyy(2y3x3)=dz9z(x3y3)\dfrac{dx}{x(y^3-2x^3)}=\dfrac{dy}{y(2y^3-x^3)}=\dfrac{dz}{9z(x^3-y^3)}.

First integral (multipliers (y,x,0)(y,x,0)): yP+xQ=3xy(x3y3)yP+xQ=-3xy(x^3-y^3); numerator =d(xy)=d(xy). Pairing with dz/Rdz/R and cancelling (x3y3)(x^3-y^3): d(xy)/(3xy)=dz/(9z)d(xy)/(-3xy)=dz/(9z), integrating gives u1=z(xy)3=c1u_1=z(xy)^3=c_1.

Second integral (multipliers (3x2,3y2,0)(3x^2,3y^2,0)): denominator =6(x3y3)(x3+y3)=-6(x^3-y^3)(x^3+y^3); numerator =d(x3+y3)=d(x^3+y^3). Same cancel gives u2=z2(x3+y3)3=c2u_2=z^2(x^3+y^3)^3=c_2.

General solution: Φ(z(xy)3,  z2(x3+y3)3)=0\Phi(z(xy)^3,\;z^2(x^3+y^3)^3)=0.

Integral surface. On the curve (x,y,z)=(t,t2,1)(x,y,z)=(t,t^2,1): c1=t9c_1=t^9, c2=(t3+t6)3=t9(1+t3)3c_2=(t^3+t^6)^3=t^9(1+t^3)^3. Eliminate tt: c11/3=t3c_1^{1/3}=t^3, then c21/3=t3(1+t3)=c11/3(1+c11/3)c_2^{1/3}=t^3(1+t^3)=c_1^{1/3}(1+c_1^{1/3}). After cubing:

  z(x3+y3x2y2)3=(xy)3.  \boxed{\;z(x^3+y^3-x^2y^2)^3=(xy)^3.\;}


2019 Paper 2, 2019-P2-Q6a (15 marks)

xux+(uxy)uy=x+2yxu_x+(u-x-y)u_y=x+2y; u=1+yu=1+y on x=1x=1.

Auxiliary: dxx=dyuxy=dux+2y\dfrac{dx}{x}=\dfrac{dy}{u-x-y}=\dfrac{du}{x+2y}.

First integral: add dy+dudy+du: numerator d(u+y)d(u+y), denominator (u+y)(u+y). Pairing with dx/xdx/x: d(u+y)/(u+y)=dx/xd(u+y)/(u+y)=dx/x, giving C1=(u+y)/xC_1=(u+y)/x.

Second integral: along characteristic with x=esx=e^s, yy satisfies the linear ODE y+2y=(C11)esy'+2y=(C_1-1)e^s; solution y=A/x2+(C11)x/3y=A/x^2+(C_1-1)x/3. Then C2=x2(x+2yu)C_2=x^2(x+2y-u).

Applying u=1+yu=1+y on x=1x=1: C1=1+2yC_1=1+2y, C2=yC_2=y. Eliminate yy: C1=1+2C2C_1=1+2C_2. Solving for uu:

  u(x,y)=2x4+4x3y+xy1+2x3.  \boxed{\;u(x,y)=\dfrac{2x^4+4x^3y+x-y}{1+2x^3}.\;}


2020 Paper 2, 2020-P2-Q6a (15 marks)

(xy)y2zx+(yx)x2zy=(x2+y2)z(x-y)y^2z_x+(y-x)x^2z_y=(x^2+y^2)z; integral surface through xz=a3,y=0xz=a^3,y=0.

First integral: cancel the shared factor (xy)(x-y) from P,QP,Q: dx/y2=dy/(x2)dx/y^2=dy/(-x^2), giving x2dx+y2dy=0x^2\,dx+y^2\,dy=0, so u1=x3+y3=c1u_1=x^3+y^3=c_1.

Second integral: test v=z/(xy)v=z/(x-y); verify Pvx+Qvy+Rvz=0P v_x+Q v_y+R v_z=0 ✓, so u2=z/(xy)=c2u_2=z/(x-y)=c_2.

Applying curve y=0y=0, z=a3/xz=a^3/x: c1=x3c_1=x^3, c2=z/x=a3/x2c_2=z/x=a^3/x^2. Eliminate xx: c1=x3x=c11/3c_1=x^3\Rightarrow x=c_1^{1/3}, so c2c12/3=a3c_2 c_1^{2/3}=a^3, i.e. c23c12=a9c_2^3c_1^2=a^9.

  (x3+y3)2z3=a9(xy)3.  \boxed{\;(x^3+y^3)^2 z^3=a^9(x-y)^3.\;}


2024 Paper 2, 2024-P2-Q7a (15 marks)

(yϕ)ϕx+(ϕx)ϕy=xy(y-\phi)\phi_x+(\phi-x)\phi_y=x-y; through ϕ=0,xy=1\phi=0,xy=1 and circle x+y+ϕ=0,x2+y2+ϕ2=a2x+y+\phi=0,x^2+y^2+\phi^2=a^2.

First integral: add all three ratios — denominator =(yϕ)+(ϕx)+(xy)=0=(y-\phi)+(\phi-x)+(x-y)=0, so dx+dy+dϕ=0dx+dy+d\phi=0, giving C1=x+y+ϕC_1=x+y+\phi.

Second integral: multiply numerators by x,y,ϕx,y,\phi: denominator x(yϕ)+y(ϕx)+ϕ(xy)=0x(y-\phi)+y(\phi-x)+\phi(x-y)=0, so xdx+ydy+ϕdϕ=0x\,dx+y\,dy+\phi\,d\phi=0, giving C2=x2+y2+ϕ2C_2=x^2+y^2+\phi^2.

General solution: x2+y2+ϕ2=G(x+y+ϕ)x^2+y^2+\phi^2=G(x+y+\phi).

Through ϕ=0,xy=1\phi=0,xy=1: C1=x+yC_1=x+y, C2=x2+y2=(x+y)22xy=C122C_2=x^2+y^2=(x+y)^2-2xy=C_1^2-2. So G(t)=t22G(t)=t^2-2.

Surface: x2+y2+ϕ2=(x+y+ϕ)22x^2+y^2+\phi^2=(x+y+\phi)^2-2. Expand and simplify:

  xy+yϕ+xϕ=1.  \boxed{\;xy+y\phi+x\phi=1.\;}

(The second prescribed condition — the circle — gives a2=2a^2=-2, impossible for real aa; flag this inconsistency.)

Common Traps


lagrange-linear (3 question(s); 2015, 2015, 2016)

Recognition Cues

Solution Template

  1. Write auxiliary equations dx/P=dy/Q=dz/Rdx/P=dy/Q=dz/R.
  2. Find the two most accessible first integrals by: (a) pairing two ratios directly (dy/y=dz/zdy/y=dz/z etc.); (b) using multipliers to build exact differentials.
  3. State the general integral F(u1,u2)=0F(u_1,u_2)=0 (or u2=f(u1)u_2=f(u_1)).

Worked Example(s)

2015 Paper 2, 2015-P2-Q5a (10 marks)

(y2+z2x2)p2xyq+2xz=0(y^2+z^2-x^2)p-2xyq+2xz=0.

Standard form: Pp+Qq=RPp+Qq=R with R=2xzR=-2xz. Auxiliary: dxy2+z2x2=dy2xy=dz2xz\dfrac{dx}{y^2+z^2-x^2}=\dfrac{dy}{-2xy}=\dfrac{dz}{-2xz}.

First integral: pair dy/(2xy)=dz/(2xz)dy/(-2xy)=dz/(-2xz): dy/y=dz/zdy/y=dz/z, giving u1=y/z=c1u_1=y/z=c_1.

Second integral: multipliers (x,y,z)(x,y,z): xP+yQ+zR=xz2+xy2x32xy22xz2=x(x2+y2+z2)xP+yQ+zR=xz^2+xy^2-x^3-2xy^2-2xz^2=-x(x^2+y^2+z^2); numerator xdx+ydy+zdz=12d(x2+y2+z2)x\,dx+y\,dy+z\,dz=\tfrac{1}{2}d(x^2+y^2+z^2). Pairing with dz/(2xz)dz/(-2xz) and cancelling x-x: gives u2=(x2+y2+z2)/z=c2u_2=(x^2+y^2+z^2)/z=c_2.

  Φ ⁣(yz,  x2+y2+z2z)=0.  \boxed{\;\Phi\!\left(\frac{y}{z},\;\frac{x^2+y^2+z^2}{z}\right)=0.\;}


2015 Paper 2, 2015-P2-Q6a (15 marks)

pcos(x+y)+qsin(x+y)=zp\cos(x+y)+q\sin(x+y)=z.

Let u=x+yu=x+y, v=xyv=x-y. Auxiliary: du/(cosu+sinu)=dv/(cosusinu)=dz/zdu/(\cos u+\sin u)=dv/(\cos u-\sin u)=dz/z.

First integral (from dudu and dzdz): integrate dz/z=du/(cosu+sinu)dz/z=du/(\cos u+\sin u). Use cosu+sinu=2sin(u+π/4)\cos u+\sin u=\sqrt{2}\sin(u+\pi/4); cscvdv=lntan(v/2)\int\csc v\,dv=\ln|\tan(v/2)|: C1=[tan((x+y)/2+π/8)]1/2z.C_1=\frac{[\tan((x+y)/2+\pi/8)]^{1/\sqrt{2}}}{z}.

Second integral (from dudu and dvdv): dv/du=(cosusinu)/(cosu+sinu)=d(lncosu+sinu)/dudv/du=(\cos u-\sin u)/(\cos u+\sin u)=d(\ln|\cos u+\sin u|)/du, so: C2=xylnsin(x+y)+cos(x+y).C_2=x-y-\ln|\sin(x+y)+\cos(x+y)|.

  Φ ⁣(C1,  xylnsin(x+y)+cos(x+y))=0.  \boxed{\;\Phi\!\Big(C_1,\;x-y-\ln|\sin(x+y)+\cos(x+y)|\Big)=0.\;}


2016 Paper 2, 2016-P2-Q5e (10 marks)

(y+zx)p(x+yz)q=x2y2(y+zx)p-(x+yz)q=x^2-y^2.

First integral (multipliers (x,y,z)(x,y,z)): xP+yQ+zR=xz2yz2+zx2zy2+z(x2y2)=xP+yQ+zR=xz^2-yz^2+zx^2-zy^2+z(x^2-y^2)= simplifies to z(x2y2)=zRz(x^2-y^2)=zR. Numerator xdx+ydy+zdz=12d(x2+y2+z2)x\,dx+y\,dy+z\,dz=\tfrac12 d(x^2+y^2+z^2); relation d(x2+y2z2)=0d(x^2+y^2-z^2)=0 gives u1=x2+y2z2=c1u_1=x^2+y^2-z^2=c_1.

Wait — re-examine: xP+yQ+zR=x(y+zx)y(x+yz)+z(x2y2)=x2zy2z=z(x2y2)=zRxP+yQ+zR=x(y+zx)-y(x+yz)+z(x^2-y^2)=x^2z-y^2z=z(x^2-y^2)=zR. So xdx+ydy+zdzz(x2y2)=dzx2y2\frac{x\,dx+y\,dy+z\,dz}{z(x^2-y^2)}=\frac{dz}{x^2-y^2}, giving xdx+ydy=zdzzdz=0x\,dx+y\,dy=z\,dz-z\,dz=0. Actually more directly: xdx+ydy+zdz=zdzx\,dx+y\,dy+z\,dz=z\,dz, i.e. xdx+ydyzdz=0x\,dx+y\,dy-z\,dz=0, integrating to u1=x2+y2z2u_1=x^2+y^2-z^2.

Second integral (multipliers (y,x,1)(y,x,1)): yP+xQ+R=(y2x2)=RyP+xQ+R=(y^2-x^2)=-R. So ydx+xdy=dzy\,dx+x\,dy=-dz, giving d(xy+z)=0d(xy+z)=0, u2=xy+z=c2u_2=xy+z=c_2.

  F(x2+y2z2,  xy+z)=0.  \boxed{\;F(x^2+y^2-z^2,\;xy+z)=0.\;}

Common Traps


orthogonal-surfaces (2 question(s); 2013, 2016)

Recognition Cues

Solution Template

  1. Eliminate the parameter CC from F=CF=C by differentiating partially or substituting CC as a function of x,y,zx,y,z.
  2. Compute the normal direction (ϕx,ϕy,ϕz)(\phi_x,\phi_y,\phi_z) of the family (after eliminating CC).
  3. The orthogonal surface satisfies ϕxp+ϕyq=ϕz\phi_x p+\phi_y q=\phi_z (or ϕxdx=ϕydy=ϕzdz\phi_x\,dx=\phi_y\,dy=\phi_z\,dz, equivalently the surface contains the family’s normal).
  4. Write Lagrange’s auxiliary: dx/ϕx=dy/ϕy=dz/ϕzdx/\phi_x=dy/\phi_y=dz/\phi_z.
  5. Find two first integrals and write the general surface F(u1,u2)=0F(u_1,u_2)=0.
  6. If a specific surface is required, apply an initial curve condition.

Worked Example(s)

2013 Paper 2, 2013-P2-Q6b (15 marks)

Orthogonal to z(x+y)=C(3z+1)z(x+y)=C(3z+1); through circle x2+y2=1,z=1x^2+y^2=1,z=1.

Family: ϕ=z(x+y)/(3z+1)=C\phi=z(x+y)/(3z+1)=C. Normal: ϕx=ϕy=z/(3z+1)\phi_x=\phi_y=z/(3z+1), ϕz=(x+y)/(3z+1)2\phi_z=(x+y)/(3z+1)^2.

Auxiliary (multiply through by (3z+1)(3z+1)): dxz=dyz=dz(3z+1)x+y\dfrac{dx}{z}=\dfrac{dy}{z}=\dfrac{dz(3z+1)}{x+y}.

First integral: dx=dydx=dyu1=xy=c1u_1=x-y=c_1.

Second integral: with y=xc1y=x-c_1, x+y=2xc1x+y=2x-c_1; cross-multiply: (2xc1)dx=(3z2+z)dz(2x-c_1)\,dx=(3z^2+z)\,dz. Integrate: x2c1x=z3+z2/2+c2x^2-c_1x=z^3+z^2/2+c_2, i.e. u2=xyz3z2/2=c2u_2=xy-z^3-z^2/2=c_2.

Through x2+y2=1,z=1x^2+y^2=1,z=1: c1=xyc_1=x-y, c2=xy3/2c_2=xy-3/2. Using (xy)2=12xy(x-y)^2=1-2xy: xy=(1c12)/2xy=(1-c_1^2)/2, so c2=(1c12)/23/2=1c12/2c_2=(1-c_1^2)/2-3/2=-1-c_1^2/2, i.e. c12+2c2+2=0c_1^2+2c_2+2=0.

Substituting back:   x2+y2=2z3+z22.  \boxed{\;x^2+y^2=2z^3+z^2-2.\;}


2016 Paper 2, 2016-P2-Q5a (10 marks)

Surfaces orthogonal to x2+y2+z2=czx^2+y^2+z^2=cz.

Eliminate c=(x2+y2+z2)/zc=(x^2+y^2+z^2)/z: normal =(2x,2y,(z2x2y2)/z)=(2x,2y,(z^2-x^2-y^2)/z), equivalently (2xz,2yz,z2x2y2)(2xz,2yz,z^2-x^2-y^2).

Auxiliary: dx2xz=dy2yz=dzz2x2y2\dfrac{dx}{2xz}=\dfrac{dy}{2yz}=\dfrac{dz}{z^2-x^2-y^2}.

First integral: dx/x=dy/ydx/x=dy/yu1=y/x=c1u_1=y/x=c_1.

Second integral (multipliers (x,y,z)(x,y,z)): numerator xdx+ydy+zdz=12d(r2)x\,dx+y\,dy+z\,dz=\tfrac12 d(r^2) where r2=x2+y2+z2r^2=x^2+y^2+z^2; denominator z(x2+y2+z2)=zr2z(x^2+y^2+z^2)=zr^2. Pairing with dy/(2yz)dy/(2yz): d(r2)2r2=dy2y\frac{d(r^2)}{2r^2}=\frac{dy}{2y}, integrating to u2=(x2+y2+z2)/y=c2u_2=(x^2+y^2+z^2)/y=c_2.

  Φ ⁣(yx,  x2+y2+z2y)=0.  \boxed{\;\Phi\!\left(\frac{y}{x},\;\frac{x^2+y^2+z^2}{y}\right)=0.\;}

Common Traps


Marks-Aware Writing

10-mark questions (2015-Q5a, 2016-Q5a, 2016-Q5e): Two to three lines of multiplier setup, one boxed general solution. Show the multiplier calculation (denominator =0=0, numerator =d()=d(\ldots)) explicitly — that is the method step worth marks.

15-mark questions (2013, 2015-Q6a, 2018, 2019, 2020, 2024): For integral-surface questions: four phases — auxiliary equations, two first integrals, general solution, surface through the curve. Each phase is ~3 marks; don’t collapse them. For orthogonal-surface questions: show the gradient computation and the resulting Lagrange PDE before solving.

Practice Set

YearPaper/QMarksOne-line hint
2022P2-Q5a10Lagrange: find two first integrals from the auxiliary equations; straightforward multiplier approach
2020P2-Q5d10Lagrange compulsory-part question; apply “add ratios” or direct pairing
2020P2-Q7a15Integral surface; two first integrals then eliminate curve parameter
2017P2-Q6a15Integral surface through given initial curve
2025P2-Q7a15Integral surface; follow the same four-phase template

We've mapped all 13 years of this exam. Get new chapters, tools, and solutions as we release them — free.