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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
" d" {3 j7 _8 K* ?3 b1 B% u. {( Z0 tfunction J=jac(A,b,u0,eps)
if nargin==3
8 z# u+ j0 q" \    eps=1.0e-8
elseif nargin<3
    'error'
    return
1 Z/ }0 d8 c; {9 u/ bend
4 ~3 Z5 ~8 F; I+ J' q* {6 [" i
%定义内部节点矩阵u0
' d! W, c7 a4 I/ T% a- d- Ih=1;k=1;
x=0:h:17;y=0:k:10;
& V, S* r2 F9 @% s8 s" @e=length(x)-2;f=length(y)-2;
0 M+ }* g4 H& c4 A% ju=zeros(e,f);
. f# O8 J' W: D" z& d3 Au0=u;
  E& _6 X7 m, Y6 s6 a! U+ D
+ q. y4 r$ }. I: v3 V- {, s% H) H%定义外部节点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;

%定义系数矩阵A
% e: k! K* F- `* q( h5 ^A=zeros(e*f,e*f);
( @$ H7 x$ W* ]) DB=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);
" B( H$ m) K# u3 ^8 }N=-eye(e);
B{1,1}=M;B{e,e}=M
7 O5 X8 c9 R. b2 wfor i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
9 N* ~% ~- r- e; M! A9 x6 Oend
+ k/ Z+ I6 U8 cA=cell2mat(B);  
这里总是显示
??? function J=jac(A,b,u0,eps)
! W- v* T7 ~  X. M2 p5 j/ Y    |
Error: Function definitions are not permitted at the prompt or in scripts.
# F$ [7 v5 ~/ ~2 y* V
%定义b
5 H5 \$ n7 e% P3 cb=ones(e*f,1);
2 h, K9 L$ B# N4 Vfor i=1:e
    for j=1:f
, }( B3 f( C, N! l/ L        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
    end
. M# C8 B8 _9 H/ x9 l# ^end
%运用Jacobi迭代法计算
D=diag(diag(A));
D=inv(D);
L=tril(A,-1);
U=triu(A,1);
4 G, T* h$ h4 y# Y' |9 I' kB=-D*(L+U);
f=D*b;
% Z% n& S7 c; x- R8 AJ=B*u0+f;
8 \* E/ y1 k' o4 C1 M" Y& n  uwhile norm(J-u0)>=eps
x0=J;
J=B*u0+f;
& \4 I2 _% `5 p% k. a" Xend
return

舒米牛牛

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

BenCam

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

madio

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