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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
if nargin==3
    eps=1.0e-8
6 I; b# C$ T/ a. Relseif nargin<3
# h* K8 C7 J. m9 ^) k# r    'error'
& k8 V; |, n$ p+ i' o    return
' P1 e7 S% u# }/ w2 ^end

%定义内部节点矩阵u0
h=1;k=1;
x=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
: f& K+ b4 y3 \# `; @u0=u;

# a# W# d% e- S# |, |%定义外部节点p
p=zeros(e+2,e+2);
" T# L$ h' M& E- O% ip(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);
B=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
/ F' ?4 q: S* O% ad1=ones(e,1);d2=ones(e-1,1);
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
B{1,1}=M;B{e,e}=M
2 ], V0 ~* b% o+ S0 D) I$ Afor i=2:e-1
3 D. a3 t* g, g    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
6 e1 G2 A$ D3 y& R3 qA=cell2mat(B);  
这里总是显示
* \1 ~% G5 q- D3 ?1 L0 y  e0 `??? function J=jac(A,b,u0,eps)
    |
Error: Function definitions are not permitted at the prompt or in scripts.
3 N1 G. }+ D1 t6 d  i) X
%定义b
( o: X" T. I. Q& Y3 q) h5 Sb=ones(e*f,1);
for i=1:e
    for j=1:f
5 T$ f# w3 S5 x5 {& M        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
    end
end
8 t% g3 H; l% r" E( e' G%运用Jacobi迭代法计算
1 L( C& D1 s( U+ }  ~% `D=diag(diag(A));
D=inv(D);
L=tril(A,-1);
. u8 }2 S0 S$ C+ O! LU=triu(A,1);
B=-D*(L+U);
' t9 V" V' K: u5 h! K2 Of=D*b;
J=B*u0+f;
while norm(J-u0)>=eps
; z0 U, u! {2 G7 n3 cx0=J;
, k, |+ R' E7 a, F1 T- |! YJ=B*u0+f;
! P# n+ i; T3 }0 \( yend
return

舒米牛牛

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

BenCam

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

madio

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