/* A is matrix, FILENAME is file name */

def asir2bin_charpoly_minipoly(A,FILENAME){

	N=size(A)[0]$

	Z=0$
	for(I=0;I<N;I++){
		for(J=0;J<N;J++){
			if(A[I][J]!=0)
				Z++$
		}
	}

	output(FILENAME)$

	print("Z",2);
	print(Z)$

  	print("N",2)$
	print(N)$

	for(I=0;I<N;I++){
		for(J=0;J<N;J++){

			if(A[I][J]==0){
				continue$
			}

			print("I",2)$
			print(I,2)$
			print(",",2)$
			print(J,2)$
			print(",",2)$
			print(A[I][J])$
		}
	}

	output()$
}

end$