anchor.diff:      double precision                        :: error,x
       double precision                        :: t1,t0
       integer                                 :: i,j,k
  -    integer                                 :: quote,my_id
  +    integer                                 :: ratio,my_id
       integer,dimension(nr_threads)           :: start_it,stop_it
       character(len=20)                       :: str
       double precision,dimension(0:n-1,0:nr_threads-1)   :: padding_before,padding_after
  
positive.code:     T(0:n+1 , n+1)   = 1.0D0
      T(n+1   , 0:n+1) = 2.0D0
  
      quote = n/nr_threads
  
      do i=0,nr_threads-2
          start_it(i) = (i*quote)+1
          stop_it(i) = (i+1)*quote
      end do
      start_it(nr_threads-1) = ((nr_threads-1)*quote)+1
      stop_it(nr_threads-1) = n
  
      call omp_set_num_threads(nr_threads)
  
negative.code:     ! Uncomment the next part if you want to write the whole solution
      ! to a file. Useful for plotting.
  
      ! open(unit=7,action='write',file='result.dat',status='unknown')
      ! write(unit=str,fmt='(a,i6,a)') '(',N,'F10.6)'
      ! do i=0,n+1
      !   write (unit=7,fmt=str) T(i,0:n+1)
      ! end do
      ! close(unit=7)
  
  end program laplsolv
  