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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
& w( `  |1 o: pif nargin==3
/ m4 j1 z7 p: A' [- j# z6 Y, j    eps=1.0e-8
5 f# a' K5 C6 w$ e4 n1 {  M9 J3 helseif nargin<3
    'error'
& e, b( `& ^- h5 L4 C; B  t    return
, Y! m  L( E7 Q( Kend
$ [8 I! t  d, L+ C& j  ?
%定义内部节点矩阵u0
, m( u" G4 X! z2 h6 ?& M$ M( o. ~$ Eh=1;k=1;
x=0:h:17;y=0:k:10;
+ {( U/ U% K6 x2 {e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
u0=u;
& J; {3 k# @3 l) s% x6 J
%定义外部节点p
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;

/ {. E8 H. i) K  ~( e: O; ?9 N%定义系数矩阵A
A=zeros(e*f,e*f);
B=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
7 M& t' U9 J  O1 B+ iM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
5 K# {# K0 K4 i4 l& X* _( qB{1,1}=M;B{e,e}=M
) E5 e5 B, K" g. Mfor i=2:e-1
3 J" V( o* I4 H6 |    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
A=cell2mat(B);  
这里总是显示
" Y0 d8 F0 ^' V. q$ w/ _??? function J=jac(A,b,u0,eps)
7 g! S- X8 P  q3 J* @2 o    |
: }: S& }7 M9 S& gError: Function definitions are not permitted at the prompt or in scripts.
2 ?% v# I% P4 b" w7 ?
( K; B: W7 z0 h%定义b
b=ones(e*f,1);
for i=1:e
    for j=1:f
9 l! M6 D% [$ {* u/ R- k        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
    end
1 w- }4 W+ Y% ]0 Dend
%运用Jacobi迭代法计算
D=diag(diag(A));
D=inv(D);
L=tril(A,-1);
U=triu(A,1);
6 D! E( l9 r* P; MB=-D*(L+U);
f=D*b;
3 J1 }! P- U& G1 Z: I( B3 j- f, @J=B*u0+f;
, i: F8 F9 p) w$ Q6 ywhile norm(J-u0)>=eps
2 p( t, p0 S% [3 s, k! f5 Ox0=J;
: P& f' B; {: a% a# rJ=B*u0+f;
end
return

舒米牛牛

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

BenCam

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

madio

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