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

justlikethisis

function L2=lunwen(N,M)
lx=0;rx=1;lt=0;rt=1;
# ]7 R, S+ X; p- O& A8 E% Jh=(rx-lx)/N;
% E* G/ Y% d/ i' xth=(rt-lt)/M;
+ K  Z  z; D! _8 d+ A* N: V8 {2 wx=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);
% o! p6 [8 J+ S: n& [$ R; r4 MF=zeros(N+1,1);
" M* c3 Z" T6 v) V! su3=zeros(N+1,1);
5 Z  w: K* Z0 F% T; Z0 g0 ]8 W0 JP=zeros(2*N+2,1);
X=zeros(2*N+2,1);
X1=zeros(2*N+2,1);
1 H( \9 n8 h2 S+ a6 c' ?0 iX2=zeros(2*N+2,M+1);
for i=0:N
+ _, \: K; b/ X4 j    u3(i+1,1)=sin(pi*i*h);
3 Q% y% ?- A7 Mend
X(N+2:2*N+2,1)=u3(1:N+1,1);
for i=1:N
   e=[x(i),x(i+1)];
; ^4 `5 g+ [& U0 ?* w    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);
/ ]5 o2 j; V" I9 k/ {# c9 kend
% E(1:N+1,1:N+1)=A;
% E(1:N+1,N+2:2*N+2)=-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;
, g  t2 Z5 `3 w0 ?8 u0 K* y5 h! mE=[A,-B;-th.*B,B-th.*B+th.*D];
9 `: ^* T- t& }2 ?for j=1:M
* ~! G5 n. U9 ^6 w' q2 g, o8 V% K    F=zeros(N+1,1); %%记得归零
    y=exp(-t(j+1));
    for i=1:N
        e=[x(i),x(i+1)];
        F([i,i+1],1)=F([i,i+1],1)+y*Hezai3(e);
2 b) e/ w; y0 Y: P0 M    end
- |1 U) s/ b% s8 e    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
end
5 L, C0 b. ^1 z4 BX1=[A\B*X(N+2:2*N+2,1);X(N+2:2*N+2,1)];
X2=[X1;X];
" ?1 |; _& h  w0 ~2 _错误使用 vertcat
要串联的数组的维度不一致。

/ r% w5 w% O- k1 z$ u' G$ P4 a出错 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)];
检查半天也不知道哪错了
8 s% V5 V# T5 S; Z6 w
) C6 o6 w% A( V  c; N: n- P

qq666888qqw

你的程序没传完吧

somethingabc

程序好像没传完