Log base 10 -the mantissa

By Reinaldo Baretti Machín
www.oocities.org/serienumerica
e-mail: reibaretti2004@yahoo.com


c calculates log base 10   mantissa
      real mantisa
      dimension pt(15),mantisa(15)
      pwr=1./2.
      do 10 i=1,12
      if(i.eq.1)pt(i)=sqrt(10.)
      if(i.ge.2)pt(i)=sqrt(pt(i-1))
      mantisa(i)=pwr**i
c      print*,'i,mantisa(i),pt(i)=',i,mantisa(i),pt(i)
10    continue
      print *,'enter positive number    1.