  Analytical expressions for spectrum, eigenfunctions and dipole matrix
elements of a square double quantum well (DQW) are presented for a general case
when the potential in different regions of the DQW has different heights and
effective masses are different. This was achieved by Gr\"obner basis algorithm
which allows to disentangle the resulting coupled polynomials without
explicitly solving the transcendental eigenvalue equation.
