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

justlikethisis

function L2=lunwen(N,M)
- a& x7 A" n  Q# K8 y6 llx=0;rx=1;lt=0;rt=1;
h=(rx-lx)/N;
. l2 {6 Q0 ]! M7 h, T, q9 N0 bth=(rt-lt)/M;
x=0:h:1;
. d% F" j0 m6 rt=0:th:1;
3 j8 D& L4 |4 O. o# e  ~/ ?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);
8 ]6 C- h. z% Q" AF=zeros(N+1,1);
  M1 i3 {. q: s; y8 zu3=zeros(N+1,1);
P=zeros(2*N+2,1);
# f3 w8 y' ?/ q7 n  V* R( S4 MX=zeros(2*N+2,1);
X1=zeros(2*N+2,1);
. m7 w7 V) {2 N( z$ u3 h, h2 JX2=zeros(2*N+2,M+1);
) M: c1 m* P1 H5 ffor i=0:N
7 o6 v- G, P6 G  y    u3(i+1,1)=sin(pi*i*h);
* u  A' D% n5 z0 Q7 r" oend
7 h6 ?! {& e9 w! x  BX(N+2:2*N+2,1)=u3(1:N+1,1);
for i=1:N
5 v3 Y- T# T: E. W0 [! U   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);
end
+ Z4 m( B3 [+ l( Z% 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;
E=[A,-B;-th.*B,B-th.*B+th.*D];
! r" D9 }' X" C3 y: {for j=1:M
    F=zeros(N+1,1); %%记得归零
. W! j/ E( F$ X  W. b: b6 B5 X. q    y=exp(-t(j+1));
    for i=1:N
3 a8 ^' [2 p+ s        e=[x(i),x(i+1)];
( s: S* d! `4 e/ r9 C9 Z# V4 f        F([i,i+1],1)=F([i,i+1],1)+y*Hezai3(e);
& T( R. Z1 X  T( [$ |- }# R8 v% R    end
' G& R) t# c# Z/ ^/ m: p    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
end
' A. p" m* e( [; A6 vX1=[A\B*X(N+2:2*N+2,1);X(N+2:2*N+2,1)];
X2=[X1;X];
4 [* S! a- ^6 \& {错误使用 vertcat
7 G2 t: o  t% ^. F, G要串联的数组的维度不一致。

出错 lunwen (line 39)
6 ~% u. D( S4 S) Y# ?- {- v: l9 f    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
检查半天也不知道哪错了
9 C8 R. X  f9 u1 |; U2 l  y1 n

" S, w* n  _+ P( u

qq666888qqw

你的程序没传完吧

somethingabc

程序好像没传完