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

justlikethisis

function L2=lunwen(N,M)
lx=0;rx=1;lt=0;rt=1;
) M# j9 _/ ?) C6 m$ Th=(rx-lx)/N;
th=(rt-lt)/M;
6 u# ]" s8 Y# l; _2 ]x=0:h:1;
t=0:th:1;
A=zeros(N+1,N+1);
B=zeros(N+1,N+1);
D=zeros(N+1,N+1);
E=zeros(2*N+2,2*N+2);
F=zeros(N+1,1);
u3=zeros(N+1,1);
P=zeros(2*N+2,1);
X=zeros(2*N+2,1);
6 y% ~3 H( R8 T# C" |X1=zeros(2*N+2,1);
X2=zeros(2*N+2,M+1);
( o! t' H* K( m7 d$ E" F! _" Qfor i=0:N
; o. Q1 J( J4 z9 O4 `    u3(i+1,1)=sin(pi*i*h);
" t$ i  f# n) `8 _8 o6 ]  Dend
' l9 H! I  t( f6 d0 ~& M. yX(N+2:2*N+2,1)=u3(1:N+1,1);
for i=1:N
   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);
    D([i,i+1],[i,i+1])=D([i,i+1],[i,i+1])+gangdu33(e);
$ P! F  k8 o  _; l1 Jend
$ d' X! |  M. B4 h: w% E(1:N+1,1:N+1)=A;
% E(1:N+1,N+2:2*N+2)=-B;
" V, z; P9 Q7 [: A+ K2 w% E(N+2:2*N+2,1:N+1)=-th*B;
% 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];
# P8 F$ S2 R1 @" J. o" Gfor j=1:M
- G4 n# W8 S; a9 {9 i, d  Y" T: M    F=zeros(N+1,1); %%记得归零
    y=exp(-t(j+1));
    for i=1:N
$ I5 ~% {. Q0 Q' S        e=[x(i),x(i+1)];
        F([i,i+1],1)=F([i,i+1],1)+y*Hezai3(e);
" o/ k' x: ?1 }    end
    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
end
1 y' u0 E2 y4 H1 L( e: i+ kX1=[A\B*X(N+2:2*N+2,1);X(N+2:2*N+2,1)];
  G% X# q6 b; z) p4 Y' |, fX2=[X1;X];
9 ^7 S. u8 r$ N5 F1 f, R错误使用 vertcat
/ \! w6 X% B' c2 c- L8 L1 G要串联的数组的维度不一致。

# E% s) q# O$ S& N( y0 F9 Y出错 lunwen (line 39)
8 ?! e) K, c# `! _( X7 ~- T    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
检查半天也不知道哪错了
: s' H% g- ]1 V/ ~! B' W9 r9 P, U

1 [1 ]. F/ t- ]' A7 c

qq666888qqw

你的程序没传完吧

somethingabc

程序好像没传完