binary.m

function a = binary（ik）
% return the coefficients of the binary expansion of i:i的二的指数表达式系数
% i = a（1）*2^（k-1） + a（2）*2^（k-2） + ... + a（k）

if i>=2^k
   error（‘i must be such that i