matlab使用pde解三维偏微分方程的程序

Olivier1

t(x,y,0)=300
t(x,y,100)=300
t(0,y,z)=t(120,y,z)
t(x,0,z)=t(x,100,z)
  0≤  x≤120
  0≤  y≤100
  0≤  z≤10[/math]

写的程序如下，不知道怎么改：

%pdefun1

function [c,f,s]=pdefun1(x,y,z,dt)

  c=[1];

  f=[-dt,-dt];

  s=0

%pdebc1

function [pa,qa,ka,pb,qb,kb]=pdebc1(xa,ya,ta,xb,yb,tb,z)

  pa=[0];

  qa=[0];

  ka=[300];

  pb=[100];

  qb=[120];

  kb=[300];

%pdeic1

function t0=pdeic1(x,y)

  u0=[300 300];

主调动程序：

clear all;clc;

x=0:1:100;

y=0:1:120;

z=0:1:100;

m=0;

sol=pdepe(m,@pdefun1,@pdeic1,@pdebc1,x,y,z);

figure('pderunwindows')

suft(x,y,z,sol(:,:,:,1)

title('the solution of t')

xlabel('x')

ylabel('y')

zlabel('z')

请求高手解答.[math](∂^2 t)/(∂x^2 )+(∂^2 t)/(∂y^2 )+(∂^2 t)/(∂z^2 )=0