part2.m

来自「利用矩阵分析的方法对线性方程组进行求解分析」· M 代码 · 共 10 行

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