Quasi-classical modeling of a spin-boson system H=sum(pb^2+qb^2)+c*qb*(p1^2+x1^2+p2^2+x2^2)+delta*(p1^2+x1^2-p2^2-x1^2)