I changed this:

sigmoid + sigmoid + sigmoid + meanSquaredError

to:

relu + relu + sigmoid + binaryCrossentropy

And seems that the prediction is quite better:

Tensor

[[0.9998419, 4e-7 ],

[0.9998978, 8e-7 ],

[0 , 0.9999998],

[0 , 0.9999999]]