FORTRAN数值分析01求定积分与微分(VScode环境下)

最新推荐文章于 2025-02-04 22:04:34 发布

River Chandler

最新推荐文章于 2025-02-04 22:04:34 发布

阅读量4.9k

点赞数 8

分类专栏： # 科学计算原理文章标签：其他

本文链接：https://blog.csdn.net/Chandler_river/article/details/124359777

版权

科学计算原理专栏收录该内容

24 篇文章

订阅专栏

1.矩形法求定积分

program definiteintegralbyrectangular
    implicit none
    real,parameter::pi=3.1415926
    integer,parameter::N=100!调整精度
    real lowbound,upbound,len,x,sum,f
    lowbound = 0 !输入下界
    upbound  = pi!输入上界
    len = (upbound - lowbound)/N 
    sum = 0.0
    x = lowbound
    do while(x<upbound)
        f = cos(x)-x**2 !输入要积分的函数
        sum = sum + f * len
        x = x + len
    end do
    write(*,*)'The final result = ',sum
    pause
    end

2.梯形法求定积分

program definiteintegralbyrectangular
    implicit none
    real,parameter::pi=3.1415926
    integer,parameter::N=100!调整精度
    real lowbound,upbound,len,x,sum,f1,f2
    lowbound = 0 !输入下界
    upbound  = pi!输入上界
    len = (upbound - lowbound)/N 
    sum = 0.0
    x = lowbound
    do while(x<upbound)
        f1 = cos(x)-x**2 !输入要积分的函数
        f2 = cos(x+len)-(x+len)**2
        sum = sum + (f1+f2)*len/2
        x = x + len
    end do
    write(*,*)'The final result = ',sum
    pause
    end

3.辛普森法

$\int_{a}^{b}f(x)dx\approx \sum_{x=a}^{b}[f(x)+4f(x+\frac{\Delta x}{2})+f(x+\Delta x)]*\Delta x/6$

program definiteintegralbyrectangular
    implicit none
    real,parameter::pi=3.1415926
    integer,parameter::N=100!调整精度
    real lowbound,upbound,len,x,sum,f1,f2,f3
    lowbound = 0 !输入下界
    upbound  = pi!输入上界
    len = (upbound - lowbound)/N 
    sum = 0.0
    x = lowbound
    do while(x<upbound)
        f1 = cos(x)-x**2 !输入要积分的函数
        f2 = (cos(x+len/2)-(x+len/2)**2)*4
        f3 = cos(x+len)-(x+len)**2
        sum = sum + (f1+f2+f3)*len/6
        x = x + len
    end do
    write(*,*)'The final result = ',sum
    pause
    end