分别用4单元和16单元有限元法求解左端固定、右端自由的等直杆(截面积为 ,弹性模量为 ,密度为 ,杆长为 )在1/4杆长处的稳态响应并与精确解进行对比(假设1/4杆长处集中载荷 激励)(要求MATLAB编程计算,附程序清单)