part3_1_2.m

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

%A的第三函数，sa(x)的几次求导，与w,B都无关

function[y]=part3_1_2(b,a);
k=-1;
sym x;
y=sinint(x);
while(a>k)
    y=diff(y,x);
    k=k+1;
end
x=b;
y;

  %求sa(x)的几次求导