494 2
2
12. Rozwiązania zadań
8. (a) c(h)=ci{H=—■ sinnh—2n—^n2h2 -r~6nsh4 -... h
(b) c2 =4, c3 = 3n/3=5.I962, c6-6. Ekstrapolacja iterowana Richardsona daje c(0)% «6.2823, tzn. lepszą wartość od tej, którą Archimedes otrzymał za pomocą 96-bcku.
(c) c2b“4« sin =4łt \J — ^1 — cos — ^=^Snł — 4n \>4n2 — c2.
(d) Znoszenie się składników. Po zmianie instrukcji
x: •=sqrt(8 x n|2—4xnxsqrt(4 xC[0]|2)) na
x: = 2x C [0] x s?rt(0.5/(l + s*r/(l — C[0]/(n/2)T2)» otrzymuje się o wicie lepsze wyniki:
6 |
6.000000000 |
|
|
|
12 |
6.211657083 |
6.282209444 |
|
|
24 |
6.265257227 |
6.283123942 |
6.283184908 |
|
48 |
6.278700406 |
6.283181466 |
6.283185301 |
6.283185308 |
% |
6.282063902 |
6.283185068 |
6.283185308 |
6.283185308 |
§ 7.3
1. (a) Wybieramy x0 = 2, xt =4 (punkty wspólne dla obu wariantów).
/(*)» c0 + c x - 2)+ea(x - 2) (x - 4). c0 i Cj nie zależą od trzeciego punktu.
.\-c:2 ^ 2=c0 —Cq—2,
x=4 => I2«c0+c1(4—2) => c, =f 5.
Wartość c2 zależy od wyboru trzeciego punktu.
x = l b*. 0=2 + 5(l—2)+c2(l —2)(1 —4) =► c2 = l,
/(3)=2+5(3—2) + l *(3—2)(3—4)=6, x = 5 => 21 =2+5(5—2) + e2(5—2)(5—4) => c2=$.
/(3) = 2 + 5(3-2)+f(3-2)(3-4)«5.667.
Możemy stąd oszacować, że w przybliżeniu/(3) = 5.83+G.l7. Te same wyniki można oczy-wiście otrzymać z wzoru interpolacyjnego Newtona.
(b) Korzystając z wyniku (a), możemy przyjąć, te
/(x)=2 + 5(x-2)+l(x-2)(x-4)+c3(x-2)(*-4)(x-l),
/(5)=2I c3=^ *»/(3)= 5.833.
Wyszukiwarka
Podobne podstrony:
450 2 450 12. Rozwiązania zadań p(x) jcsi ostatnią wartością s. Ten sposób wymaga 2{n+ !) mnożeń i /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)454 2 454 12, Rozwiązania zadań (b) Z f=xyjz wynika, że Af Ax Av Az -T«— +---- f x y z Wprowadzamy456 2 456 12. Rozwiązania zadań 8. (a) c=(a2 + b2 — lab cos O1 2. c - wyznacza się w przybliżen458 2 458 12. Rozwiązania zadań 5. 0.5-10- °. 6. (a) 1.0«4”2łi p460 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 Aitk464 2 464 12. Rozwiązania zadań 4. (aj [ fj, j (0<y. fc^n). Jest to tzw. macier466 2 466 12. Rozwiązania zadań Dla/(x)=exp(x) na [- 1, 1] i A/=20 błąd maksymalny wielomianu interp468 2 468 12. Rozwiązania, zadań S(N) = YP/PP RETURN END (b) Poniższy program używa podprogramu470 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łasn476 2 476 12. Rozwiązania zadań 9 DO 2 1 = 2, N 10 IM 1=1-1 11 DO 2 K = 1,478 2 47B 12. Rozwiązania zadań Jest to równanie różnicowe o stałych współczynnikach, więc jego480 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żdej484 2 484 12. Rozwiązania zadań 1. Wybierzmy xN— 14 i x0—16 jako niewiadome. Równania486 2 486 12. Rozwiązaaia zadań (b) Ponieważ ATA jest macierzą symetryczną, więc w488 2 488 12. Rozwiązania zadań L: x3:~x2->2x(x2-x)}(j2—yI); >3: =/(x3); ifwięcej podobnych podstron