Как повысить точность вещественного деления в среде C++?
Может вы выводите на экран с точностью в 4 знака? У double явно побольше точность
Нормально все считается. У меня вышло

sizeof(double) = 8 bytes
z = 16.66666666666666785090

gcc version 9.4.0