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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
) h  ]4 \3 e9 J5 `function J=jac(A,b,u0,eps)
7 [$ X( i" h  o  t+ _. J. Oif nargin==3
9 v, x+ s3 [8 v- P3 |    eps=1.0e-8
) A' S  G4 C) J; N' e* `$ z* oelseif nargin<3
7 y' y2 n: S0 v- R6 t    'error'
, `9 U! n4 c$ W! Q# T* M  P) L    return
end
' Z+ o" w& C' D" Z# H; c9 k6 |
$ u% _% K* p* f2 y$ g%定义内部节点矩阵u0
h=1;k=1;
x=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
u0=u;
) W; L) M8 \9 C0 @
%定义外部节点p
p=zeros(e+2,e+2);
p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
p(1:e+2,1)=100;p(1:e+2,f+2)=100;
% T3 m6 t: m, F9 n8 O6 j
%定义系数矩阵A
A=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);
3 e. e# S( v* u9 Z6 V2 m. a/ VM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
* _& y$ N1 ]; u5 }0 f6 PN=-eye(e);
6 U: e1 ]9 K( a  H! ^B{1,1}=M;B{e,e}=M
for i=2:e-1
. I( F% d) a! p4 ?$ |+ M    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
% J& U/ H! ^5 k4 j4 gend
2 J" M9 \# g  r6 X3 z. N% L( W* gA=cell2mat(B);  
- i7 Y' A4 k2 l4 B) m& D$ k7 Z这里总是显示
, K- \: E6 B! O0 d/ a??? function J=jac(A,b,u0,eps)
    |
$ z& D; {; c, }! w' nError: Function definitions are not permitted at the prompt or in scripts.

9 v7 Y! h/ ?; s2 D) ]%定义b
  z8 |: f/ K1 L# a9 ?b=ones(e*f,1);
for i=1:e
- q7 Q3 o. _8 j5 w) p    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)
% P7 V3 [' i4 g3 O    end
$ `3 U, P  P( D7 o, `* wend
%运用Jacobi迭代法计算
D=diag(diag(A));
D=inv(D);
7 p6 ~% v" K! j3 b5 U/ LL=tril(A,-1);
U=triu(A,1);
( U" I/ }9 m7 s& H, TB=-D*(L+U);
f=D*b;
8 t! l8 ?! t( {( }' x, g0 eJ=B*u0+f;
# ~. O6 w3 w, L/ C1 x; ^6 ]! Fwhile norm(J-u0)>=eps
x0=J;
% W* [" \- e3 y" s' V4 Q* O, C, o1 y9 ?J=B*u0+f;
' e% C6 _, S6 @end
% K; `2 w$ Z; J# X. |; L- {return

舒米牛牛

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

BenCam

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

madio

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