%求取導數公式的系數
m=input('m=');					%鍵盤輸入m
k=input('k=');					%鍵盤輸入k，k-1次多項式
n=2*m+1;
x=zeros(n,k);
for I=-m:m							%構造X矩陣行循環
    for J=0:k-1						%構造X矩陣列循環
	x(I+m+1,J+1)=I^J;		
    end
end
d=inv(x'*x)*x'