load("./asir2bin_det_simple.rr")$

def seki_res(C,A,B){

	M=deg(A,C)$
	N=deg(B,C)$

	if(M>=N){
		L=M$
		S=N$
	}
	else{
		L=N$
		S=M$
		T=A$
		A=B$
		B=T$
	}

	MAT=newmat(L,L)$

	for(I=0;I<(L-S);I++)
		for(J=0;J<=S;J++)
			MAT[I][J+I]=coef(B,J,C)$

	B=C^(L-S)*B$

	F=srem(A,C^L)$
	G=srem(B,C^L)$
	H=coef(A,L,C)*G-coef(B,L,C)*F$
	for(J=0;J<L;J++){
		MAT[2*L-S-L][J]=coef(H,J,C)$
	}

	for(I=L-1;I>L-S;I--){
		H=C*srem(H,C^(L-1))+coef(A,I,C)*G-coef(B,I,C)*F$
		for(J=0;J<L;J++){
			MAT[2*L-S-I][J]=coef(H,J,C)$
		}
	}

	return MAT$
}

P=a0+a1*x+a2*x^2+a3*x^3+a4*x^4+a5*x^5+a6*x^6+a7*x^7+x^8+x^9$
MAT=seki_res(x,P,diff(P,x))$
V=[a7,a6,a5,a4,a3,a2,a1,a0]$
W1=[2,3,4,5,6,7,8,9]$
W2=[7,6,5,4,3,2,1,9]$
FILENAME="./9test3.dat"$
asir2bin_det_simple(MAT,V,[W1,W2],FILENAME)$

FILENAME="./9test4.dat"$
asir2bin_det_simple(MAT,V,[0],FILENAME)$

end$