{equation}\left\{\begin {aligned} &f (x+ y, z)+ f (xy, z)= 2f (x, z)+ f (y, z)+ f (-y,-z)\nonumber\\&f
(x, y+ z)+ f (x, yz)= 2f (x, y)+ f (x, z)+ f (-x,-z)\nonumber\end {aligned}\right.\end {equation} for
all $ x, y, z\in X $. We will also investigate the Hyers-Ulam stability of the system of bi-Drygas
functional equations.