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

舒米牛牛

%运用Jacobi迭代拟合出来的u关于x和y的矩阵
function J=jac(A,b,u0,eps)
9 Y( P! w3 e# z# ^if nargin==3
    eps=1.0e-8
elseif nargin<3
* T! J+ f! p1 y1 I* J5 p' O    'error'
9 K1 T" ~% a, R! j2 T- X    return
  p* H& @0 E/ J8 k$ ?6 I! E/ Gend
1 |6 X; m& k& k. b& P5 }1 V+ t
%定义内部节点矩阵u0
h=1;k=1;
x=0:h:17;y=0:k:10;
e=length(x)-2;f=length(y)-2;
u=zeros(e,f);
u0=u;

%定义外部节点p
0 R7 W* w0 R7 q! Q$ K. Qp=zeros(e+2,e+2);
  a5 ?4 m, g" D7 x/ Mp(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;
7 ~/ f' J' C# m- G
$ @2 t4 w( A6 _6 D%定义系数矩阵A
A=zeros(e*f,e*f);
% s4 Q4 R/ F) ]4 v& X' YB=mat2cell(A,ones(e*f/e,1)*e,ones(e*f/e,1)*e);
d1=ones(e,1);d2=ones(e-1,1);
9 t0 Q# W8 @' e+ D( E: fM=4*diag(d1)-diag(d2,1)-diag(d2,-1);
; T6 A0 c# P# gN=-eye(e);
  f4 Z" z4 j# w! jB{1,1}=M;B{e,e}=M
% c% y7 w- \( [8 ]- Hfor i=2:e-1
7 `* f$ n- `- n2 }$ e9 p! ~    B{i,i}=M;B{i-1,i}=N;B{i+1,i}=N
" l3 A, R! o6 |5 j$ lend
6 i0 F; E' l5 D$ I" i! W2 ^A=cell2mat(B);  
这里总是显示
% q5 o" Q& \2 C" s/ E# _??? function J=jac(A,b,u0,eps)
" d. A0 E2 ?, m- ~    |
& w1 w# Y8 f2 x0 MError: Function definitions are not permitted at the prompt or in scripts.

& M3 e* k' T) g" I( M+ g2 \%定义b
b=ones(e*f,1);
for i=1:e
9 d: G' u3 j1 I0 V7 m' o, d    for j=1:f
) l1 `. R( p& \4 \7 n. H2 z& _        b(i+j)=p(i,j+1)+p(i+1,j)+p(i+2,j+1)+p(i+1,j+2)
/ k5 V7 i; f! @; f    end
6 Q) ^: E! |/ {7 K9 fend
, s# b" i0 h& j) N/ ~& ]) ]%运用Jacobi迭代法计算
* w2 |; S: [6 k8 s' B8 k. m5 pD=diag(diag(A));
( w5 P3 W( z; r. L4 g1 U  cD=inv(D);
L=tril(A,-1);
U=triu(A,1);
B=-D*(L+U);
f=D*b;
. L1 q% |/ J9 A; qJ=B*u0+f;
4 c+ y# B! H' @: y& lwhile norm(J-u0)>=eps
x0=J;
J=B*u0+f;
& Q- P/ H3 `0 s& oend
/ h. H2 I; H. E7 O6 [1 Z9 preturn

舒米牛牛

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

BenCam

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

madio

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