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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
. A+ T" W0 M4 U7 zfunction J=jac(A,b,u0,eps)
* U& T. J* ~8 J0 i' }if nargin==3
    eps=1.0e-8
1 `2 j* X$ O5 d9 a( W& n# e* \elseif nargin<3
  I' _1 \) y6 R, E: p7 u    'error'
7 X+ ^0 K6 L+ ^7 h/ o0 n    return
, W6 w, h7 x/ t- Xend
0 f5 u- ?0 S9 H" R
%定义内部节点矩阵u0
7 j. k; d# r4 v) Q: W* Y4 Uh=1;k=1;
! E9 h6 @1 L- ?. nx=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
) \+ z2 Y! F! ~9 b% e( S0 r8 E% Du0=u;

%定义外部节点p
p=zeros(e+2,e+2);
, W4 D4 L, R: [+ Yp(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
A=zeros(e*f,e*f);
! P$ f- [/ o5 j/ ^B=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
. _) Y( x& B" f6 y) }* ^. rM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
6 W& u" e! y' E5 Q, f' S8 W9 n% XN=-eye(e);
( Z) T) L3 }3 q9 sB{1,1}=M;B{e,e}=M
6 Z2 o% y/ Z3 g9 d$ K$ x& ^for i=2:e-1
3 B8 t% h8 H  o3 E) P4 m    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
3 p8 c' ]3 ]! k: ]  R. dend
) {* f! ~2 a: IA=cell2mat(B);  
6 N7 }9 `4 S6 r; }5 v这里总是显示
??? function J=jac(A,b,u0,eps)
% z! E  E$ n" |) `! X    |
Error: Function definitions are not permitted at the prompt or in scripts.

%定义b
+ k+ I: i( t$ w0 M, p+ F3 j8 A" Fb=ones(e*f,1);
! _% P2 u0 ~$ N0 a  sfor i=1:e
    for j=1:f
: \# m8 C6 V  A        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
' Q3 |7 w3 d9 p1 \: o* E& s3 A: j    end
end
* B: X0 i* G# ?2 O7 N) O%运用Jacobi迭代法计算
D=diag(diag(A));
$ S5 i. w1 ?6 y: OD=inv(D);
' j5 k4 i: g( QL=tril(A,-1);
7 w; F. x0 m* [, D1 |U=triu(A,1);
B=-D*(L+U);
f=D*b;
0 F1 W% K; q3 X% lJ=B*u0+f;
) M( [1 a1 u# w' k9 k' X% Mwhile norm(J-u0)>=eps
7 w, o# r: [8 u8 l* O! F  {7 j0 fx0=J;
J=B*u0+f;
end
/ t) h, g# f! o$ N0 }% D, c, Breturn

舒米牛牛

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

BenCam

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

madio

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