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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
0 N( q: D" ^' afunction J=jac(A,b,u0,eps)
if nargin==3
    eps=1.0e-8
elseif nargin<3
# `$ u" S) m1 O  R    'error'
    return
end

1 ~" M) ~. d- H% _7 q; p3 f; ?%定义内部节点矩阵u0
7 _  G: H3 g0 q, o2 eh=1;k=1;
2 z% f' G. r* ]x=0:h:17;y=0:k:10;
; z2 \/ b( n; E" u, d; N9 qe=length(x)-2;f=length(y)-2;
u=zeros(e,f);
) p) \/ X9 w) Cu0=u;

  N# u* F5 x& N  m1 a%定义外部节点p
p=zeros(e+2,e+2);
6 `6 T. d. @; B' l8 [" j. m' Tp(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;

- w$ j& U7 {% r2 p) Q%定义系数矩阵A
: ?: R# U+ U3 g* G$ Y! [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);
8 _1 M" l, x, ~8 x* x; b; u( h( _M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
2 i. @, M" W; H* lN=-eye(e);
) l) s# X* z* w1 `B{1,1}=M;B{e,e}=M
9 }; K4 Z: g/ I3 O9 ]5 `7 Lfor i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
A=cell2mat(B);  
5 E- P% \8 A! P7 m: e: k7 t. D这里总是显示
??? function J=jac(A,b,u0,eps)
8 j# `" R' c3 v    |
Error: Function definitions are not permitted at the prompt or in scripts.

%定义b
b=ones(e*f,1);
for i=1:e
    for j=1:f
  _$ T/ D- N4 n: I( v        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
( N# S# Z$ h  c& M7 V    end
end
& K. I7 r! n$ D7 o+ ~4 g: }1 E0 V%运用Jacobi迭代法计算
# `. j3 ]' v! g& f9 x% W* |D=diag(diag(A));
! h8 w& R/ p  p8 P: S  K* ^D=inv(D);
4 V( |5 N0 c; q3 G4 vL=tril(A,-1);
) ~0 c" _* a8 R2 Z) q0 qU=triu(A,1);
B=-D*(L+U);
, Z4 Z" R$ a! zf=D*b;
% w, ?2 e8 ]2 Z5 s  ?J=B*u0+f;
6 B- u- b% A' J4 D8 R. Zwhile norm(J-u0)>=eps
x0=J;
J=B*u0+f;
end
return

舒米牛牛

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

BenCam

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

madio

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