求助：关于分块矩阵的还原，为什么实行不了

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
9 ~# e  @, m7 B! e& Sfunction J=jac(A,b,u0,eps)
; {0 y6 b7 c3 ^* }! U, a! `% m+ `$ uif nargin==3
9 R7 n  _( K! _    eps=1.0e-8
+ p6 ~2 s2 z0 q/ celseif nargin<3
0 f; D+ `8 h* r. Y6 n    'error'
    return
- I5 _0 O5 D) v$ u3 hend
0 Y8 J/ j7 K# Z8 b& ?% T6 A
5 w5 p: b) T5 y: g2 q) m) Z%定义内部节点矩阵u0
h=1;k=1;
- u. A7 A# j, Z5 `% R+ E" i( M! ex=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
. N7 W" g( _6 su=zeros(e,f);
u0=u;
! Y% Q5 q& ~3 X) Y; h. `  m' E! X
%定义外部节点p
p=zeros(e+2,e+2);
p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
6 P* W1 l! q7 o" ^9 m4 mp(1:e+2,1)=100;p(1:e+2,f+2)=100;

%定义系数矩阵A
& N$ ~/ c- [8 y5 WA=zeros(e*f,e*f);
B=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
; M, Y: t1 `/ \6 {/ W) iM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
B{1,1}=M;B{e,e}=M
: Y& s9 ^; a6 Q6 u/ T0 G; f( v; p% Ufor i=2:e-1
* j4 x& Y. D8 @: e% @) T* l  p    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
A=cell2mat(B);  
这里总是显示
# u9 l' @' j' h5 }. j1 I??? function J=jac(A,b,u0,eps)
    |
/ i( ^- v) U; Q$ U2 a% _+ k/ eError: Function definitions are not permitted at the prompt or in scripts.
( P: t: }% s$ F" i" g% @  K( N
%定义b
b=ones(e*f,1);
for i=1:e
    for j=1:f
" e( F9 M& j' h5 q. `        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
! J) ]# f- [' C    end
end
%运用Jacobi迭代法计算
D=diag(diag(A));
D=inv(D);
6 Y) O! x. x! O, N3 t2 LL=tril(A,-1);
( O: r, _) G7 D3 R5 G7 qU=triu(A,1);
B=-D*(L+U);
f=D*b;
5 ?. s( _5 V- k3 O8 M. jJ=B*u0+f;
while norm(J-u0)>=eps
6 E( B  j" l6 q5 \6 Kx0=J;
J=B*u0+f;
0 t! t# }7 x$ o. Yend
" x( d) u4 n4 w, L4 Wreturn

舒米牛牛

自己顶一下，拜托哪位高手指点一下，纠结这个矩阵的还原，想了好多方法还是不行~~实在想不出哪里出错了

BenCam

对不起,我也不知道,帮不了忙!

madio

你是不是函数的定义没有放在M文件中？