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

justlikethisis

function L2=lunwen(N,M)
+ x2 H3 [; @! k0 U' olx=0;rx=1;lt=0;rt=1;
3 ^( q: d; k' D0 `7 J$ E) uh=(rx-lx)/N;
th=(rt-lt)/M;
x=0:h:1;
t=0:th:1;
) ^3 g: T' [* z' w( i- S6 T1 X" ^A=zeros(N+1,N+1);
B=zeros(N+1,N+1);
D=zeros(N+1,N+1);
/ W* n4 i: Y8 @* OE=zeros(2*N+2,2*N+2);
& j0 U: }$ ~5 C- J7 {2 k% BF=zeros(N+1,1);
8 Z- f& r9 A2 j; b  o9 ^( W  Au3=zeros(N+1,1);
# @( X7 s* q( e' {! \% v8 p/ ~P=zeros(2*N+2,1);
X=zeros(2*N+2,1);
7 }  A4 {9 I1 T. h3 IX1=zeros(2*N+2,1);
X2=zeros(2*N+2,M+1);
, S; _) a0 M" gfor i=0:N
8 u/ `2 Y, ^/ b8 u! Y- Q. f    u3(i+1,1)=sin(pi*i*h);
end
" ^; ?. x5 S% z7 w& M: ZX(N+2:2*N+2,1)=u3(1:N+1,1);
for i=1:N
( L2 i* g# h$ b, p: h+ B, w. T   e=[x(i),x(i+1)];
8 V, T4 B  Z2 Q. j- }/ p- I) T" |    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);
1 R' E$ B' K& C2 i4 t, x    D([i,i+1],[i,i+1])=D([i,i+1],[i,i+1])+gangdu33(e);
1 D+ X! s! G8 tend
% E(1:N+1,1:N+1)=A;
( R; `/ g4 G- K* i0 t% E(1:N+1,N+2:2*N+2)=-B;
2 @! @6 L0 p  P3 ?8 d# x# p# F9 P% 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;
3 q0 G- }) ^% {! @2 ~2 x( Z/ i/ TE=[A,-B;-th.*B,B-th.*B+th.*D];
6 j/ z1 M5 F+ S0 y7 Rfor j=1:M
    F=zeros(N+1,1); %%记得归零
6 t) `6 ^$ P% w% N5 Y( c( q    y=exp(-t(j+1));
    for i=1:N
- ~+ r! X5 W3 R, W) `: J" o        e=[x(i),x(i+1)];
7 `/ F( L5 T) ~1 j        F([i,i+1],1)=F([i,i+1],1)+y*Hezai3(e);
    end
4 Y  U4 m+ M3 `2 J3 h( o    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
3 L' q8 b8 Y! }end
$ i8 B+ u8 G' z2 B. N/ c7 q3 X, xX1=[A\B*X(N+2:2*N+2,1);X(N+2:2*N+2,1)];
7 D6 S% e' F% a( AX2=[X1;X];
4 t2 k6 ?. ^; i7 A3 q6 [- ?错误使用 vertcat
要串联的数组的维度不一致。
" V' u1 E; E% Y7 r
4 S( s5 p7 X+ _) U出错 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)];
检查半天也不知道哪错了
) j: s) [3 f! ?( u. L9 f4 d

' K7 e2 a* j( L; d

qq666888qqw

你的程序没传完吧

somethingabc

程序好像没传完