برنامه اصلاح شده:
% function pdexFR
global ps e uz dh N L m Kf n
close all
ps=1.1;
Kf=0.991;
n=0.782;
e=0.25;
dh=3e-3;
uz=100;
N=50;
L=20;
m=0;
x=linspace(0,L,N);
t=linspace(0,1,51);
sol = pdepe(m, @pdeFRENDLICH, @icFRENDLICH, @bcFRENDLICH ,x ,t);
sol=real(sol);
u = sol(:,:,1);
surf(x,t,u)...