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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
if nargin==3
' t% t: \! B5 L/ z/ S    eps=1.0e-8
( [! P& p- `1 ^3 qelseif nargin<3
, k7 P# G8 D+ ]% a! [    'error'
* u1 Z7 Q* O& V% q- z    return
" m& c! I# U" ^2 f9 n% Uend

%定义内部节点矩阵u0
h=1;k=1;
5 C' z3 L4 o: [) c! ?/ x6 D: ?+ ~x=0:h:17;y=0:k:10;
. L; ~4 {7 [, R+ q+ R, ^- g' pe=length(x)-2;f=length(y)-2;
2 @6 X7 ?1 n  fu=zeros(e,f);
4 F( I$ U  D2 B0 V  O/ Lu0=u;
) r9 S0 W: R; F7 i
%定义外部节点p
/ C3 Q- [) J/ m7 ^5 P8 [p=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;
/ k  L  y/ v9 q
( a2 N7 a5 C4 |: U%定义系数矩阵A
A=zeros(e*f,e*f);
4 T: c! ~, c4 |' O" p& gB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
' [$ P( q0 a7 P2 Nd1=ones(e,1);d2=ones(e-1,1);
/ E- R# b* y: i$ Q' H, W3 R. EM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
B{1,1}=M;B{e,e}=M
for i=2:e-1
  C5 F, U$ k5 C2 e    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
4 [+ l1 [+ j% I( R9 x9 K& u6 zA=cell2mat(B);  
- d2 x8 T' |) Y' ?8 _  G! U: R这里总是显示
??? function J=jac(A,b,u0,eps)
    |
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
3 x1 [5 ^! S( Z7 O+ M1 v/ t        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
4 y8 g) f$ Q& J* U8 e    end
' p/ J: U  w. o/ s, A) x6 i& _! F5 yend
%运用Jacobi迭代法计算
. ?1 q. h; P, ~4 d" [+ L+ gD=diag(diag(A));
3 @! A; }7 i9 M9 y/ I- I# ?4 z$ S; wD=inv(D);
8 C! O+ `# u4 dL=tril(A,-1);
4 b/ r2 t! n. _' \" I" F8 wU=triu(A,1);
B=-D*(L+U);
- {0 m4 b. l! B& D8 af=D*b;
J=B*u0+f;
while norm(J-u0)>=eps
x0=J;
J=B*u0+f;
' m3 V2 {( [* X7 send
return

舒米牛牛

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

BenCam

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

madio

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