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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
2 N6 d) j# j$ ?, O8 L( o; B+ efunction J=jac(A,b,u0,eps)
7 h- _: d: P; Z. {' P# ]8 ?if nargin==3
    eps=1.0e-8
elseif nargin<3
    'error'
    return
2 m  j$ H2 u1 x9 eend

7 q- @2 t' f! \8 L%定义内部节点矩阵u0
2 E, ?) d- o- ?  l& Q& Sh=1;k=1;
8 F# x3 s/ N: |4 f, v9 fx=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
u0=u;

%定义外部节点p
p=zeros(e+2,e+2);
: `- o4 l# ^2 {8 F9 Mp(1,1:f+2)=0;p(e+2,1:f+2)=0;      
1 t  D% \/ ]! ]- T4 w3 a0 zp(1:e+2,1)=100;p(1:e+2,f+2)=100;

$ O4 O6 O  {4 z0 W3 W" d4 K%定义系数矩阵A
# a" }1 I7 L2 k. GA=zeros(e*f,e*f);
) t. I9 P! p* O4 i8 o7 u% I  lB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
B{1,1}=M;B{e,e}=M
for i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
5 R6 z' C  y- b% G9 J& J# @A=cell2mat(B);  
3 v3 D  Q9 V" a/ x3 E& _! e6 Y, c这里总是显示
??? function J=jac(A,b,u0,eps)
    |
* G3 ~% T& Q6 ~3 x) g2 \& cError: Function definitions are not permitted at the prompt or in scripts.

1 z  ]* J6 ?- M5 J%定义b
1 r* G& F; d; H# pb=ones(e*f,1);
for i=1:e
0 e  B' i( F. q. k. \9 V    for j=1:f
        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
    end
7 n* N# o* M( J6 b: `5 f. iend
%运用Jacobi迭代法计算
D=diag(diag(A));
  t8 N* E; n. K+ Q& x2 n; B' l0 k* ED=inv(D);
L=tril(A,-1);
U=triu(A,1);
B=-D*(L+U);
  z$ H5 f5 O$ v+ O1 W+ P. `f=D*b;
4 O# o" c; O" ]4 C: C. @J=B*u0+f;
2 x2 K: m; x, ?0 ~: f6 }8 cwhile norm(J-u0)>=eps
x0=J;
5 j' f; k- C. J6 @J=B*u0+f;
end
5 d1 v, n" u' s3 p+ v, Q. sreturn

舒米牛牛

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

BenCam

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

madio

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