File: intrinsics_221.f90

package info (click to toggle)
lfortran 0.45.0-1
  • links: PTS, VCS
  • area: main
  • in suites: sid, trixie
  • size: 46,332 kB
  • sloc: cpp: 137,068; f90: 51,260; python: 6,444; ansic: 4,277; yacc: 2,285; fortran: 806; sh: 524; makefile: 30; javascript: 15
file content (16 lines) | stat: -rw-r--r-- 490 bytes parent folder | download 
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
 
program intrinsics_221

    implicit none
    integer, parameter :: N = 10
    integer, parameter :: sp = kind(1.0)
    integer, parameter :: dp = kind(1.d0)
    real(dp), parameter :: pi = 2*asin(1._dp)
    real(dp), parameter :: a = 0, b = pi
    real(dp), parameter :: dx = (b-a)/N
    integer :: i
    real(dp), parameter :: X(N) = [(sin(a+(b-a)*i/N), i = 1, N)]
    real(dp), parameter :: S = sum(X)*dx
    print*, S
    if (abs(S - 1.983523537509) > 1e-7) error stop
    
end program