aes - Converting a multivariate equation into a Binary Matrix in Python or Sagemath -

i have multivariate quadratic sum

x0 + x1*y1 + x0*y0 + y1

the variables in case are

{x0,x1,y0,y1,x0*y0,x1*y1,x0*y1,x1*y0}

so giving previous quadratic sum input , how can take matrix [1,0,0,1,1,1,0,0] according previous basis.

i'm working in sagemath.