subroutine le_disk
! Variables
integer :: N
integer :: flag
double precision, dimension (1:N) :: d
double precision, dimension (1:N) :: u
integer :: i
double precision :: c1
double precision :: c2
double precision :: c3
end subroutine le_disk
integer :: N
integer :: flag
double precision, dimension (1:N) :: d
double precision, dimension (1:N) :: u
integer :: i
double precision :: c1
double precision :: c2
double precision :: c3