Hello all,(adsbygoogle = window.adsbygoogle || []).push({});

I am working on a Fortran90 code that will give me an approximation of a few mathematical constants. My code will compile and run, however, my outputs are given as 0.000000 for all constants that I am calculating in the loops. Any help would be greatly appreciated. My code is as follows:

PROGRAM HomeWorkOne_ProbOne

IMPLICIT NONE

! Computes an estimation of the exponential function, pi, and the Euler Constant

REAL :: e, EC, pi, n_factorial, x_factorial, gamma

INTEGER :: i, n, x

x = n + 3.0/2.0

n_factorial = 1.0

DO WHILE (n < 10)

DO i = 1, n

n_factorial = n_factorial * i

END DO

e = 1/n_factorial

END DO

DO WHILE (n < 10)

EC = (1.0/n) - (LOG(REAL(n)+1.0)-LOG(REAL(n)))

END DO

DO WHILE (x < 10)

DO i = 1, x

x_factorial = x * i

gamma = (2.0 * x_factorial)/(2.0 * n + 3.0)

END DO

DO i = 1, n

n_factorial = n_factorial * i

END DO

pi = (n_factorial)/((2**(n-(1.0/4.0)*Cos(3.14*(2.0*n+1.0))))+(3.0/4.0)*(3.14**((1.0/4.0)*Cos(3.14*(2*n+1))-(1.0/4.0)))*gamma)

END DO

PRINT *, "Approximate values of:"

PRINT *, "e =", e

PRINT *, "gamma =", EC

PRINT *, "Pi =", pi

PRINT *, "Actual values of :"

PRINT *, "e =", 2.71828182

PRINT *, "gamma =", 0.57721566

PRINT *, "Pi =", 3.14159265

END PROGRAM HomeWorkOne_ProbOne

Thanks for any help!!

**Physics Forums | Science Articles, Homework Help, Discussion**

Join Physics Forums Today!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

The friendliest, high quality science and math community on the planet! Everyone who loves science is here!

# Fortran Help - Do loops not calculating.

**Physics Forums | Science Articles, Homework Help, Discussion**