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

justlikethisis

function L2=lunwen(N,M)
3 g# F7 p* Q" j, p' U& e, dlx=0;rx=1;lt=0;rt=1;
( z8 x! m* U( A6 L3 W$ }h=(rx-lx)/N;
th=(rt-lt)/M;
' z2 \" |, x. ]# Sx=0:h:1;
t=0:th:1;
# j1 e+ A! \0 ?% d* wA=zeros(N+1,N+1);
B=zeros(N+1,N+1);
8 @" l' H* C' M$ cD=zeros(N+1,N+1);
* v; F3 X8 s; aE=zeros(2*N+2,2*N+2);
F=zeros(N+1,1);
! {" s0 q& R4 O- W% M. Lu3=zeros(N+1,1);
$ a# g2 `8 [  ]7 g5 t: r, ?- t- uP=zeros(2*N+2,1);
. s( V" l& {& X# v9 H. n/ O4 _X=zeros(2*N+2,1);
X1=zeros(2*N+2,1);
X2=zeros(2*N+2,M+1);
for i=0:N
    u3(i+1,1)=sin(pi*i*h);
7 ~( n! g9 K( Q  M! rend
2 I8 s) \, G  ?4 q, D+ CX(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);
( Y) [" Q* [. Q( j    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);
( f& d8 ^" k/ ^9 Yend
" a0 M; q6 p" U6 @/ x3 n" J% E(1:N+1,1:N+1)=A;
% E(1:N+1,N+2:2*N+2)=-B;
$ N/ i* t# e# ^6 B. p; Q" h% 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;
. N- ~- z% b. x6 I; @E=[A,-B;-th.*B,B-th.*B+th.*D];
for j=1:M
    F=zeros(N+1,1); %%记得归零
% P! A. P; \: q) M7 {0 m    y=exp(-t(j+1));
    for i=1:N
9 E2 Y9 e' E) D/ d) {- S# r        e=[x(i),x(i+1)];
% z1 V* a% k6 r1 l        F([i,i+1],1)=F([i,i+1],1)+y*Hezai3(e);
" T7 I& \/ i; z    end
* W  I9 I7 g3 i0 R    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
5 K8 y+ F( J8 Q) V& S7 }+ @" |end
X1=[A\B*X(N+2:2*N+2,1);X(N+2:2*N+2,1)];
X2=[X1;X];
& X; \1 k) E0 O2 k  n错误使用 vertcat
要串联的数组的维度不一致。

出错 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)];
检查半天也不知道哪错了

" o& B+ k* ]9 [4 v

qq666888qqw

你的程序没传完吧

somethingabc

程序好像没传完