home *** CD-ROM | disk | FTP | other *** search
/ Power-Programmierung / CD1.mdf / fortran / compiler / bcf / mathtest.f < prev    next >
Text File  |  1990-04-25  |  4KB  |  119 lines
  1.       program mathte
  2.       CHARACTER c
  3.       CALL r
  4.       print '(40x,"...Return")'
  5.       read '(a)',c
  6.       CALL d
  7.       END
  8.  
  9.  
  10.       SUBROUTINE r
  11.       logical l1,l2,l3,l4,l5,l6,l7,l8
  12.       print *,'REAL'
  13.         a=0.1
  14.         b=0.
  15.         c=-a
  16.         apb=a+b
  17.         cm1=c-1.
  18.         l1=a.gt.b
  19.         l2=b.gt.c
  20.         l3=a.gt.c
  21.         l7=c.gt.a
  22.         l4=apb.lt.c
  23.         l5=c.lt.cm1
  24.         l6=a.le.b
  25.         l8=c.lt.a
  26.         print *,' t t t t f f f f'
  27.         print *,l1,l2,l3,l8,l4,l5,l6,l7
  28.       a=0.5
  29.       b=10.
  30.       print *,a,'+',b,'=',a+b
  31.      1 ,a,'-',b,'=',a-b
  32.      1 ,a,'*',b,'=',a*b
  33.      1 ,a,'/',b,'=',a/b
  34.       print *,' abs(-1.)   =',abs(-1.)
  35.      1     ,'   abs(1.)    =',abs(1.)
  36.       print *,' nint(-1.2) =',nint(-1.2)
  37.      1     ,'   nint(1.2)  =',nint(1.2)
  38.       print *,'sign(-1.,1.)=',sign(-1.,1.)
  39.      1     ,'  sign(1.,-1.)=',sign(1.,-1.)
  40.       print *,' dim(-1.,1.)=',dim(-1.,1.)
  41.      1     ,'   dim(1.,-1.)=',dim(1.,-1.)
  42.       b=100.
  43.       d=asin(sin(a))
  44.       print '(a,f4.2,a,g16.7)','    asin(sin(',a,'))=',d
  45.       d=acos(cos(a))
  46.       print '(a,f4.2,a,g16.7)','    acos(cos(',a,'))=',d
  47.       d=atan(tan(a))
  48.       print '(a,f4.2,a,g16.7)','    atan(tan(',a,'))=',d
  49.       d=tan(atan2(a,1.))
  50.       print '(a,f4.2,a,g16.7)','tan(atan2(',a,',1.))=',d
  51.       d=log(exp(a))
  52.       print '(a,f4.2,a,g16.7)','     log(exp(',a,'))=',d
  53.       d=log10(b)
  54.       print '(a,f4.0,a,g16.7)','        log10(',b,')=',d
  55.       d=b**2.
  56.       d=sqrt(d)
  57.       print '(a,f4.0,a,g16.7)','     sqrt(',b,'**2.)=',d
  58.       d=cosh(a)**2.-sinh(a)**2.
  59.       print '(2(a,f3.1),a,g16.7)','cosh(',a,')**2-sinh(',a,')**2=',d
  60.       d=tanh(a)
  61.       print '(a,f3.1,a,g16.7)','         tanh(',a,')=',d
  62.       end
  63.  
  64.       SUBROUTINE d
  65.       logical l1,l2,l3,l4,l5,l6,l7,l8
  66.       DOUBLE PRECISION a,b,d
  67.       print *,'DOUBLE PRECISION '
  68.         a=0.1d0
  69.         b=0.d0
  70.         c=-a
  71.         apb=a+b
  72.         cm1=c-1.d0
  73.         l1=a.gt.b
  74.         l2=b.gt.c
  75.         l3=a.gt.c
  76.         l7=c.gt.a
  77.         l4=apb.lt.c
  78.         l5=c.lt.cm1
  79.         l6=a.le.b
  80.         l8=c.lt.a
  81.         print *,' t t t t f f f f'
  82.         print *,l1,l2,l3,l8,l4,l5,l6,l7
  83.       a=0.5d0
  84.       b=10d0
  85.       print *,a,'+',b,'=',a+b
  86.      1 ,a,'-',b,'=',a-b
  87.      1 ,a,'*',b,'=',a*b
  88.      1 ,a,'/',b,'=',a/b
  89.       print *,'      abs(-1.d0)=',abs(-1.d0)
  90.      1     ,'         abs(1.d0)=',abs(1.d0)
  91.       print *,'   nint(-1.2d0) =',nint(-1.2d0)
  92.      1     ,'     nint(1.2d0)  =',nint(1.2d0)
  93.       print *,'sign(-1.d0,1.d0)=',sign(-1.d0,1.d0)
  94.      1     ,'  sign(1.d0,-1.d0)=',sign(1.d0,-1.d0)
  95.       print *,' dim(-1.d0,1.d0)=',dim(-1.d0,1.d0)
  96.      1     ,'   dim(1.d0,-1.d0)=',dim(1.d0,-1.d0)
  97.       b=100d0
  98.       d=asin(sin(a))
  99.       print '(a,f4.2,a,g25.14)','    asin(sin(',a,'))=',d
  100.       d=acos(cos(a))
  101.       print '(a,f4.2,a,g25.14)','    acos(cos(',a,'))=',d
  102.       d=atan(tan(a))
  103.       print '(a,f4.2,a,g25.14)','    atan(tan(',a,'))=',d
  104.       d=tan(atan2(a,1.D0))
  105.       print '(a,f4.2,a,g25.14)','tan(atan2(',a,',1.))=',d
  106.       d=log(exp(a))
  107.       print '(a,f4.2,a,g25.14)','     log(exp(',a,'))=',d
  108.       d=log10(b)
  109.       print '(a,f4.0,a,g25.14)','        log10(',b,')=',d
  110.       d=b**2.D0
  111.       d=sqrt(d)
  112.       print '(a,f4.0,a,g25.14)','     sqrt(',b,'**2.)=',d
  113.       d=cosh(a)**2.-sinh(a)**2.
  114.       print '(a,2(f3.1,a),g25.14)','cosh(',a,')**2-sinh(',a,')**2=',d
  115.       d=tanh(a)
  116.       print '(a,f3.1,a,g25.14)','         tanh(',a,')=',d
  117.       end
  118.  
  119.