Implementing Simpson's Rule in Fortran for Accurate Numerical Integration

[Daniel1992]

I am trying to program Simpson's Rule in Fortran

I=(1/3)*h*(y0 + 4*y1 + 2*y2 + 4*y3 + 2*y4 + ... 2*yn-2 + 4*yn-1 + yn)

n=100 (number of strips)

I have generated the y values but I don't know how to get Fortran to follow the pattern in the brackets after h to add up the y values.

Any help would be appreciated.

 

[SteamKing]

Well, do you notice any patterns about which y-values take which Simpsons's multipliers?

 

[Daniel1992]

Well, do you notice any patterns about which y-values take which Simpsons's multipliers?

A part from the first and last y values the odd ones are multiplied by 4 and the even are multiplied by 2.

 

[SteamKing]

Does this pattern suggest how you might go about accumulating products of the y-values and the corresponding Simpson's multipliers?

 

[Daniel1992]

Does this pattern suggest how you might go about accumulating products of the y-values and the corresponding Simpson's multipliers?

I can't work it out :(

 

[Redbelly98]

When you write a do-loop, you can make the loop variable increase by 2 each time through the loop.