一维抛物有限元编程，检查半天也不知道哪错了，大家可以帮忙看看吗，谢谢了

justlikethisis

function L2=lunwen(N,M)
lx=0;rx=1;lt=0;rt=1;
h=(rx-lx)/N;
th=(rt-lt)/M;
2 Y* B! a. A5 P  ]' H- p4 q; a/ }x=0:h:1;
t=0:th:1;
* d# V6 v& L3 }. g" H2 fA=zeros(N+1,N+1);
8 L- c: F3 s/ O  ?, P( \B=zeros(N+1,N+1);
D=zeros(N+1,N+1);
# v4 K+ g4 z! ]% |1 W. \/ eE=zeros(2*N+2,2*N+2);
. k3 x2 r4 ~2 ^. s9 F" _F=zeros(N+1,1);
& W) ~" I8 V/ Yu3=zeros(N+1,1);
P=zeros(2*N+2,1);
X=zeros(2*N+2,1);
9 B2 E5 E0 ?3 \! {6 `X1=zeros(2*N+2,1);
6 b2 V9 h2 {3 T- z' C8 Q7 BX2=zeros(2*N+2,M+1);
/ M- U1 C; b  G9 a5 X7 Tfor i=0:N
    u3(i+1,1)=sin(pi*i*h);
  \( n) b5 @. F" d, mend
X(N+2:2*N+2,1)=u3(1:N+1,1);
( G8 D; o5 R7 [$ o- tfor i=1:N
$ R+ R7 s* L( y   e=[x(i),x(i+1)];
    A([i,i+1],[i,i+1])=A([i,i+1],[i,i+1])+ganddu1(e);
    B([i,i+1],[i,i+1])=B([i,i+1],[i,i+1])+gangdu22(e);
3 |, @* w9 X' j( l8 P2 V, X. z/ l    D([i,i+1],[i,i+1])=D([i,i+1],[i,i+1])+gangdu33(e);
end
% E(1:N+1,1:N+1)=A;
( a  @( A' {8 t  b# {4 h' l5 l% E(1:N+1,N+2:2*N+2)=-B;
# k! p$ O6 f( n; d1 C% E(N+2:2*N+2,1:N+1)=-th*B;
4 l/ N. ?0 s2 _. _1 y" \% E(N+2:2*N+2,N+2:2*N+2)=B-th*B+th*D;
E=[A,-B;-th.*B,B-th.*B+th.*D];
for j=1:M
2 Y) e- ~/ Q$ b7 h    F=zeros(N+1,1); %%记得归零
    y=exp(-t(j+1));
    for i=1:N
        e=[x(i),x(i+1)];
, ?0 L9 r7 c7 y8 s) S# L& H        F([i,i+1],1)=F([i,i+1],1)+y*Hezai3(e);
. @% ]" E' y% k$ i    end
5 d1 ]3 Y, }, _  Y    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
end
5 V  ~7 _& B5 W4 |' zX1=[A\B*X(N+2:2*N+2,1);X(N+2:2*N+2,1)];
X2=[X1;X];
错误使用 vertcat
) @5 D. b, E* o% M( n要串联的数组的维度不一致。
% l+ P/ ?, Y2 R0 g
出错 lunwen (line 39)
    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
检查半天也不知道哪错了

' ^* I' A- i: j! r$ z2 V( k' R
1 H( Y# a8 S9 h4 a

qq666888qqw

你的程序没传完吧

somethingabc

程序好像没传完