clear all

% digits(4)
x    = 0.54617;
y    = 0.54601;
fl_x = vpa(x,4);
fl_y = vpa(y,4);
z1   = fl_x - fl_y;
z    = double(z1)