452 2

452 2



452


12. Rozwiązania zadań

Jeśli z, =0. to nwd (r_łt r0) = nwd (x, y)=r0-y. Zauważmy,że

(i)    dowolna liczba dzieląca r_x i r0 musi dzielić także r,,

(ii)    odwrotnie, dowolna liczba dzieląca r, i r0 musi dzielić także r_ ,.

To sugeruje następujący algorytm (Euklidesa). Tworzymy ciągi skończone {r„} \    ^

+ *■«    (&><>. r,n = 2, 3.....M.

Obliczenia są skończone, gdyż /'B+l<rll dla «>0. czyli musi być /-„-O. np. dla * = >-. Jeśli jednak rv = 0. to *>,,[/•*_ 2 i r,v|r^_ j (w istocie jest rv_x=nwd (rA_2, /•*_,)). (cź 3 wobec (ii). Rozumując tak dalej, otrzymujemy, że rK_v|r0 i rv_,|r_1. To, żc rs-1 jest największym wspólnym dzielnikiem r0 i r_. wynika stąd, że dzięki (i) dowolny dzielnik r0 i r, musi dzielić rN_,.

Poniższy program w języku Basic jest rozwiązaniem ostatniej części zadania. Działa on dla dowolnego ułamka 11/12.

10 PR1NT "WPROWADZIĆ LICZBY CAŁKOWITE Tl I 12, GDZIE 12#0"

20 INPUT II, 12 30 S-SGN(I1)*SGN (12)

40 RI = ABS(Il)

50 R2= ABS(I2)

60 Q=INT(RI/R2)

70 PRINT "CZESC CAŁKOWITA =", S*Q 80 R1 = RI — Q*R2

90 REM STOSUJE SIE TERAZ ALGORYTM EUKLIDESA (WIERSZE 130- 170) 100 REM DO UPROSZCZENIA CZĘŚCI UŁAMKOWEJ (=S*(RI/R2))

110 Tl = Rl 120 T2=R 2

130 R3 = RI — INT (R1/R2)*R2 140 IF R3 = 0 TI1EN ISO 150 Rl — R2 160 R2 = R3 170 GOTO 130

180 PRrNT "CZESC UŁAMKOWA-', S»T!/R2, 7*\ T2/R2 190 END

(Funkcja INT(X/Y) daje największą liczbę całkowitą nic większą od X/Y.) Czytelnik może wykonać ten program na maszynie lub ręcznie, aby upewnić się, że daje on nastę pujące wyniki.

11

12

Część całkowita

Część ułamkowa

5

4

1

1/4

162

54

3

0/1

172

22

7

9/11

-44

154

0

— 2/7

-7

7

-l

on

— 65

— 3

21

2/3


Wyszukiwarka

Podobne podstrony:
478 2 47B 12. Rozwiązania zadań Jest to równanie różnicowe o stałych współczynnikach, więc jego
464 2 464 12. Rozwiązania zadań 4. (aj [ fj,    j (0<y. fc^n). Jest to tzw. macier
500 2 5f*U 12. Rozwiązania zadań (ft) 7hieżność jest szybsza, jeśli a -1 2 sin tzn. jeśli (I -+ 2>
504 2 504 12. Rozwiązania zadań. if abs (R)<e.px then go to out; T[k, /+!]:« H*. r]+i?/(2tf/+l)-l
450 2 450 12. Rozwiązania zadań p(x) jcsi ostatnią wartością s. Ten sposób wymaga 2{n+ !) mnożeń i /
454 2 454 12, Rozwiązania zadań (b) Z f=xyjz wynika, że Af Ax Av Az -T«— +---- f x y z Wprowadzamy
456 2 456 12. Rozwiązania zadań 8. (a) c=(a2 + b2 — lab cos O1 2. c - wyznacza się w przybliżen
458 2 458 12. Rozwiązania zadań 5.    0.5-10- °. 6.    (a) 1.0«4”2łi p
460 2 460 12. Rozwiązania zadań§ 3.2 1.    (a) 0.693: (b) około 1000. 2.
462 2 462 12. Rozwiązania zadań Iloraz kolejnych błędów jest więc stały i dlatego ekstrapolacja Aitk
466 2 466 12. Rozwiązania zadań Dla/(x)=exp(x) na [- 1, 1] i A/=20 błąd maksymalny wielomianu interp
468 2 468 12. Rozwiązania, zadań S(N) = YP/PP RETURN END (b) Poniższy program używa podprogramu
470 2 470 12. Rozwiązania zadań § 4.5 1.    Dla £(/)=$>,/00 wybieramy jako/dowolną
474 2 474 12. Rozwiązania zadań 4. (a) Utworzyć i porównać Ax i )jg. Dla wielokrotnej wartości własn
476 2 476 12. Rozwiązania zadań 9 DO 2 1 = 2, N 10 IM 1=1-1 11 DO 2 K = 1,
480 2 480 12. Rozwiązania zadań (b)    <łi=/<-‘<SA,
482 2 482 12. Rozwiązania zadań i używamy metody Gaussa-Seidela, tj. ostatniego przybliżenia każdej
484 2 484 12. Rozwiązania zadań 1. Wybierzmy xN— 14 i x0—16 jako niewiadome. Równania
486 2 486 12. Rozwiązaaia zadań (b)    Ponieważ ATA jest macierzą symetryczną, więc w

więcej podobnych podstron