3 (a + b) - Sqrt[(3 a + b) (a + 3b)] = 1, a = 1/2*3/2, 2b