XProfan

Français

p.specht

Eines qui ersten Programme qui monde beschäftigte sich avec qui Berechnung qui sog. Bernoulli-payons [...]

, qui u.a. dans höherdimensionalen Reihenformeln une grand rôle spielen. qui Algorithmus stammt de 1843, aus einer Veröffentlichung de Lady Ada Lovelace, qui ersten Programmiererin qui monde, et était pour den ersten echten (Turing-mächtigen) ordinateur qui monde gedacht, qui "Analytical Engine" de Charles Babbage (... malheureusement nie fertiggestellt).

Kompilieren Marque Séparation

Titre de la fenêtre "Die BERNOULLI-ZAHLEN jusqu'à B_28 berechnen":Fenêtre Style 24:CLS
Set("decimals",17)
'*********************************************************************
'* This program computes Bernoulli number Bn using subroutine BERNOA *
'* ----------------------------------------------------------------- *
'* SAMPLE RUN:                                                       *
'*                                                                   *
'*  Compute Bernoulli number Bn for n = 0,1,...,10.                  *
'*                                                                   *
'*  Please enter Nmax: 10                                            *
'*                                                                   *
'*    n            Bn                                                *
'*  --------------------------                                       *
'*    0    0.100000000000D+01                                        *
'*    1   -0.500000000000D+00                                        *
'*    2    0.166666666667D+00                                        *
'*    4   -0.333333333333D-01                                        *
'*    6    0.238095238095D-01                                        *
'*    8   -0.333333333333D-01                                        *
'*   10    0.757575757576D-01                                        *
'*  --------------------------                                       *
'*                                                                   *
'* ----------------------------------------------------------------- *
'* REFERENCE: "Fortran Routines for Computation of Special Functions,*
'*             jin.ece.uiuc.edu/routines/routines.html".             *
'*                                 F90 Release By J-P Moreau, Paris. *
'*                                         (www.jpmoreau.fr)         *
'*                                                                   *
'* XProfan 11.2a-Version by P.Specht, Vienna/Austria                 *
'* No warranty whatsoever!                                           *
'*********************************************************************
'    PROGRAM MBERNOA
Déclarer B![200],n&,k&
Imprimer
nochma:
Imprimer "\n s'il te plaît Nmax eingeben: ";:Contribution N&

Si n&>29:Imprimer " trop grand!":goto "nochma":endif

    BERNOA(N&,B![])
    Imprimer "\n  n            Bn    "
    Imprimer "---------------------------"
    Imprimer "  ";0,Tab(11);Format $("%g",B![0])
    Imprimer "  ";1,Tab(11);Format $("%g",B![1])

    WhileLoop 2,N&,2:k&=&Boucle

        Imprimer "  ";k&,Tab(11);Format $("%g",B![k&])

    Endwhile

    Imprimer "---------------------------"
    WaitInput
    FIN

    Proc BERNOA :Paramètres N&,BN![]

        ' ======================================
        '  Purpose: Compute Bernoulli number Bn
        '  Contribution :  n --- Serial number
        '  Output:  BN(n) --- Bn
        ' ======================================
        '    IMPLICIT DOUBLE PRECISION (A-H,O-Z)
        Déclarer S!,R!,j&,k&,m&
        BN![0]=1
        BN![1]=-0.5

        WhileLoop 2,n&:m&=&Boucle

            S!=-(1.0e0/(M&+1)-0.5)

            WhileLoop 2,M&-1:k&=&Boucle

                R!=1

                WhileLoop 2,K&:j&=&Boucle

                    R!=R!*(J&+M&-K&)/J&

                Endwhile

                S!=S!-R!*BN![K&]

            Endwhile

            BN![M&]=S!

        Endwhile

        WhileLoop 3,N&,2:m&=&Boucle

            BN![M&]=0

        Endwhile

        Retour

    ENDPROC

    ' end of file mbernoa.prf