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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
if nargin==3
    eps=1.0e-8
* g0 }6 @! e% {# {( ?, d6 C' b* uelseif nargin<3
  n& H' ]( N# A, E3 L    'error'
    return
end

/ M4 w( n# ^% `: b; k' p9 X%定义内部节点矩阵u0
h=1;k=1;
) z/ d5 V+ N% [6 ^1 }3 d; ^x=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
. @- ~, w* G+ c# h& wu=zeros(e,f);
u0=u;
0 m2 {+ ]" M- |: f( J9 p% G: @+ c
  |2 J# A2 E$ _1 Q9 v%定义外部节点p
p=zeros(e+2,e+2);
, L! v' H7 `# F) V' }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;

%定义系数矩阵A
7 `" N2 K; o' g" r* l1 }+ J3 DA=zeros(e*f,e*f);
4 A1 Z5 h" @4 z3 r1 rB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
! g: ~- j. q! I' PM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
B{1,1}=M;B{e,e}=M
( Z9 q7 F& m# T" F, wfor i=2:e-1
, i% D6 t( c$ Z$ I    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
# v2 ^& p6 v7 \, ^3 rend
, @9 i3 C9 T1 R/ p2 I! T; ~, S+ mA=cell2mat(B);  
这里总是显示
??? function J=jac(A,b,u0,eps)
    |
Error: Function definitions are not permitted at the prompt or in scripts.
9 Q; S- \, r8 U* A8 z
%定义b
( o: S/ K' L( G$ @b=ones(e*f,1);
  {0 |9 r- S4 y6 {4 Ffor 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)
6 ~2 }" T, d+ B    end
end
* y; i4 \3 k/ \& e%运用Jacobi迭代法计算
+ M$ k* n+ p2 I: i6 E% t3 ?D=diag(diag(A));
! k1 w+ B6 J. n/ m. X6 I4 _; bD=inv(D);
: y9 V  ^7 _. j% o9 m4 eL=tril(A,-1);
: }5 |; \5 \5 G5 k, O7 OU=triu(A,1);
- D: `7 V  u6 i' L3 F# Q( xB=-D*(L+U);
6 W, Z( B  X% k& qf=D*b;
, a. z, C) W* S5 uJ=B*u0+f;
while norm(J-u0)>=eps
. `# U, ~% ?7 K8 {- c% q6 C2 P6 jx0=J;
J=B*u0+f;
9 K; H, N, p1 h  D7 `2 }end
return

舒米牛牛

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

BenCam

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

madio

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