Fortran Duda Sumatorios Fortran

[Becquerel]

Hola!

Esto es una chorrada para los que tienen práctica programando, pero yo soy poco hábil y me atasco en estas cosas.

Estoy haciendo integración numérica con Fortran y tengo que hacer un programa que me evalúe una integral de resultado conocido en un intervalo. La fórmula para calcular la integral es:

$\frac{h}{2}\left[ f(a) + 2\sum_{i = 1}^{n-1} f(a+i\cdot b ) + f(b)\right]$

Donde f(a) y f(b) son los valores de la función en los extremos de los intervalos.

Mi problema es que no sé cómo hacer el sumatorio. Por ejemplo, si tengo la función $cos x$, ¿cómo implemento el sumatorio con bucles?

Seguro que es una tontería, pero no me sale. Muchas gracias,

Un saludo

[Al2000]

Re: Duda Sumatorios Fortran

Inicializa la suma (digamos S) a cero, acontinuación pon un bucle DO I=1,n-1. En el cuerpo del DO incrementa el valor de la suma, S = S + f(a + I * b). A continuación del bucle multiplica por 2 y suma los valores extremos de la función, S = 2 * S + f(a) + f(b).

Te pusiera el código, pero ya tengo el FORTRAN algo olvidado. Queda algo así, pero corrige la sintaxis:

S = 0
DO I = 1,N-1
  S = S + F(A + I * B)
END
S = 2 * S + F(A) + F(B)

Saludos,

Al

PD. Cuando te respondí iba saliendo a trabajar y no me percaté que tienes un error en el argumento de la función dentro de la sumatoria. Debería decir $f(a + i \cdot h)$.

[Becquerel]

Re: Duda Sumatorios Fortran

Cierto, es verdad, me he colado escribiendo eso.

Muchas gracias Al!! Ya he resuelto las dudas que tenía. Esto de calcular los sumatorios numéricamente me trae un poco de cabeza.

Un saludo!!