Application Center - Maplesoft

Contact Maplesoft Request Quote

Submit your work

Application Center Applications Boundary value problems for 3-D elliptic equations (updated to Maple 7) Preview

App Preview:

Boundary value problems for 3-D elliptic equations (updated to Maple 7)

You can switch back to the summary page by clicking here.

Learn about Maple

Download Application

elliptic3d.mws

1 Example
2 Example
3 Example
4 Example
5 Example
6 Example
7 Example
8 Example

Boundary value problems for three-dimensional elliptic equations

by Aleksas Domarkas

Vilnius University, Faculty of Mathematics and Informatics,

Naugarduko 24, Vilnius, Lithuania

aleksas@ieva.mif.vu.lt

NOTE: In this session we find solutions of boundary value problems for three-dimensional

Laplace's , Poisson or Helmholtz equations in regions r<R, r>R and R1<r<R2.

1 Example

Problem

Solve boundary-value problem

, r1< r < r2,

u| r=r 1 = u1, u| r=r 2 = u2

here r = sqrt(x^2+y^2+z^2), u = u(x,y,z) ; -- constants.

Solving method

We consider problem

, r1< r < r2,

u| r=r1 =0, u| r=r2 =0

and problem

, r1< r < r2,

w| r=r1 =u1, w| r=r2 =u2.

Then u+w is solution of our problem. Solutions of these problems are sought in the form

u = v(r)*cos(phi)*sin(theta)*cos(theta),

Solving problem

`>`	restart;

`>`	with(linalg,laplacian):

`>`	*l:=expand(laplacian(u(r,theta,phi),[r,theta,phi], coords=spherical))-xz;**

