Try to approximate l1 solution (doesn't look good at the moment).
Also tried to switch to scipy.special import binom coef, but that caused n=6 numeric solution to fail