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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
$ K5 p7 q- Z; |/ p0 |function J=jac(A,b,u0,eps)
; o/ M& J4 H8 P+ ]if nargin==3
    eps=1.0e-8
! c& z& z4 U$ j! p3 Helseif nargin<3
; M. s2 \% c3 C& D! r0 d" p    'error'
; u0 w2 ?& A. P7 M5 ^, b5 h+ _    return
4 D/ m, j3 d7 t* v5 I1 l1 kend

+ P/ b9 W1 ]2 y5 L  x; N%定义内部节点矩阵u0
h=1;k=1;
x=0:h:17;y=0:k:10;
% q5 W3 O" b2 R/ m, w0 _e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
& u4 _5 f7 f# o8 H' ]1 n/ Zu0=u;

%定义外部节点p
% j. Y* ]: m9 J$ S; C. g0 X6 D4 ~. Pp=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;

%定义系数矩阵A
/ u0 n8 H7 Q2 x7 i9 VA=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);
M=4*diag(d1)-diag(d2,1)-diag(d2,-1);
& e4 h9 E: g8 C* t0 I9 F( uN=-eye(e);
, @3 l$ h7 ^& f- NB{1,1}=M;B{e,e}=M
for i=2:e-1
9 m# {& ?' P, M6 P2 R) x6 N7 ~    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
$ R. v$ ]' F" Y3 p9 Z; t  N0 e& Qend
0 A7 u& |  z1 X# t' f% l, SA=cell2mat(B);  
这里总是显示
??? function J=jac(A,b,u0,eps)
; ]% x& U. k% T/ {    |
Error: Function definitions are not permitted at the prompt or in scripts.

+ a  Y0 A) x7 @' ?: v$ z%定义b
$ Q* w# V: `# |3 s5 bb=ones(e*f,1);
: `* w/ a5 p$ f+ p# m7 u1 pfor i=1:e
1 c7 K8 A! j8 V" u; f, G4 \    for j=1:f
        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
, i* B5 H/ T1 P    end
6 y5 Z# z1 q" }! q" j- N6 X' U% y8 Jend
%运用Jacobi迭代法计算
D=diag(diag(A));
0 G8 F; H4 c* I2 K, ^D=inv(D);
: |9 c$ j& u( y* D9 m2 c' x6 K* U7 k5 _L=tril(A,-1);
) j3 U0 D$ b9 ]5 M& G4 c( O& x3 t- `U=triu(A,1);
& p! }5 q3 N6 g  C6 YB=-D*(L+U);
( |2 b! d8 m7 N8 S  Qf=D*b;
/ F# b+ |% c5 jJ=B*u0+f;
6 x/ P: N, Y1 x1 \7 F3 [while norm(J-u0)>=eps
$ t9 I- o* K* w8 q  t2 Gx0=J;
+ t- Z; ^8 l2 F! g7 |( [J=B*u0+f;
3 j! `  b2 v! p% k% Dend
return

舒米牛牛

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

BenCam

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

madio

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