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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
  m) v2 W% V$ E7 ]; ifunction J=jac(A,b,u0,eps)
if nargin==3
' k3 r8 ]3 f! [8 a" b" ?' i  z7 [    eps=1.0e-8
7 M% b  N; @( E- s# celseif nargin<3
  t% B8 h8 B% }) q, l+ L6 U7 v4 K    'error'
) }$ s( T. A8 g( L( A* m% K! w    return
- v% r. r" x' \3 |" F3 x( Mend

%定义内部节点矩阵u0
2 P6 j7 ~& z8 H1 r8 y8 sh=1;k=1;
* K5 s. H+ Z" @1 n& k0 Cx=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
; ?9 Z1 d* U2 `' tu0=u;

%定义外部节点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;

3 c/ c8 h. X4 Q" p) H4 d%定义系数矩阵A
A=zeros(e*f,e*f);
) c+ w2 n8 R0 H5 RB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
$ {3 Z* p& r; E# p1 R! Od1=ones(e,1);d2=ones(e-1,1);
9 b5 S% P, ]  M2 }6 {9 i) gM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
' g4 W6 k  z% c0 M5 A' p6 [N=-eye(e);
" z7 ~6 H0 i' K& TB{1,1}=M;B{e,e}=M
! V( Y. |5 r, Ufor i=2:e-1
, `+ Q) v' I2 o1 B8 r3 f. ~    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
, m+ H1 [1 }) A4 Mend
; {4 x/ |5 |/ ~/ Y" i4 I; fA=cell2mat(B);  
这里总是显示
0 l, S: A( T) Y, y0 X6 I??? function J=jac(A,b,u0,eps)
    |
Error: Function definitions are not permitted at the prompt or in scripts.

% f3 X8 I& H; I, j* ]3 f* d%定义b
b=ones(e*f,1);
for i=1:e
    for j=1:f
/ L" H4 S. h2 y+ @        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
    end
$ x4 |0 ~8 _4 `2 Xend
%运用Jacobi迭代法计算
D=diag(diag(A));
3 U9 a3 W# X# \3 p! \D=inv(D);
L=tril(A,-1);
1 S/ g7 O& d9 l/ i: `6 }2 SU=triu(A,1);
B=-D*(L+U);
' Q! d* r+ j9 h* L" tf=D*b;
3 G$ |8 }' K6 v' U' mJ=B*u0+f;
- l# l5 R# G6 q% y: Q4 Cwhile norm(J-u0)>=eps
x0=J;
* G& R' P, W4 R, S) M$ p6 @, m9 XJ=B*u0+f;
% N# L( s1 F5 y- ]% X9 `/ G! V+ wend
+ \  G3 r$ z4 `6 ~" {4 `+ o, hreturn

舒米牛牛

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

BenCam

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

madio

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