Unbalanced or unexpected parenthesis or bracket急求！！

function [c,z,t]=PDEadsorptionbed(Z,T,cbc,IZ,IT)
%边界条件：c(0,t)=cbc
%输出参数c,q-解矩阵，第一行表示初值，第一列和最后一列表示边界值，第二行表示第二层。
%z-空间变量
%t-时间变量
%输入参数：Z—空间变量z的取值上限
% T-时间变量t的取值上限
% cbc-变量c的入口浓度
% IZ-沿z轴的等分区间数
% IT-沿t轴的等分区间数
%计算步长
dz=Z/IZ;%z的步长
dt=T/IT;%t的步长
%计算初值和边界
c=zeros(IZ+1,IT+1);q=zeros(IZ+1,IT+1);p=zeros(IZ+1,IT+1);qm=zeros(IZ+1,IT+1);
for i=1:IZ+1
c(i,1)=0;q(i,1)=0;
p(i,1)=0;
qm(i,1)=0;
end
for j=1:IT+1
c(1,j)=cbc;
p(1,j)=8.314*3.23*1.09*c(1,j);
qm(1,j)=2.983*100*(0.336*p(1,j)/(1+0.336*p(1,j)+0.0000000005046*power(p(1,j),12.608)/(1+0.0000000005046*power(p(1,j),12.608))));
end
for j=1:IT
for i=2:IZ
c(i,j)=(-0.661*c(i,j+1)/dt-0.1*c(i+1,j)/dz+0.0051*(c(i+1,j)+c(i-1,j))/dz/dz-0.097*(qm(i,j)-q(i,j))/(0.0102/dz/dz-0.661/dt-0.1/dz));
p(i,j)=8.314*3.23*1.09*c(i,j);
qm(i,j)=2.983*100*(0.336*p(i,j)/(1+0.336*p(i,j)+0.0000000005046*power(p(i,j),12.608)/(1+0.0000000005046*power(p(i,j),12.608))));
q(i,j)=(q(i,j+1)-0.01*qm(i,j))/(1-0.01*dt);
end
end
return

Unbalanced or unexpected parenthesis or bracket.
这个是括号不匹配的意思，少了半边括号

你这段得问题在于
x=0:0.1:1;
plot(x,T(x,11),':*b')
x是小数，怎么能做T数组的下表呢

y=0.79*(x)^0.048(x>0&x<=150)+exp(-0.021*(x-150))(x>150&x<=270); % 这行有问题 % 0.048(x>0&x<=150)，中间是*的关系吗？直接这样写不行啊 % (-0.021*(x-150))(x>150&x<=270)，括号和括号中间是不是也少了乘*号。