New double Wronskian exact solutions for a generalized (2+1)-dimensional nonlinear system with variable coefficients

New generalized (2+1)-dimensional Boussinesq system with variable coefficients has been introduced. A double Wronskian solutions has been formulated to the new system under certain constraints on the variable coefficients. Hirota differential operator and its properties have been employed to transform the system into a bilinear form. Moreover, in order to better understand the dynamic behavior, the characteristics of the double Wronskian solutions are discussed through some diverting graphics under different parameters choices.