function [wynik] = solvesin(a, b, c)
  k=asin(c);
  l=a-modulo(a, 2*%pi) - 2*%pi;
  wynik1 = l + k;
  wynik2 = l + %pi - k;
  n=1;
  if sin(wynik1) ~= sin(wynik2)
    while wynik1(n) <= b+2*%pi
      n = n + 1;
      wynik1(n) = wynik1(n-1) + 2*%pi;
    end
    n=1;
    while wynik2(n) <= b
      n = n + 1;
      wynik2(n) = wynik2(n-1) + 2*%pi;
    end
    wynik = [wynik1; wynik2];
  else
    while wynik1(n) <= b+2*%pi
      n = n + 1;
      wynik1(n) = wynik1(n-1) + 2*%pi;
    end
    wynik=wynik1;
  end
  wynik = wynik(find(wynik <= b & wynik >= a));
  wynik = gsort(wynik, 'g', 'd');
endfunction

function [wynik]=solvecos(a, b, c)
  k=acos(c);
  l=a-modulo(a, 2*%pi)-2*%pi;
  wynik1 = l - k;
  wynik2 = l - 2*%pi + k;
  n=1;
  if cos(wynik1) ~= cos(wynik2)
    while wynik1(n) <= b+2*%pi
      n = n + 1;
      wynik1(n) = wynik1(n-1) + 2*%pi;
    end
    n=1;
    while wynik2(n) <= b
      n = n + 1;
      wynik2(n) = wynik2(n-1) + 2*%pi;
    end
    wynik = [wynik1; wynik2];
  else
    while wynik1(n) <= b+2*%pi
      n = n + 1;
      wynik1(n) = wynik1(n-1) + 2*%pi;
    end
    wynik=wynik1;
  end
  wynik = wynik(find(wynik >=a & wynik <= b));
  wynik = gsort(wynik, 'g', 'd');
endfunction

function [wynik] = solvetan(a, b, c)
  k=atan(c);
  l=a-modulo(a, %pi)- %pi;
  wynik1 = l - %pi + k;
  n=1;
    while wynik1(n) <= b + %pi
    n = n + 1;
    wynik1(n) = wynik1(n-1) + %pi;
  end
  wynik=wynik1;
  wynik = wynik(find(wynik <= b & wynik >= a));
  wynik = gsort(wynik, 'g', 'd');
endfunction