[X,Y]=taylorova_metoda('sin(x^2*y)-y', 1, 0, 4, 0.1, 4, 1);