Numerical Analysis Programs Using Fortran 90 - University of ...

Problem: Write a Fortran program to find the value of the integral
π
2
dx
∫ 2
sin x + cos
0
dimension y(20)
f(x)=1/((sin(x))**2+1./4*(cos(x))**2)
pi=22./7
a=0; b=pi/2; n=10
h=(b-a)/n
do i=1,n+1
x=a+h*(i-1)
y(i)=f(x)
print*,y(i)
enddo
t=0
do i=1,n
t=t+h*(y(i)+y(i+1))/2
enddo
print*,'Integration by Trapezoidal rule =',t
end
1
4
2
x
using Trapezoidal rule with n = 10
^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
Integration by Trapezoidal rule = 3.142227
34