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

justlikethisis

function L2=lunwen(N,M)
1 k* {7 n* |  l6 Klx=0;rx=1;lt=0;rt=1;
, w/ D6 z/ }# A5 G0 A6 q2 W: q1 fh=(rx-lx)/N;
; i' _% j/ ~  c# N. ^th=(rt-lt)/M;
x=0:h:1;
t=0:th:1;
1 U' h4 J3 d6 T% gA=zeros(N+1,N+1);
B=zeros(N+1,N+1);
- S5 T9 N+ b3 Y, u" vD=zeros(N+1,N+1);
  B, I! [3 \8 D" A/ Q8 fE=zeros(2*N+2,2*N+2);
F=zeros(N+1,1);
* d  C, p$ f1 ju3=zeros(N+1,1);
0 f- U# s6 N# Y8 V; S) V( EP=zeros(2*N+2,1);
X=zeros(2*N+2,1);
* K& e+ I, d" A' fX1=zeros(2*N+2,1);
1 P6 B! e; w9 m3 jX2=zeros(2*N+2,M+1);
for i=0:N
    u3(i+1,1)=sin(pi*i*h);
  J/ j5 |4 k% U' K% o) F/ h# Mend
X(N+2:2*N+2,1)=u3(1:N+1,1);
, e7 _, n: t4 g4 x8 g4 `for i=1:N
   e=[x(i),x(i+1)];
# e; x2 o( L0 f/ _1 R% D    A([i,i+1],[i,i+1])=A([i,i+1],[i,i+1])+ganddu1(e);
% W6 ^4 h/ ]) {% K    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);
end
2 C8 O; e* ^5 w9 y0 Z, i% E(1:N+1,1:N+1)=A;
% E(1:N+1,N+2:2*N+2)=-B;
. w/ L6 E. T* L" f, P: d5 p: B% 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];
0 B- n. A6 y9 O7 v4 ^: W, n: Dfor j=1:M
    F=zeros(N+1,1); %%记得归零
3 ?! q+ \! b4 {  R/ ^; r0 j    y=exp(-t(j+1));
    for i=1:N
# x+ |$ @, o+ G1 T$ |3 m) j  z7 R! [        e=[x(i),x(i+1)];
; M+ Z# ~# Q3 q* |, {        F([i,i+1],1)=F([i,i+1],1)+y*Hezai3(e);
4 X: E+ {+ F# B6 r! B7 w    end
: X5 L7 f4 P+ i" y    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
3 C; l: R8 h4 p/ m: Y+ Jend
0 d# I* j* i, F! I" P. OX1=[A\B*X(N+2:2*N+2,1);X(N+2:2*N+2,1)];
2 J  N5 g7 s. v: I1 d$ `: nX2=[X1;X];
错误使用 vertcat
& E! M& T$ [2 n3 z, B8 m5 h( a9 R要串联的数组的维度不一致。

  E7 r3 ?! N+ q9 d. B5 N出错 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)];
. {6 T1 n0 d- t% v检查半天也不知道哪错了
1 E& B" y3 e  ~/ _# g4 W2 i
! R6 C# [* w8 }

qq666888qqw

你的程序没传完吧

somethingabc

程序好像没传完