Sprawozdanie – Jerzy Jurczak nr albumu 13578

Polecnie do zadania 1 było następujące:

Wykorzystujac operatory arytmetyczne Scilaba rozwiaz układ równan

x + 2 y + 3 z = 1

2 x + 3 y + z = 3

3 x + y + 2 z = 2

Wartości stojące przy niewiadomych wypisujemy i tworzymy z nich macierz A:

A =[1,2,3

2,3,1

3,1,2]

Analogicznie postępujemy z wartościami będącymi wynikami równań tworząc macierz B:

B = [1,

3,

2]

Wynik układu równań otrzymujemy korzystając z polecenia: X=A/B

Zadanie 2

Należało napisać funkcje MojePi o argumencie n, zwracającą przybliżenie wartości liczby za pomocą:

a) iteracyjnego obliczania sumy szeregu,

b) utworzenia wektora kwadratów odwrotności kolejnych liczb naturalnych, a następnie zsumowania elementów

tego wektora (polecenie sum).

Iteracyjnego obliczenia szeregu sumy dokonujemy poprzez utworzenie funkcji wykorzystującej pętlę for:

function x=MojePi(n)

suma=0;

for i=1:n do

suma=suma+1/(i^2);

end

x=sqrt(6*suma);

Podpunkt drugi opiera się na utworzeniu zmiennej wektor (wek) o zakresie od 1 do n a następnie obliczenia sumy opisanej w zadaniu:

function x=MojePi1(n)

wek = 1:n;

suma=sum(wek.^-2);

x=sqrt(6*suma);

W zadaniu 3 należało utworzyć skrypt wykonujący wykres pierwszych 100 liczb harmonicznych, liczonych według wzoru podanego w zadaniu.

Tworzymy funkcję H odpowiadającą za zliczanie sumy:

function suma=H(n)

suma=0;

for k=1:n do

suma=suma+1/(k);

end

endfunction

Następnie za pomocą pętli for tworzymy ciąg pierwszych 100 liczb harmonicznych:

for i=1:100 do

x(i) =i

harm(i)=H(i)

logar(i)=log(i)

end

Wykresy kolejno liniowo-liniowy i logarytmiczno-liniowy tworzymy poprzez polecenia:

plot2d(x,[harm,logar],[2,3])

plot2d('ln',x,[harm,logar],[2,3])