function comparison theta = pi/6; x = tan(theta); y = sin(theta)/cos(theta); if x == y fprintf('x and y are equal\n'); else fprintf('x and y are not equal: x-y = %e\n', x-y); end