\(\left(1+x^2\right)\left(1+y^2\right)+4xy+2\left(x+y\right)\left(1+xy\right)\)\(\)
\(=1+y^2+x^2+x^2y^2+4xy+2\left(x+y\right)\left(1+xy\right)\)
\(=x^2+2xy+y^2+x^2y^2+2xy+1+2\left(x+y\right)\left(1+xy\right)\)
\(=\left(x+y\right)^2+\left(xy+1\right)^2+2\left(x+y\right)\left(xy+1\right)\)
\(=\left(x+y+xy+1\right)^2\)