integrate [1+b^2 + 4ab