Немножко занимательной математики

13 views
Skip to first unread message

Leo B.

unread,
Dec 11, 2020, 8:53:58 PM12/11/20
to БЭСМ-6
Редкий случай, когда точности плавающей точки БЭСМ-6 хватает для впечатляющего результата (листинг слегка сокращен для удобства чтения):

 Ф O P T P A H
   /16.07.73/
                REAL FUNCTION MAGIC(F)
                I = F
       2        R = ALOG(I+6.)*(ALOG(I+6.)-1.)
       3        MAGIC = (F-I)*INT(R)+F
       4        RETURN
                END

                    BЫЗЫBAEMЫE ФYHKЦИИ И ПOДПPOГPAMMЫ

       INT         ALOG

                PROGRAM MAIN
                C = 2.635859741455
       2        DO 10 I =1,15
       3        N = C
       4        PRINT 1,N
       5    10  C = MAGIC(C)
             1  FORMAT(1XI5)
                END

                    BЫЗЫBAEMЫE ФYHKЦИИ И ПOДПPOГPAMMЫ

       MAGIC

           *EXECUTE
     2
     3
     5
     7
    11
    13
    17
    19
    23
    29
    31
    37
    41
    43
    47                                                      

А значащих цифр в мантиссе - всего 40.

Leo

Leo B.

unread,
Dec 11, 2020, 10:18:45 PM12/11/20
to БЭСМ-6
Если с двойной точностью, то в Фортране-Дубна и ГДР серьёзные ошибки, и даже форекс не способен самостоятельно ввести необходимую константу,
но работает ещё более впечатляющим образом:

F О R Е Х   ИПМ АН СССР           ВЕРСИЯ  2.02 ОТ  15.07.80

                 1       DОUВLЕ РRЕСISIОN FUNСТIОN МАGIС(F)                                            1
                 2       DОUВLЕ РRЕСISIОN F                                                            2
                 3       I = F                                                                         3
                 4       R = АLОG(I+6.)*(АLОG(I+6.)-1.)                                                4
                 5       МАGIС = (F-I)*INТ(R)+F                                                        5
                 6       RЕТURN                                                                        6
                 7       ЕND                                                                           7

F О R Е Х   ИПМ АН СССР           ВЕРСИЯ  2.02 ОТ  15.07.80

                 1       РRОGRАМ МАIN                                                                  1
                 2       DОUВLЕ РRЕСISIОN С                                                            2
                 3       RЕАD 22,С                                                                     3
                 4   22  FОRМАТ(D40)                                                                   4
                 5       РRINТ 2,С                                                                     5
                 6    2  FОRМАТ(1Х’С = ’D40.26)                                                        6
                 7       DО 10 I =1,25                                                                 7
                 8       N = С                                                                         8
                 9       РRINТ 1,N                                                                     9
                10   10  С = МАGIС(С)                                                                 10
                11    1  FОRМАТ(1ХI5)                                                                 11
                12       ЕND                                                                          12

           *ЕХЕСUТЕ
 С =         2.63585974145479126553562569+000
     2
     3
     5
     7
    11
    13
    17
    19
    23
    29
    31
    37
    41
    43
    47
    53
    59
    61
    67
    71
    73
    79
    83
    89
    97


Leo
Reply all
Reply to author
Forward
0 new messages