clear all;
clc;

A = randi([-9, 9], 20, 3);

for i = A'
  a = i(1);
  b = i(2);
  c = i(3);

  disp('');f
  disp(['a: ' num2str(a) ' b: ' num2str(b) ' c: ' num2str(c)]);

  if (a ~= 0)
    delta = b^2 - 4*a*c;
    if (delta > 0)
      %obliczanie pierwiastków
      x1 = (- b + sqrt(delta))/(2*a);
      x2 = (- b - sqrt(delta))/(2*a);
      disp(['x1: ' num2str(x1)]);
      disp(['x2: ' num2str(x2)]);

      %testowanie
      disp(['Wynik dla x1: ' num2str(a * x1^2 + b * x1 + c)]);
      disp(['Wynik dla x2: ' num2str(a * x2^2 + b * x2 + c)]);
    elseif (delta == 0)
      %obliczanie pierwiastka
      x0 = -b/(2*a);
      disp(['x0: ' num2str(x0)]);

      %testowanie
      disp(['Wynik dla x0: ' num2str(a * x0^2 + b * x0 + c)]);
    else
      disp(['Brak rozwiazan.']);
    endif
  elseif (b ~= 0)
    %obliczanie pierwiastka
    x0 = -c/b;
    disp(['x0: ' num2str(x0)]);

    %testowanie
    disp(['Wynik dla x0: ' num2str(b*x0 + c)]);
  elseif (c == 0)
    disp('Nieskoniecznie wiele miejsce zerowych.');
  else
    disp(['Funkcja stala y = ' num2str(c) '. Brak miejsc zerowych.']);
  endif
endfor