fortran：使用IMSL函数库计算第一/二类完全椭圆积分

IMSL函数库是一款强大的商业函数库，高校学生或是老师可以使用教育邮箱进行申请。下面给出IMSL的官网
https://www.roguewave.com/products-services/imsl-numerical-libraries

接下来给出Fortran使用IMSL函数库计算第一类椭圆积分的示例代码
program test
        include 'link_fnl_shared.h'
        use elk_int
        use ele_int
        
        real*8 :: x, vaule 
        
        x     = 5d-1
        value = d_elk(x)  !// 计算第一类完全椭圆积分
        write(*,'(g0)') value
        
        x     = 5d-1
        value = d_ele(x)  !// 计算第二类完全椭圆积分
        write(*,'(g0)') value
        
        !// 与matlab计算结果相同
        !// M = 0.5;
        !// [K,E] = ellipke(M)
        !// K = 1.8541
        !// E = 1.3506
        
end program test