program initialisation
! This program creates a file inputfile.dat
! 

use randvectors
integer :: iprint,n,choix,choixsol
real*8 ::choixtol
print*,'Would you print results ? : 1-No / 2-All / 3-Time and final residual norm'
read*,iprint
print*,'Enter matrix size : '
read*,n
print*,'Enter matrix choice : (more details in README.txt) '
print*,'1- A(i,j)=2j-1 / n-i+j for j=1,k   '
print*,'          2i-1 / n-i+j for j=i+1,n '
print*,'2- A(i,j)=i for j=1,i '
print*,'          j for j=i+1,n '
print*,'3-helsing.f90'
print*,'4-rbf.f90'
print*,'5- Matrix File (copy in matrixfile.dat) '
read*,choix
if (choix==3) then
    call rand0(n)
    call vecrand(n)
    call helsing(n)
end if
if (choix==4) then
    call rand0(n)
    call vecrand(n)
    call rbf(n)
end if


print*,'Enter solution choice : 1- [1,...,1]^T  / 2- [n,...,1]^T '
print*,'3-rand vector'
read*,choixsol

print*,'Enter tol :'
read*,choixtol

open(unit=10,file='inputfile.dat',status='unknown')
write(10,*) iprint,n,choix,choixsol,choixtol
close(10)
end  program initialisation