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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
# j9 w7 _) t4 k( D+ q, afunction J=jac(A,b,u0,eps)
6 p/ V5 b: Z4 _* v. o4 r% J7 Dif nargin==3
' B. n) O+ d( ?9 ]6 M    eps=1.0e-8
elseif nargin<3
    'error'
; a' q3 s. b: \    return
end

%定义内部节点矩阵u0
h=1;k=1;
x=0:h:17;y=0:k:10;
) u9 [1 L8 G8 }# t( B9 z1 Y; we=length(x)-2;f=length(y)-2;
u=zeros(e,f);
, h  C4 a; H; P* {! \: F9 m" s: Hu0=u;
7 Q# q6 l' m2 T0 f0 [; t! k
/ L2 u7 g9 ~) V, U%定义外部节点p
$ m$ p7 Z. [7 P6 K: cp=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;
8 T$ {& V( B# Z8 a5 x: F
%定义系数矩阵A
A=zeros(e*f,e*f);
9 r' u/ o+ B9 k3 l( p# Z" Y- IB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
. T3 e' X8 j3 R0 H( r8 \: Wd1=ones(e,1);d2=ones(e-1,1);
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
$ b2 d, W4 x8 yN=-eye(e);
B{1,1}=M;B{e,e}=M
for i=2:e-1
: D8 Z& `* N# J4 ]    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
  W( b( f# R% _0 nend
A=cell2mat(B);  
这里总是显示
??? function J=jac(A,b,u0,eps)
    |
Error: Function definitions are not permitted at the prompt or in scripts.
/ b  [% T2 y8 E- i- W, S* h' T
; J: V# K/ L  E9 D; i& K8 S& ]%定义b
/ |# a1 {! Y& M7 |b=ones(e*f,1);
9 Z' z3 B  A3 S- m) hfor 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)
( G$ N& c. ~3 t$ p    end
8 n& ~7 d$ N# V- x# D/ j; o4 aend
%运用Jacobi迭代法计算
( T# R4 Y/ i. Z0 wD=diag(diag(A));
/ \3 {' e  E8 d: N, J4 JD=inv(D);
L=tril(A,-1);
U=triu(A,1);
0 x$ o8 q* F% m2 D4 T$ h5 IB=-D*(L+U);
f=D*b;
J=B*u0+f;
; h* M* i# R  c+ r7 Ywhile norm(J-u0)>=eps
x0=J;
9 j$ n/ Z9 W  uJ=B*u0+f;
end
( ~3 Q4 n0 C& F/ [. ?& u0 r* vreturn

舒米牛牛

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

BenCam

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

madio

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