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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
if nargin==3
' N3 m/ @* Y0 W/ c# W( o* e" s    eps=1.0e-8
# D4 L$ G/ V! C( lelseif nargin<3
    'error'
    return
end
9 G5 O8 s- K( y) a) g( {( G. `
% Z. F8 ]  Z* j! V" g%定义内部节点矩阵u0
  W  y, K8 ~2 `( R0 ?h=1;k=1;
7 ?0 `4 t; E; W% i9 @x=0:h:17;y=0:k:10;
, Q9 F7 ^& ]/ s7 R3 Ce=length(x)-2;f=length(y)-2;
* p( s/ }& b. j0 j0 L0 I5 |' k7 _# Hu=zeros(e,f);
u0=u;
* h, P! C7 s, C0 b
- O0 c7 y% p+ e%定义外部节点p
2 t0 @. a  W/ W7 R6 [p=zeros(e+2,e+2);
p(1,1:f+2)=0;p(e+2,1:f+2)=0;      
9 F0 o# N! s+ e$ a- O& @4 X) S! Tp(1:e+2,1)=100;p(1:e+2,f+2)=100;

%定义系数矩阵A
$ m; }' L/ m& L& M' J/ ?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);
5 l$ u" z- r3 N. q/ {4 M* J2 WM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
N=-eye(e);
9 @7 Z' [9 K3 aB{1,1}=M;B{e,e}=M
8 n1 `) l0 N# sfor i=2:e-1
    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
end
A=cell2mat(B);  
这里总是显示
, V7 K  s6 b) k2 \, w??? function J=jac(A,b,u0,eps)
    |
* W- y# u9 v: q- u; `( m& NError: Function definitions are not permitted at the prompt or in scripts.
% c  b$ b) ?! X. M1 X) o; `+ r
%定义b
* }1 I2 I% T8 w* @8 Rb=ones(e*f,1);
" `8 T8 q; [7 \6 [) G+ ]! u* ufor i=1:e
    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)
    end
end
%运用Jacobi迭代法计算
D=diag(diag(A));
D=inv(D);
L=tril(A,-1);
U=triu(A,1);
B=-D*(L+U);
f=D*b;
J=B*u0+f;
while norm(J-u0)>=eps
x0=J;
J=B*u0+f;
end
1 f1 f& w# J* L+ j  i, h/ }2 mreturn

舒米牛牛

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

BenCam

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

madio

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