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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
if nargin==3
' b  M! G0 [9 x) b6 k    eps=1.0e-8
elseif nargin<3
    'error'
    return
7 ?* o1 O2 m( `$ `) B" ^8 f/ vend
1 Q* m4 Q; f+ G# I) c- L* [
! |0 R) P) C( [3 m( h$ N%定义内部节点矩阵u0
h=1;k=1;
x=0:h:17;y=0:k:10;
' {: R4 U" A0 s# W$ s( n+ [e=length(x)-2;f=length(y)-2;
, ~  k* b, s2 M* R1 u6 B: Iu=zeros(e,f);
  x: Q+ q* l) `. `/ x( ]u0=u;
- ?* h; u" Y3 \: p/ C
  V) c0 U/ s2 O' x5 f6 d%定义外部节点p
! W1 s$ i2 c# c/ C7 N6 u; [% |p=zeros(e+2,e+2);
p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
/ s2 e; D& i7 S6 ?' t0 @- Y6 Lp(1:e+2,1)=100;p(1:e+2,f+2)=100;

! @0 b! B5 c1 q! b%定义系数矩阵A
7 c4 g  |5 Y* W7 [A=zeros(e*f,e*f);
8 T3 Z0 ~( j& ~' S& }. WB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
! y! A, u  C, DM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
8 b7 Q, e  y1 o+ R6 a* RN=-eye(e);
  ]' h3 M4 w8 N  ~3 mB{1,1}=M;B{e,e}=M
for i=2:e-1
2 S6 U( i0 ~; V% u* _, I  Q+ ?    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
A=cell2mat(B);  
这里总是显示
??? function J=jac(A,b,u0,eps)
    |
Error: Function definitions are not permitted at the prompt or in scripts.

%定义b
5 x9 {6 T, @: ^( z9 l& pb=ones(e*f,1);
4 ]; u# b  k8 y& w5 f3 Cfor i=1:e
/ D8 K5 |4 g0 a* b4 L  T+ v+ _    for j=1:f
0 ?6 y! f+ k  X1 m# O) z        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
6 T' G" s) d% N$ R. t7 j    end
end
%运用Jacobi迭代法计算
D=diag(diag(A));
/ x" j; d4 \) p- K7 D0 rD=inv(D);
6 x1 k( [1 F7 D! rL=tril(A,-1);
U=triu(A,1);
B=-D*(L+U);
( V: W8 i& K% p# f6 z2 yf=D*b;
* b, _0 r7 ]3 C8 T( R! JJ=B*u0+f;
$ E# ]) F! P8 a& ^4 _9 {while norm(J-u0)>=eps
x0=J;
J=B*u0+f;
end
3 r# p  Y# S3 f8 Breturn

madio

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

BenCam

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

舒米牛牛

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