1 |       program testprof
   2 | 
   3 |       implicit none
   4 |       
   5 |       integer n,i
   6 |       parameter (n= 10)
   7 |       double precision r(0:n+1)
   8 |       double precision rayy, black.f" HREF="../src/profblack.f.html">profblack, black
   9 | 
  10 | 
  11 | 
  12 |       do i=0, n+1
  13 |          rayy =((dble(i)/dble((n+1))*(50-25)))+25
  14 |          black= profblack(25.d+0,50.d+0, rayy)
  15 |          print *, rayy , black
  16 |       enddo
  17 | 
  18 |       end
  19 | 
  20 | 
  21 | 
  22 | 
  23 |