XProfan

Español

p.specht

Eines el ersten Programas el Welt beschäftigte se con el Berechnung el sog. Bernoulli-Pagar [...]

, el u.a. en höherdimensionalen Reihenformeln una große Papel spielen. Der Algorithmus stammt de 1843, de uno Veröffentlichung de Lady Ada Lovelace, el ersten Programmiererin el Welt, y war para el ersten echten (Turing-mächtigen) Computer el Welt pensamiento, el "Analytical Engine" de Charles Babbage (... desafortunadamente nie fertiggestellt).

Kompilieren Marca Separación

Título de la ventana "Die BERNOULLI-ZAHLEN a B_28 berechnen":Ventana de Estilo 24:CLS
Conjunto("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-Versión by P.Pájaro carpintero, Vienna/Austria                 *
'* No warranty whatsoever!                                           *
'*********************************************************************
'    PROGRAM MBERNOA
Declarar B![200],n&,k&
Imprimir
nochma:
Imprimir "\n Bitte Nmax eingeben: ";:Entrada N&

If n&>29:Imprimir " Zu groß!":goto "nochma":endif

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

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

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

    EndWhile

    Imprimir "---------------------------"
    WaitInput
    FIN

    Proc BERNOA :Parámetros N&,BN![]

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

        WhileLoop 2,n&:m&=&Loop

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

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

                R!=1

                WhileLoop 2,K&:j&=&Loop

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

                EndWhile

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

            EndWhile

            BN![M&]=S!

        EndWhile

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

            BN![M&]=0

        EndWhile

        Volver

    ENDPROC

    ' end of file mbernoa.prf