l := 2/r*diff(u(r,theta,phi),r)+diff(u(r,theta,phi),`$`(r,2))+1/r^2/sin(theta)*cos(theta)*diff(u(r,theta,phi),theta)+1/r^2*diff(u(r,theta,phi),`$`(theta,2))+1/r^2/sin(theta)^2*diff(u(r,theta,phi),`$`(p...

`>`	*tr:={x=rcos(phi)sin(theta),y=rsin(phi)sin(theta),z=rcos(theta)};**

`>`	l := subs(tr,l);

`>`	*subs(u(r,theta,phi)=v(r)cos(phi)sin(theta)cos(theta),l):**

`>`	*l1 := simplify(%/(cos(phi)sin(theta)cos(theta)));*

l1 := (diff(v(r),`$`(r,2))*r^2-6*v(r)+2*diff(v(r),r)*r-r^4)/r^2

`>`	spr := dsolve(l1,{v(r)});

spr := v(r) = r^2*_C2+1/r^3*_C1+1/14*r^4

`>`	solve({subs(r=r1,rhs(spr)),subs(r=r2,rhs(spr))},{_C1,_C2});

${_C2 = -1/14*(r1^6+r2*r1^5+r2^2*r1^4+r2^3*r1^3+r2^4*r1^2+r2^5*r1+r2^6)/(r1^4+r1^3*r2+r1^2*r2^2+r1*r2^3+r2^4), _C1 = 1/14*r2^5*r1^5*(r1+r2)/(r1^4+r1^3*r2+r1^2*r2^2+r1*r2^3+r2^4)}$
${_C2 = -1/14*(r1^6+r2*r1^5+r2^2*r1^4+r2^3*r1^3+r2^4*r1^2+r2^5*r1+r2^6)/(r1^4+r1^3*r2+r1^2*r2^2+r1*r2^3+r2^4), _C1 = 1/14*r2^5*r1^5*(r1+r2)/(r1^4+r1^3*r2+r1^2*r2^2+r1*r2^3+r2^4)}$

`>`	v(r) := subs(%,rhs(spr));

v(r) := -1/14*r^2*(r1^6+r2*r1^5+r2^2*r1^4+r2^3*r1^3+r2^4*r1^2+r2^5*r1+r2^6)/(r1^4+r1^3*r2+r1^2*r2^2+r1*r2^3+r2^4)+1/14/r^3*r2^5*r1^5*(r1+r2)/(r1^4+r1^3*r2+r1^2*r2^2+r1*r2^3+r2^4)+1/14*r^4

`>`	subs(u(r,theta,phi)=w(r),laplacian(u(r,theta,phi),[r,theta,phi], coords=spherical));

1/r^2/sin(theta)*(2*r*sin(theta)*diff(w(r),r)+r^2*sin(theta)*diff(w(r),`$`(r,2))+cos(theta)*diff(w(r),theta)+sin(theta)*diff(w(r),`$`(theta,2))+1/sin(theta)*diff(w(r),`$`(phi,2)))

`>`	dsolve(%,w(r));

`>`	solve({subs(r=r1,rhs(%))=u1,subs(r=r2,rhs(%))=u2},{_C1,_C2});

${_C1 = (-u2*r2+u1*r1)/(-r2+r1), _C2 = -r2*r1*(u1-u2)/(-r2+r1)}$

`>`	w(r) := subs(%,rhs(%%));

`>`	*sol:=v(r)cos(phi)sin(theta)cos(theta)+w(r):**

Solution

`>`	u=sol;

u = (-1/14*r^2*(r1^6+r2*r1^5+r2^2*r1^4+r2^3*r1^3+r2^4*r1^2+r2^5*r1+r2^6)/(r1^4+r1^3*r2+r1^2*r2^2+r1*r2^3+r2^4)+1/14/r^3*r2^5*r1^5*(r1+r2)/(r1^4+r1^3*r2+r1^2*r2^2+r1*r2^3+r2^4)+1/14*r^4)*cos(phi)*sin(th...

Checking the Solution

`>`	*simplify(laplacian(sol,[r,theta,phi], coords=spherical))-subs(tr,xz);**

`>`	simplify(subs(r=r1,sol));

`>`	simplify(subs(r=r2,sol));

Solution in Cartesian coordinates

`>`	solve(tr,{r,theta,phi}):

`>`	allvalues(%):

`>`	subs(%[1],sol):

`>`	sol1:=simplify(%):

This is solution in Cartesian coordinates.

Checking the Solution:

`>`	*simplify(laplacian(sol1,[x,y,z])-xz);**

`>`	subs(tr,r=r1,sol1):simplify(%,radical,symbolic):combine(%,trig);

`>`	subs(tr,r=r2,sol1):simplify(%,radical,symbolic):combine(%,trig);

>

2 Example

Problem

`>`	restart;with(linalg,laplacian):with(orthopoly,P):

Solve three-dimensional interior boundary-value problem

, r < R,

u| r=R = g .

Here r=sqrt(x^2+y^2+z^2), u=u(x,y,z) or u=u(r,phi,theta).

The function g may be polynomial with variables z and x^2+y^2 or with cos(theta) .

`>`	g:=randpoly([z,x^2+y^2],degree=rand(2..5)(),terms=3);

`>`	#g:=randpoly([cos(theta)]);

>

Solving problem

`>`	n:=degree(g);

`>`	*tr:={x=rcos(phi)sin(theta),y=rsin(phi)sin(theta),z=rcos(theta)};**

`>`	subs(tr,r=R,g);

`>`	f:=simplify(expand(%));

`>`	*us:=sum(A[k]r^kP(k,cos(theta)),k=0..n);*

`>`	*k:='k': ug:=sum(A[k]r^kP[k],k=0..n);*

`>`	koef:=solve(identity(subs(r=R,us)=f,theta),{seq(A[k],k=0..n)});

$koef := {A[3] = 0, A[2] = 100/3+158/21*R^2, A[1] = 0, A[4] = -632/35, A[0] = 97+50/3*R^2+158/15*R^4}$

`>`	subs(koef,us):sol:=expand(%);

`>`	sol1:=subs(koef,ug);

In Cartesian coordinates:

`>`	solve(tr,{r,theta,phi}):allvalues(%): subs(%[1],sol): solc:=simplify(%);

Solution

Problem:

`>`	print(laplacian(u(x,y,z),[x,y,z])=0,` if r<R`);u=g,` if r`=R;

diff(u(x,y,z),`$`(x,2))+diff(u(x,y,z),`$`(y,2))+diff(u(x,y,z),`$`(z,2)) = 0, ` if r<R`

Solution:

> sol;

> sol1;

where is the nth Legendre (spherical) polynomial.

In Cartesian coordinates:

> solc;

>

Checking the Solution

`>`	simplify(laplacian(sol,[r,theta,phi],coords=spherical));

`>`	simplify(laplacian(solc,[x,y,z]));

`>`	simplify(subs(tr,solc-g),{x^2+y^2+z^2=R^2,r=R}):

`>`	combine(expand(%));

> sol1;

`>`	subs(seq(P[k]=P(k,cos(theta)),k=0..n),sol1):

`>`	simplify(%-sol);

Then sol=sol1 .

>

Note

In similar way can solve boundary-value problems in domains r>R and R1<r<R2.

3 Example

Problem

`>`	restart;with(linalg,laplacian):

Solve three-dimensional interior boundary-value problem

Delta*u+a^2*u = 0 , r < R,

| r=R = g .

Here r = sqrt(x^2+y^2+z^2) , or .

`>`	g:=cos(theta);

Solving problem

`>`	with(linalg,laplacian);

`>`	assume(a>0);

`>`	eq:=laplacian(u(r,theta,phi),[r,theta,phi],coords=spherical)+ a^2*u(r,theta,phi)=0:

`>`	*subs(u(r,theta,phi)=v(r)cos(theta),eq):**

`>`	l:=combine(%/cos(theta));

l := (2*r*diff(v(r),r)+r^2*diff(v(r),`$`(r,2))-2*v(r)+a^2*v(r)*r^2)/r^2 = 0

`>`	dsolve({l,v(0)<infinity,D(v)(R)=1},v(r));

v(r) = -R^3/(2*R*a*cos(R*a)-2*sin(R*a)+a^2*sin(R*a)*R^2)/r^2*(r*a*cos(r*a)-sin(r*a))

`>`	*sol:=subs(%,v(r)cos(theta));**

sol := -R^3/(2*R*a*cos(R*a)-2*sin(R*a)+a^2*sin(R*a)*R^2)/r^2*(r*a*cos(r*a)-sin(r*a))*cos(theta)

Solution

Problem:

`>`	*Deltau+a^2u=0,` if r<R`;u=g,` if r`=R;*

Solution:

`>`	u=sol;

u = -R^3/(2*R*a*cos(R*a)-2*sin(R*a)+a^2*sin(R*a)*R^2)/r^2*(r*a*cos(r*a)-sin(r*a))*cos(theta)

>

Checking the Solution

`>`	simplify(subs(u(r,theta,phi)=sol,eq));

`>`	diff(sol,r):

`>`	simplify(subs(r=R,%));

>

4 Example

Problem

u| r=1 =f .

We use package orbitals from share libary R5 or R4.

`>`	restart;with(share):readshare(orbitals,science): with(orbitals):

`>`	assume(theta,real);assume(phi,real);

See ?share and ?share,contents for information about the share library

`>`	#f:=randpoly([x,y,z],terms=3);

`>`	*f:=xyz;*

Solving problem

`>`	*tr:={x=rcos(phi)sin(theta),y=rsin(phi)sin(theta),z=rcos(theta)}:**

`>`	fs:=subs(tr,r=1,f);

`>`	d:=degree(f);R:=1:

We search solution in the form:

`>`	*u=Sum(Sum('a(i,j)(r/R)^iY(i,j)','j'=-'i'..'i'),'i'=0..d);*

u = Sum(Sum(a(i,j)*(r/R)^i*Y(i,j),j = -i .. i),i = 0 .. 3)

Here is from orbitals .

`>`	value(rhs(%));

`>`	Y:=(l,m)->ComplexSurfaceHarmonic(l,m,theta,phi):

We compute coefficients:

`>`	*a:=(i,j)->Int(Int(fsconjugate(Y(i,j))sin(theta),theta=0..Pi),phi=0..2Pi):**

`>`	convert(%%%,trig):value(%):

`>`	evalc(%):

`>`	expand(%,trig):

> %;

Solution:

`>`	sol:=map(simplify,collect(%,r));

>

Checking the Solution

`>`	linalg[laplacian](sol,[r,theta,phi],coords=spherical):combine(%,trig);

`>`	simplify(fs-subs(r=R,sol)):combine(%,trig);

>

5 Example

Problem

=0, r<1,

u| r=1 = cos(2*phi+Pi/3)*sin(theta)^2 .

Solving problem

`>`	restart;

`>`	Order:=2;

We search solution in the form:

`>`	Sum(r^n(A[0,n]LP(n,0)+ Sum((A[k,n]cos(kphi)+B[k,n]sin(kphi))*LP(n,k),k=1..n)),n=0..Order);

Sum(r^n*(A[0,n]*LP(n,0)+Sum((A[k,n]*cos(k*phi)+B[k,n]*sin(k*phi))*LP(n,k),k = 1 .. n)),n = 0 .. 2)

Here LP(n,k)=P(n, k, cos(theta)), P(n, k, x) is associated Legendre function. Function LP we define later.

`>`	so:=value(%);

`>`	fc:=combine(subs(r=1,so));

`>`	*f:=cos(2phi+Pi/3)sin(theta)^2;*

`>`	applyop(expand,1,f);

`>`	applyop(combine,1,%);

`>`	*collect(%,[cos(2phi),sin(2phi)]);*

`>`	dp:=collect(%,[sin(phi),cos(phi)]);

`>`	*sk:=collect(fc-dp,[sin(phi),cos(phi),cos(2phi),sin(2phi)]);*

`>`	LP:=proc(n::nonnegint,k::nonnegint) local i; orthopoly[P](n,t); for i to k do diff(%,t) end do; sin(theta)^k*subs(t=cos(theta),%); RETURN(%); end;

> sk;

Next we compute coefficients:

`>`	K:={}:

`>`	coeff(sk,sin(phi));

`>`	indets(%,indexed);

`>`	solve(identity(%%,theta),%);K:=K union %;

`>`	coeff(sk,cos(phi));

`>`	indets(%,indexed);

`>`	solve(identity(%%,theta),%);K:=K union %;

`>`	*coeff(sk,sin(2phi));**

`>`	indets(%,indexed);

`>`	solve(identity(%%,theta),%);K:=K union %;

`>`	*coeff(sk,cos(2phi));**

`>`	indets(%,indexed);

`>`	solve(identity(%%,theta),%);K:=K union %;

`>`	*subs(sin(phi)=0,cos(phi)=0,cos(2phi)=0,sin(2phi)=0,sk);*

`>`	indets(%,indexed);

`>`	solve(identity(%%,theta),%);K:=K union %;

`>`	solution:=subs(K,so);

Other form of the solution:

`>`	*sol:=r^2cos(2phi+Pi/3)sin(theta)^2;**

`>`	expand(solution-sol);

>

Checking the Solution

`>`	combine(linalg[laplacian](sol,[r,theta,phi],coords=spherical));

`>`	f-subs(r=1,sol);

>

6 Example

problem

| r=R =

theory

u = Sum(r^n*(A[0,n]*LP(n,0)/2+Sum((A[k,n]*cos(k*phi)+B[k,n]*sin(k*phi))*LP(n,k),k = 1 .. n)),n = 0 .. infinity)

where

A[k,n] = (2*n+1)*(n-k)!/2/(n+k)!/Pi/(R^n)*int(int(f*LP(n,k)*cos(k*phi)*sin(theta),phi = -Pi .. Pi),theta = 0 .. Pi)

B[k,n] = (2*n+1)*(n-k)!*int(int(f*LP(n,k)*sin(k*phi)*sin(theta),phi = -Pi .. Pi),theta = 0 .. Pi)/(2*(n+k)!*Pi*R^n)

, is associated Legendre function.

procedures

`>`	restart;

Procedure for computing LP(n,k)=P(n, k, cos(theta)), here P(n, k, x) is associated Legendre function.

`>`	LP:=proc(n::nonnegint,k::nonnegint) local i; orthopoly[P](n,t); for i to k do diff(%,t) end do; sin(theta)^k*subs(t=cos(theta),%); RETURN(%); end:

Procedure for solving problem , | r=R = :

> lap3di:=proc(f,R)
local A,B,so,lp,N;
global LP;N:=degree(f);
Sum(r^n*(A(0,n)*lp(n,0)/2+
Sum((A(k,n)*cos(k*phi)+B(k,n)*sin(k*phi))*lp(n,k),k=1..n)),n=0..N);
so:=value(%):
lp:=LP;
A:=(k,m)->(2*m+1)*(m-k)!/2/(m+k)!/Pi/R^m*
int(int(f*LP(m,k)*cos(k*phi)*sin(theta),phi=-Pi..Pi),theta=0..Pi);
B:=(k,m)->(2*m+1)*(m-k)!/2/(m+k)!/Pi/R^m*
int(int(f*LP(m,k)*sin(k*phi)*sin(theta),phi=-Pi..Pi),theta=0..Pi);
RETURN(map(factor,value(so)));
end:

examples

Example 1 ([vlad],16.20.3.)

`>`	*f:=sin(theta)(sin(phi)+sin(theta));R:=1;**

`>`	sol:=lap3di(f,R);

Checking the Solution:

`>`	simplify(linalg[laplacian](sol,[r,theta,phi],coords=spherical));

`>`	simplify(expand(f-subs(r=R,sol)),trig);

Example 2 ([vlad],16.21.2.)

`>`	*f:=sin(3phi+Pi/4)sin(theta)^3;R:='R':*

`>`	sol:=lap3di(f,R);

sol := 1/2*r^3*2^(1/2)*(cos(3*phi)+sin(3*phi))/R^3*sin(theta)^3

Checking the Solution:

`>`	simplify(linalg[laplacian](sol,[r,theta,phi],coords=spherical));

`>`	simplify(expand(f-subs(r=R,sol)),trig);

Example 3

`>`	f:=randpoly([x,y,z],degree=2);R:=3;

`>`	*tr:={x=rcos(phi)sin(theta),y=rsin(phi)sin(theta),z=rcos(theta)}:**

`>`	f:=subs(tr,r=R,f);

`>`	sol:=lap3di(f,R);

Checking the Solution:

`>`	simplify(linalg[laplacian](sol,[r,theta,phi],coords=spherical));

`>`	simplify(expand(f-subs(r=R,sol)),trig);

>

7 Example

problem

[vlad] 16.21.4.

( + )| r=R = sin(theta)^2*(sqrt(2)*cos(2*phi+Pi/4)+2*cos(phi)^2)

solving problem

`>`	restart;

`>`	*f:=sin(theta)^2(sqrt(2)cos(2phi+Pi/4)+2cos(phi)^2);*

We express in Fourier series.

`>`	select(has,f,phi);

`>`	Order:=degree(%);

`>`	fseries:=proc(f)

`>`	local k,S;

`>`	S:=int(f,phi=-Pi..Pi)/2/Pi;

`>`	for k from 1 to Order do S:=S,

`>`	*simplify(int(fcos(kphi),phi=-Pi..Pi)/Pi)cos(kphi),*

`>`	*simplify(int(fsin(kphi),phi=-Pi..Pi)/Pi)sin(kphi);od;*

`>`	RETURN(sum(S[m],m=1..nops([S])));

> end:

>

`>`	ff:=fseries(f);

`>`	simplify(expand(ff-f),trig);

Then ff=f.

We search solution in the form:

`>`	Sum(r^n(A[0,n]LP(n,0)+ Sum((A[k,n]cos(kphi)+B[k,n]sin(kphi))*LP(n,k),k=1..n)),n=0..Order);

Sum(r^n*(A[0,n]*LP(n,0)+Sum((A[k,n]*cos(k*phi)+B[k,n]*sin(k*phi))*LP(n,k),k = 1 .. n)),n = 0 .. 2)

`>`	so:=value(%);

We define LP(n,k)=P(n, k, cos(theta)), here P(n, k, x) is associated Legendre function.

`>`	LP:=proc(n::nonnegint,k::nonnegint) local i; orthopoly[P](n,t); for i to k do diff(%,t) end do; sin(theta)^k*subs(t=cos(theta),%); RETURN(%); end:

`>`	subs(r=R,so+diff(so,r))-ff:

`>`	*T:=collect(%,[seq(cos(kphi),k=1..Order),seq(sin(kphi),k=1..Order)]);*

We need from identity T=0 compute coefficients.

`>`	K:={}:

`>`	for k to Order do

`>`	*coeff(T,sin(kphi));**

`>`	indets(%,indexed);

`>`	solve(identity(%%=0,theta),%);

`>`	K:=K union %; od:

`>`	for k to Order do

`>`	*coeff(T,cos(kphi));**

`>`	indets(%,indexed);

`>`	solve(identity(%%=0,theta),%);

`>`	K:=K union %; od:

`>`	*subs(seq(cos(kphi)=0,k=1..Order),seq(sin(kphi)=0,k=1..Order),T):*

`>`	indets(%,indexed):

`>`	solve(identity(%%=0,theta),%):

`>`	K:=K union %;

Solution:

`>`	sol:=subs(K,so);

sol := 2/3+r^2*(-2/3/R/(2+R)*(3/2*cos(theta)^2-1/2)+3*(2/3/R/(2+R)*cos(2*phi)-1/3*1/R/(2+R)*sin(2*phi))*sin(theta)^2)

>

Checking the Solution

`>`	simplify(linalg[laplacian](sol,[r,theta,phi],coords=spherical));

`>`	simplify(expand(f-subs(r=R,sol+diff(sol,r))),trig);

>

8 Example

problem

[vlad] 16.24

| r=R1 = , | r=R2 =

solving problem

`>`	restart;

`>`	R1:=1;R2:=2;

> #f1:=sin(theta)*sin(phi);f2:=0;#1
#f1:=3*sin(theta)^2*sin(2*phi);f2:=3*cos(theta);#2
#f1:=7*sin(theta)*cos(phi);f2:=7*cos(theta);#3
#f1:=sin(theta)^2*(3-sin(2*phi));f2:=4*f1;#4
#f1:=12*sin(theta)*cos(theta/2)^2*cos(phi);f2:=0;#5
#f1:=sin(2*phi)*sin(theta)^2;f2:=cos(2*phi)*sin(theta)^2;#6
#f1:=cos(phi)*sin(2*theta);f2:=sin(phi)*sin(2*theta);#7
#f1:=31*sin(2*theta)*sin(phi);f2:=31*sin(theta)^2*cos(2*phi);#8
f1:=cos(theta);f2:=cos(phi)*(12*sin(theta)-15*sin(theta)^3);#9

We express 1 and f2 in Fourier series.

`>`	Order:=max(degree(expand(f1)),degree(expand(f2)));

`>`	fseries:=proc(f)

`>`	local k,S;

`>`	S:=int(f,phi=-Pi..Pi)/2/Pi;

`>`	for k from 1 to Order do S:=S,

`>`	*simplify(int(fcos(kphi),phi=-Pi..Pi)/Pi)cos(kphi),*

`>`	simplify(int(fsin(kphi),phi=-Pi..Pi)/Pi)sin(kphi); end do;

`>`	RETURN(sum(S[m],m=1..nops([S])));

`>`	end proc:

`>`	ff1:=fseries(f1);

`>`	ff2:=fseries(f2);

We search solution in the form:

`>`	Sum(r^n(A[0,n]LP(n,0)+ Sum((A[k,n]cos(kphi)+B[k,n]sin(kphi))LP(n,k),k=1..n)),n=0..Order)+ Sum(r^(-n-1)(C[0,n]LP(n,0)+ Sum((C[k,n]cos(kphi)+D[k,n]sin(kphi))LP(n,k),k=1..n)),n=0..Order);

Sum(r^n*(A[0,n]*LP(n,0)+Sum((A[k,n]*cos(k*phi)+B[k,n]*sin(k*phi))*LP(n,k),k = 1 .. n)),n = 0 .. 4)+Sum(r^(-n-1)*(C[0,n]*LP(n,0)+Sum((C[k,n]*cos(k*phi)+D[k,n]*sin(k*phi))*LP(n,k),k = 1 .. n)),n = 0 .. 4...

`>`	so:=value(%):

We define LP(n,k)=P(n, k, cos(theta)), here P(n, k, x) is associated Legendre function.

`>`	LP:=proc(n::nonnegint,k::nonnegint) local i; orthopoly[P](n,t); for i to k do diff(%,t) end do; sin(theta)^k*subs(t=cos(theta),%); RETURN(%); end:

`>`	subs(r=R1,so-ff1):

`>`	*T1:=collect(%,[seq(cos(kphi),k=1..Order),seq(sin(kphi),k=1..Order)]):*

`>`	subs(r=R2,so-ff2):

`>`	*T2:=collect(%,[seq(cos(kphi),k=1..Order),seq(sin(kphi),k=1..Order)]):*

Next we compute coefficients from identities T1=0 and T2=0 .

`>`	K:={}:

`>`	for k to Order do

`>`	*coeff(T1,sin(kphi));**

`>`	indets(%,indexed);

`>`	solve(identity(%%=0,theta),%);

`>`	K:=K union %; od:

`>`	for k to Order do

`>`	*coeff(T1,cos(kphi));**

`>`	indets(%,indexed);

`>`	solve(identity(%%=0,theta),%);

`>`	K:=K union %; od:

`>`	*subs(seq(cos(kphi)=0,k=1..Order),seq(sin(kphi)=0,k=1..Order),T1):*

`>`	indets(%,indexed):

`>`	solve(identity(%%=0,theta),%):

`>`	K1:=K union %:

`>`	T:=subs(K1,T2):

`>`	K:={}:

`>`	for k to Order do

`>`	*coeff(T,sin(kphi));**

`>`	indets(%,indexed);

`>`	solve(identity(%%=0,theta),%);

`>`	K:=K union %; od:

`>`	for k to Order do

`>`	*coeff(T,cos(kphi));**

`>`	indets(%,indexed);

`>`	solve(identity(%%=0,theta),%);

`>`	K:=K union %; od:

`>`	*subs(seq(cos(kphi)=0,k=1..Order),seq(sin(kphi)=0,k=1..Order),T):*

`>`	indets(%,indexed):

`>`	solve(identity(%%=0,theta),%):

`>`	K:=K union %;

Solution:

`>`	sol:=map(normal,subs(K1,K,so));

sol := -1/7*r*cos(theta)+48/127*r^3*cos(phi)*sin(theta)*(5*cos(theta)^2-1)+8/7/r^2*cos(theta)-48/127*cos(phi)*sin(theta)*(5*cos(theta)^2-1)/r^4

>

Checking the Solution

`>`	simplify(linalg[laplacian](sol,[r,theta,phi],coords=spherical));

`>`	combine(expand(f1-subs(r=R1,sol)),trig);

`>`	combine(expand(f2-subs(r=R2,sol)),trig);

>

While every effort has been made to validate the solutions in this worksheet, Waterloo Maple Inc. and the contributors are not responsible for any errors contained and are not liable for any damages resulting from the use of this material.

Back to contents

Maple

Maple Add-Ons

Student Success Platform

Improving Retention Rates

MapleSim

MapleSim Add-Ons

Systems Engineering

Consulting Services

Maple T.A. and Möbius

Education

Industries

Automotive and Aerospace

Robotics

Machine Design & Industrial Automation

Other

Application Areas

Product Pricing

Purchasing

Institutional Student Licensing

Maplesoft Elite Maintenance (EMP)

Support

Product Training

Online Product Help

Webinars & Events

Publications

Content Hubs

Examples & Applications

Community

About Maplesoft

Media Center

User Community

Contact

App Preview:

Boundary value problems for 3-D elliptic equations (updated to Maple 7)

Maple

Powerful math software that is easy to use

Maple Add-Ons

Student Success Platform

Improving Retention Rates

MapleSim

Advanced System Level Modeling

MapleSim Add-Ons

Systems Engineering

Consulting Services

Maple T.A. and Möbius

Education

Industries

Automotive and Aerospace

Robotics

Machine Design & Industrial Automation

Other

Application Areas

Product Pricing

Purchasing

Institutional Student Licensing

Maplesoft Elite Maintenance (EMP)

Support

Product Training

Online Product Help

Webinars & Events

Publications

Content Hubs

Examples & Applications

Community

About Maplesoft

Media Center

User Community

Contact

App Preview:

Boundary value problems for 3-D elliptic equations (updated to Maple 7)