FORTRAN数值分析04 numerov算法

最新推荐文章于 2023-03-25 11:01:00 发布

River Chandler

最新推荐文章于 2023-03-25 11:01:00 发布

阅读量410

点赞数

分类专栏： # 科学计算原理文章标签：数据分析

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

版权

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

24 篇文章 2 订阅

订阅专栏

这段代码展示了Numerov算法的实现，用于数值求解常微分方程初值问题。程序首先读取初始条件，然后利用CPU时间计算并输出运行时间。尽管未提供完整解释，但可以看出算法的核心部分在于迭代更新变量y。最后，虽然注释提到会对比实际解，但实际比较部分被注释掉了。

摘要由CSDN通过智能技术生成

program numerov
    implicit none
    real f(10000000)
    real y(10000000)
    real reallyy(10000000)
    real,parameter :: h = 0.01
    real w0
    real w1
    real y0
    real dy0
    real yn
    real fn
    real w2
    integer count 
    integer i,j,k,p 
    real time_begin
    real time_end



    write(*,*)'input y0'
    read(*,*) y0
    y(1) = y0 
    
    write(*,*)'input dy0'
    read(*,*) dy0
    y(2) = y(1) + h * dy0

    call cpu_time(time_begin)

    do j = 1,10000000
        f(j) = 1
    end do 
    
    w0 = y0 * (1 + h**2/12 * f(1))
    w1 = y(2) * (1 + h**2/12 * f(2))
    yn = y(2)
    fn = f(1)

    do i=2,10000000
        w2 = 2 * w1 - w0 - h**2 * fn * yn
        fn = f(i)
        yn = w2 / (1 + h**2/12 * fn)
        y(i) = yn 
        w0 = w1
        w1 = w2
    end do 

    call CPU_time(time_end)

    !do k=1,100
        !write(*,*) y(k)
        !reallyy(k) = cos(h*k) + 2 * sin(h*k)
        !write(*,*) abs(reallyy(k) - y(k))
    !end do    
    
    write(*,*)time_end - time_begin
    read(*,*)
end program

上代码,关于numerov算法的话

之后会补上的