titand.f

Go to the documentation of this file.

      subroutine titand ( tau, gtau, n )

c*********************************************************************72
c
cc TITAND represents a temperature dependent property of titanium.
c
c  from  * a practical guide to splines *  by c. de boor    
c  these data represent a property of titanium as a function of
c  temperature. they have been used extensively as an example in spline
c  approximation with variable knots.
c
      implicit none

      integer n,   i
      double precision gtau(49),tau(49),gtitan(49)
      data gtitan /
     &  0.644d+00, 0.622d+00, 0.638d+00, 0.649d+00, 0.652d+00,
     &  0.639d+00, 0.646d+00, 0.657d+00, 0.652d+00, 0.655d+00,
     &  0.644d+00, 0.663d+00, 0.663d+00, 0.668d+00, 0.676d+00,
     &  0.676d+00, 0.686d+00, 0.679d+00, 0.678d+00, 0.683d+00,
     &  0.694d+00, 0.699d+00, 0.710d+00, 0.730d+00, 0.763d+00,
     &  0.812d+00, 0.907d+00, 1.044d+00, 1.336d+00, 1.881d+00,
     &  2.169d+00, 2.075d+00, 1.598d+00, 1.211d+00, 0.916d+00,
     &  0.746d+00, 0.672d+00, 0.627d+00, 0.615d+00, 0.607d+00,
     &  0.606d+00, 0.609d+00, 0.603d+00, 0.601d+00, 0.603d+00,
     &  0.601d+00, 0.611d+00, 0.601d+00, 0.608d+00 /
      n = 49
      do 10 i=1,n
         tau(i) = 585.0d+00 + 10.0d+00*dble(i)
   10    gtau(i) = gtitan(i)
                                        return
      end