DataMuseum.dk

Presents historical artifacts from the history of:

RegneCentralen GIER Computer

This is an automatic "excavation" of a thematic subset of
artifacts from Datamuseum.dk's BitArchive.

See our Wiki for more about RegneCentralen GIER Computer

Excavated with: AutoArchaeologist - Free & Open Source Software. 

top - download

⟦1f0b33355⟧ Bits:30000632 Opg 698 S.Chr. 23.8.62 proceduren PK 10 (kommenteret), 8-hole paper tape

    Length: 2150 (0x866)
    Description: Bits:30000632 Opg 698 S.Chr. 23.8.62 proceduren PK 10 (kommenteret)
    Types: 8-hole paper tape
    Notes: Gier Text

GIER Text (HTML)

procedure PK10(xi,yi,uu,vu,a0,a1,a2,a3,a4,a5,a6,a7,a8,a9,a10);
_________


     value     xi,yi,      a0,a1,a2,a3,a4,a5,a6,a7,a8,a9,a10 ;
     _____


     real      xi,yi,uu,vu,a0,a1,a2,a3,a4,a5,a6,a7,a8,a9,a10 ;
     ____


     comment
     _______

     Idet z = xi + i ⨯ yi ,  w = uu + i ⨯ vu ,
     beregner PK10 værdien af polynomiet

          w  =  a10 ⨯ z∧10  +  a9 ⨯ z∧9  +  ...  +  a1 ⨯ z  +  a0 ,
                       |             |

     hvor koefficienterne er reelle og argumentet er komplext.
     Parametrene er:
     a0,...,a10    : koefficienterne i polynomiet,
     xi            : skal ved indhoppet indeholde realdelen af den uafh.var.
     yi            :  -    -      -         -     imaginærdelen -    -   - .
     uu            : indeholder ved udhoppet realdelen af polynomiets værdi.
     vu            :     -       -     -     imaginærdelen    -         -  .
     Metoden er angivet i
     C.Lanczos: Applied Analysis, Pitman, London, 1957,p.16;

begin
_____


     real h0,h1,h2,h3,h4,h5,h6,h7,h8,h9,h10,M,T;
     ____


          M   := -xi∧2-yi∧2              ;
                    |    |

          T   := 2⨯xi                    ;
          h10 := a10                     ;
          h9  := a9  + T ⨯ h10           ;
          h8  := a8  + T ⨯ h9  + M ⨯ h10 ;
          h7  := a7  + T ⨯ h8  + M ⨯ h9  ;
          h6  := a6  + T ⨯ h7  + M ⨯ h8  ;
          h5  := a5  + T ⨯ h6  + M ⨯ h7  ;
          h4  := a4  + T ⨯ h5  + M ⨯ h6  ;
          h3  := a3  + T ⨯ h4  + M ⨯ h5  ;
          h2  := a2  + T ⨯ h3  + M ⨯ h4  ;
          h1  := a1  + T ⨯ h2  + M ⨯ h3  ;
          h0  := a0            + M ⨯ h2  ;
          uu  := h1⨯xi + h0              ;
          vu  := h1⨯yi

end PK10 er testet paa DASK,
___

                             23.8.62.
                             Søren Christiansen;
[stop]