part2.m

利用矩阵分析的方法对线性方程组进行求解分析

function[y]=part2(a,t);
w=10;
if((mod(a,2))==1)
    y=((-1).^((a+1)/2))*((w).^a)*sin(w*t);
else
    y=((-1).^((a)/2))*((w).^a)*cos(w*t);
end
%A中的第二部分，有要改的w