subroutine dupa(x)
       double precision x
       read (*,*) x
       return
       end