possion-equation-finite-element.zip_naturallyvop_possion_possion

[point,pnum,element,enum]=triangle([0,pi],[0,1],33,33); A=zeros(pnum,pnum);b=zeros(1,pnum); for i=1:enum [gx,s]=squarepoint(point,element,i); [ji,djix,djiy]=hanshu(point,element,i,gx,s); for k=1:3 ii=element(i,k); b(ii)=b(ii)+s*fn(gx)*ji(k); for m=1:3 jj=element(i,m); A(ii,jj)=A(ii,jj)+s*(djix(k)*djix(m)+djiy(k)*djiy(m)); end end end for i=1:33 b(i)=0; A(i,:)=0; A(i,i)=1; end for i=33*32+1:33*33 b(i)=0; A(i,:)=0; A(i,i)=1; end u=b/A; k=1; for i=1:33 for j=1:33 z(i,j)=u(k); k=k+1; end end mesh([0:pi/32:pi],[0:1/32:1],z)