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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
if nargin==3
8 T! F; K5 }: [! }% q2 h! s4 t    eps=1.0e-8
2 u# ]4 P7 Z( jelseif nargin<3
6 S; }( E  O( H- n& k    'error'
    return
2 K0 j/ M7 a/ K3 R! V( Pend
1 U. t& S4 ~+ W  [* I% A0 f, h
6 ^9 ?3 z3 D4 U- ]%定义内部节点矩阵u0
3 I! E& {3 S% f$ z5 P% ^: x: Oh=1;k=1;
  @; d- }! r0 G9 `$ A' ~x=0:h:17;y=0:k:10;
! D9 q9 g* _9 Q( B* _* A- Ke=length(x)-2;f=length(y)-2;
u=zeros(e,f);
% y. d/ ~+ W) Fu0=u;
+ r& n8 S' o* t+ I
%定义外部节点p
p=zeros(e+2,e+2);
p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
* l4 D! e( U, ^/ F0 B* a) Np(1:e+2,1)=100;p(1:e+2,f+2)=100;

- z& j7 W( L, ^* S" e  c+ }; l%定义系数矩阵A
A=zeros(e*f,e*f);
B=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
! b  q6 ?: w3 ?: t- Dd1=ones(e,1);d2=ones(e-1,1);
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
" E) n8 l; Y1 I: e) H& IB{1,1}=M;B{e,e}=M
9 U+ d) r5 S8 E9 y( xfor i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
( s: D# V: C9 nend
A=cell2mat(B);  
" x* ]$ W0 K% _6 x# e6 {这里总是显示
* L3 [: ^! P0 n??? function J=jac(A,b,u0,eps)
. a3 T. ^; E1 ]( X( n/ m    |
Error: Function definitions are not permitted at the prompt or in scripts.
5 m( ^! t/ C& W* n3 N( ], `
1 t$ H# o4 F( I( s2 c! b: I. y; l%定义b
# W7 q, v% ]' U4 O9 H  V2 db=ones(e*f,1);
for i=1:e
    for j=1:f
# k7 F$ n  K- W* W- [$ @5 Y        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
    end
6 o& S5 q. K1 z6 m* ~- [end
%运用Jacobi迭代法计算
D=diag(diag(A));
D=inv(D);
L=tril(A,-1);
- S3 c! ?* s* _8 GU=triu(A,1);
B=-D*(L+U);
7 Q3 Q9 d1 E. e4 Y( r+ w* Z4 U6 nf=D*b;
1 K' e* X; s  C7 v( c, F# BJ=B*u0+f;
while norm(J-u0)>=eps
x0=J;
9 d, d4 q! [% S, z1 m9 xJ=B*u0+f;
end
9 m  Z7 t( |5 D4 i2 xreturn

舒米牛牛

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

BenCam

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

madio

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