boussinesq equation final version: solved through Casulli's method, by adding 1 to null diagonal entry of matrix T