Lab3

//a) parzystosc liczby
liczba=22;

if modulo(liczba, 2)==0
    disp("Liczba jest parzysta.");
else
    disp("Liczba jest nieparzysta.");
end

//b) rozwiazania rownania kwadratowego ax^2+bx+c
a=-6; //wspolczynnik przy x^2
b=3; //wspolczynnik przy x
c=-1; //wyraz wolny

//miejsca zerowe
x1=0;
x2=0;

    delta = b*b - 4*a*c;

    if (delta > 0)
        x1 = (-b + sqrt(delta)) / (2*a);
        x2 = (-b - sqrt(delta)) / (2*a);
        disp ("Pierwiastki sa rozne.");
        disp (x1,"x1 = ");
        disp (x2,"x2 = ");
    elseif (delta == 0)
        disp("Dwa pierwiastki rowne (delta=0).");
        x1 = (-b + sqrt(delta)) / (2*a);
        disp (x1,"x1 = x2 = ");
    else
        rzeczywista = -b/(2*a);
        urojona =sqrt(-delta)/(2*a);
        disp("Delta <0. Pierwiastki zespolone.");
        disp ("i",urojona,"+",rzeczywista, "x1 = ");
        disp ("i",urojona,"-",rzeczywista, "x2 = ");
    end

//c) ciag Fibonacciego
dlugosc=20; //liczba elementow ciagu
n=dlugosc-1; //elementy numeruje sie od 0

a=0; //n-2 wyraz
b=1; //n-1 wyraz

W=zeros(1, n); //wektor z elementami ciagu

for i=0:n
    W(1,i+1)=a;
    tmp=a+b;
    a=b;
    b=tmp;
end

  disp(W," wyrazow ciagu Fibonacciego: ",dlugosc,"Wektor zawierajacy ");

//d) Silnia liczby naturalnej n
function [s]=silnia(n)
    if ( n ~= round(n) || n<0)
        error("Blad! Argument nie jest liczba naturalna.")
    elseif (n<2)
        s=n;
    else
        s=n*silnia(n-1);
    end
endfunction

//e) Pole i objetosc ostroslupa
function [p,v]=ost(a,b,h)
    Pp=a*b;
    h1=sqrt((a/2)^2+h^2); //wysokosc sciany 1
    p1=1/2*b*h1; //pole sciany 1
    h2=sqrt((b/2)^2+h^2); //wysokosc sciany 2
    p2=1/2*a*h2; //pole sciany 2
    Pb=2*p1+2*p2; //pole powierzchni bocznej
    v=1/3*Pp*h;
    p=Pp+Pb;
endfunction