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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
! Z! x  Z. S3 o+ l& m; g$ h9 ]function J=jac(A,b,u0,eps)
if nargin==3
) j8 o0 i  |) d* u% t    eps=1.0e-8
0 e( l" z  p/ j; R0 a1 Belseif nargin<3
    'error'
    return
6 a/ ]5 T& Y" n0 a, S4 H+ `! Bend
8 p5 y2 V$ c: c9 p8 h4 h
%定义内部节点矩阵u0
$ F6 U$ {- D! b' ?h=1;k=1;
x=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
" B- R" \4 r6 ^u=zeros(e,f);
u0=u;
: ]+ y' f2 c  R0 Y" B4 G
# N2 |) n$ B: L+ X4 i( E%定义外部节点p
) y  T+ @. P$ L/ A2 ~9 r5 Wp=zeros(e+2,e+2);
2 }" ^3 D2 G% k4 I5 T/ Op(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;
/ D: M' ~; j6 K# D  J
+ A; f9 [/ O, ^4 C4 q$ `1 ^%定义系数矩阵A
A=zeros(e*f,e*f);
5 c, a7 m3 q3 _: d! A3 S! OB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
5 |: ~) ?! S8 {M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
! n, a8 _! `& Q  ~3 TN=-eye(e);
B{1,1}=M;B{e,e}=M
, c! D' i8 T" [. e  g' H8 b% n& dfor i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
" \$ s6 H5 i" s. ?( H( t9 lend
% w0 I2 k3 |' i* a) M' u7 jA=cell2mat(B);  
这里总是显示
??? function J=jac(A,b,u0,eps)
    |
' u- v: Y( C- z$ C# B4 Q$ JError: Function definitions are not permitted at the prompt or in scripts.
8 ]. i$ B* G1 I, F- y+ f  R4 P( i
%定义b
3 h4 t, m( m, u% V0 D% Y: \b=ones(e*f,1);
for i=1:e
    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)
* H+ S- n: }  E/ x% V) i    end
) J3 M$ |# [" `) c+ a. Yend
1 a! r4 `/ a( H/ A1 I& D%运用Jacobi迭代法计算
D=diag(diag(A));
# r8 n# X, e$ _D=inv(D);
1 }3 P( u6 k1 c* K! k8 nL=tril(A,-1);
U=triu(A,1);
B=-D*(L+U);
$ q1 K5 m9 m* @. @# rf=D*b;
J=B*u0+f;
9 A5 z, t3 r( M, Pwhile norm(J-u0)>=eps
9 D0 Q- z; d3 o; x/ Hx0=J;
J=B*u0+f;
0 N! L6 O/ }. N4 gend
return

舒米牛牛

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

BenCam

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

madio

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