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

justlikethisis

function L2=lunwen(N,M)
* H. }4 W' h% elx=0;rx=1;lt=0;rt=1;
h=(rx-lx)/N;
2 n1 v9 M- h: s1 ~: k- Zth=(rt-lt)/M;
* V8 I: Z9 Z, o( q* P- fx=0:h:1;
, @0 h1 q/ {1 D8 Ot=0:th:1;
$ c, n  _5 Q% ?! c  J* K: S( tA=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);
F=zeros(N+1,1);
u3=zeros(N+1,1);
7 i7 c3 ?, K- S0 i% ?0 U3 bP=zeros(2*N+2,1);
. i; z2 D$ C+ Y% U' rX=zeros(2*N+2,1);
% w; `9 E8 m! L# P* o+ u- JX1=zeros(2*N+2,1);
X2=zeros(2*N+2,M+1);
for i=0:N
6 ^( f! h* C' U; D6 N0 |2 a0 }1 ?    u3(i+1,1)=sin(pi*i*h);
  \5 ^0 e+ Y% |. Bend
X(N+2:2*N+2,1)=u3(1:N+1,1);
# G9 i7 M5 Z4 t+ T6 ~; ]  ffor i=1:N
+ M) M* F& o9 q  \, p- T" F   e=[x(i),x(i+1)];
    A([i,i+1],[i,i+1])=A([i,i+1],[i,i+1])+ganddu1(e);
( l/ E0 B, h. a, W    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);
( L% l/ b& l4 A$ pend
& ?% s: v" q- f% 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];
$ t5 S: _( n; s; f5 Wfor j=1:M
9 g* ^) v- O. d8 c    F=zeros(N+1,1); %%记得归零
! P; V% b. Z# C' x0 T! l    y=exp(-t(j+1));
" N/ o/ [$ h# y7 q9 ?    for i=1:N
$ H3 ^5 r# Y3 T- t' ^' b        e=[x(i),x(i+1)];
2 \$ C9 _: g' T' V$ n        F([i,i+1],1)=F([i,i+1],1)+y*Hezai3(e);
    end
    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
( O& Q- `- z# Zend
; V1 M/ z% ^4 M" I4 ?X1=[A\B*X(N+2:2*N+2,1);X(N+2:2*N+2,1)];
X2=[X1;X];
错误使用 vertcat
要串联的数组的维度不一致。

/ U* u7 G: F" f" T  q" M. i出错 lunwen (line 39)
1 x& t9 P8 u7 Z# a; q& f    X(1:2.*N+2,j+1)=E\[zeros(N+1,1);-th.*F+B.*X(N+2:2.*N+2,j)];
检查半天也不知道哪错了

somethingabc

程序好像没传完

qq666888qqw

你的程序没传完吧