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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
. S$ ~/ o6 Z( Q3 t7 r+ gfunction J=jac(A,b,u0,eps)
3 Y: Y& Y! p7 U" ]if nargin==3
    eps=1.0e-8
elseif nargin<3
    'error'
    return
end
( D% b6 U# G- c# m
%定义内部节点矩阵u0
+ [) o# T! P, f+ J  o8 ^h=1;k=1;
+ d2 ]0 @# W6 G" Z8 k3 ax=0:h:17;y=0:k:10;
# r+ r9 \! J+ O9 Ie=length(x)-2;f=length(y)-2;
7 R9 G, s" {) P( c, Su=zeros(e,f);
u0=u;
- T* i. b: }5 y! D
%定义外部节点p
+ \* B1 \; b8 p. ~% ^p=zeros(e+2,e+2);
p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
6 S0 g2 |9 c1 Op(1:e+2,1)=100;p(1:e+2,f+2)=100;
6 |" a& Y" s) Z& O% A
%定义系数矩阵A
( ]- C) D! _+ ZA=zeros(e*f,e*f);
2 D; G5 I9 K" c( OB=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);
" _' G9 o8 W( ~; G1 n7 D8 y7 g2 XN=-eye(e);
/ r2 }4 a  u3 l" \# L) N$ vB{1,1}=M;B{e,e}=M
for i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
/ @" @( n; l# B2 l' ?3 fend
' e& k; H  X1 k( N# kA=cell2mat(B);  
5 ^8 {4 @/ j) _4 y这里总是显示
+ o4 a" g' z' F2 d??? function J=jac(A,b,u0,eps)
( l+ m3 n0 Q+ [% U  e0 d$ L$ k    |
Error: Function definitions are not permitted at the prompt or in scripts.
+ H/ C1 u$ _/ @; s8 w0 H
%定义b
b=ones(e*f,1);
( K1 V6 ]# w% @' S+ Xfor i=1:e
    for j=1:f
$ |2 W' j, |% V( ]% [) q        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
% H' v( e. d% K( U$ ?2 b    end
1 H5 X( u, X0 Z8 @end
$ M* b7 k7 f8 s- k5 V/ r+ a%运用Jacobi迭代法计算
: W4 c4 o. i& GD=diag(diag(A));
+ t, r1 O. T' j/ B$ M4 H7 dD=inv(D);
$ y$ U2 |4 x  e. T: g7 LL=tril(A,-1);
! \9 r6 ?. J3 p5 ?U=triu(A,1);
$ M6 [1 r8 k- b: W( U& aB=-D*(L+U);
7 H* {  x# M  q' Cf=D*b;
J=B*u0+f;
while norm(J-u0)>=eps
* `; Q. U( X  h; u+ E! q4 Tx0=J;
J=B*u0+f;
end
  M4 g) ~( M. C- Y& |& Rreturn

舒米牛牛

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

BenCam

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

madio

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