program senos
!Este programa compara o cálculo do seno de um valor
!pelo desenvolvimento em série com a função intrinseca do FORTRAN
real x, seno, sen1, sen2, sen3
   write (*,*) "entre com o valor de x"
   read (*,*) x
   sen1 = seno(x, 3)
   sen2 = seno(x, 5)
   sen3 = sin(x)
   write (*,*) "Os tres resultados foram:"
   write (*,1) "Com 3 iteracoes:", sen1
   write (*,1) "Com 5 iteracoes:", sen2
   write (*,1) "Pela funcao sin do FORTRAN:", sen3
      1 format(A30, F19.15)
stop
end program senos

! funcao seno
!Esta função calcuma n termos da série infinita que define o seno
function seno(x, n)
real x, seno
integer n, i,fatorial,fa
   seno = 0
   do i = 0,n
     seno = seno + (pot(-1.,i)*pot(x, 1+2*i))/fatorial(1+2*i)
   end do
return
end function seno

! potencial
!Calcula o expoente exp de um número real x
function pot(x, exp)
real x, pot
integer exp, i
   pot = 1
   do i = 1, exp
     pot = pot*x
   end do
return 
end function pot

! fatorial
function fatorial(n)
integer n, fatorial, i
fatorial = 1
do i = 1, n
  fatorial = fatorial*i
  end do
return
end function fatorial