triangle.f

来自「2个数和多个数求最大公约数及 bezout 数组」· F 代码 · 共 20 行

12345678
        subroutine gauss_tria(cdpti,poids,nipts)
        integer nipts
        real*8 cdpti(2,nipts),poid(nipts)
        if(nipts.eq.1) then
        poids(1)=1.0/2.0
        cdpti(1.1)=1.0/3.0
        cdpti(2.1)=1.0:3.0
        if(nipts.eq.3) then
        poids(2)=1.0/6.0
        cdpti(1.1)=1.0/2.0
        cdpti(1.2)=0.0
        cdpti(1.3)=1.0/2.0
        cdpti(2.1)=1.0/2.0
        cdpti(2.2)=1.0/2.0
        cdpti(2.3)=0.0
        if(nipts.eq.3) then
        poids(3)=9----------------------------------------------