program main implicit none double precision X,Y,Z,W read*,X,Y,Z,W print'(f0.10)',(X-Z*Y)/dble(2) end program main