FEA的分段近似步骤

Question

我正在尝试在0到1之间的1维中求解热方程u _t +Δu = e ^u 在域上。

要解决PDE，我想从头开始编写有限元代码，以帮助我更好地理解MATLAB和有限元方法。但是，我正在努力创建可以用作网格的分段函数。

npoint=5 %points in decomposition - points on our line
nelement=npoint-1 %number of intervals on our line

x=linspace(0,1,npoint)  %creates vertices

intervals(1:nelement,1)=1:nelement %creates e2p, part 1????

(1:nelement,2)=2:npoint %creates e2p, part 2
%e2p this descirbes the sections which will be approximated, e.g. area of
%bounday 1,2

figure=plot(x,0*x, 'b-o', 'MarkerFaceColor','r') %draw composition 
figure

如果有人可以建议我如何继续使用有限元代码或为我提供良好的指导，那么将不胜感激。非常感谢。