program maillage

implicit none

integer :: N, k
real :: L, gama, pi
real, dimension(100) :: x_reg, x_irreg

! ouverture du fichier d'entree et recuperation des donnees du fichier d'entree
  open(10,file='donnees.dat')
  read(10,*) N
  read(10,*) L
  read(10,*) gama
  close(10)

! calcul des coordonnees du maillage regulier
  do k=1,N
    x_reg(k) = L * (real(k)-1) / (real(N)-1.)
  end do

! calcul des coordonnees du maillage irregulier
  pi = acos(-1.)
  do k=1,N
    x_irreg(k) = x_reg(k) + gama * L / (3.*pi) * sin(2.*pi*x_reg(k)/L)
  end do

! creation du fichier de resultats et ecriture dans ce fichier
  open(20,file='resultats.dat')
  do k=1,N
    write(20,*) k, x_irreg(k)
  end do
  close(20)

end program maillage