MATH529 L05
Boundary Value Problems

Rectangular coordinates

Separation of variables

Try u(x,y)=X(x)Y(y) in (Subscript[u, x x])=4 (Subscript[u, y])⇔(Subscript[u, x x])- 4 (Subscript[u, y])=0, homogeneous PDE

PDE = (∂_(x,x))u[x,y]==4 (∂_y)u[x,y]; vPDE = (∂_(x,x))v[x,y]==4 (∂_y)v[x,y];
u[x_,y_]=X[x] Y[y];
sPDE=Assuming[u[x,y]≠0,DivideSides[PDE,u[x,y]]]

xDE = sPDE[[1]] == -λ; yDE = sPDE[[2]] == -λ; {xDE,yDE}

ODE solutions depend on sign of λ, consider each possible case

Case   λ=0

xDE0 = xDE /. λ->0; yDE0 = yDE /. λ->0;
{xDE0,yDE0}

xsol0=DSolve[xDE0,X[x],x][[1,1]]; ysol0=DSolve[yDE0,Y[y],y][[1,1]] /. (C[1])->(C[3]) ;
sol0={xsol0,ysol0}

u0[x_,y_] = X[x] Y[y] /. sol0

Verify solution

v[x_,y_] = u0[x,y];
vPDE

Clear[v]; vPDE

Case   λ<0

xDEn = xDE /. λ->-(α^2); yDEn = yDE /. λ->-(α^2);
{xDEn,yDEn}

xsoln=DSolve[xDEn,X[x],x][[1,1]]/. {(C[1])->(C[4]),(C[2])->(C[5])}; ysoln=DSolve[yDEn,Y[y],y][[1,1]] /.(C[1])->(C[6]) ;
soln={xsoln,ysoln}

un[x_,y_] = X[x] Y[y] /. soln

Verify solution

v[x_,y_] = un[x,y];
Simplify[vPDE ]

Clear[v]; vPDE

Case   λ>0

xDEp = xDE /. λ->(α^2); yDEp = yDE /. λ->(α^2);
{xDEp,yDEp}

xsolp=DSolve[xDEp,X[x],x][[1,1]]/. {(C[1])->(C[7]),(C[2])->(C[8])}; ysolp=DSolve[yDEp,Y[y],y][[1,1]] /.(C[1])->(C[9]) ;
solp={xsolp,ysolp}

up[x_,y_] = X[x] Y[y] /. solp

Verify solution

v[x_,y_] = up[x,y];
Simplify[vPDE ]

Clear[v]; vPDE

Superposition of solutions

v[x_,y_]= a0 u0[x,y] + ap up[x,y] + an un[x,y]

Simplify[vPDE ]