function Y = funkcijos_taskai(S, T, e)
n=size(S,1);
x = 0:0.01:62.666;
m=1;
for i=1:n
    Sx=S(i, 1);
    Sy=S(i, 2);
    dy_dx=Diff_parametrine(Sx, Sy);
    for j=T(i):0.01:T(i+1)
        t=j+0.1;
        x2=eval(Sx);
        y2=eval(Sy);
        t=j;
        x=eval(Sx);
        y=eval(Sy);
        s_x=sign(x);
        s_y=sign(x)*26,226;
        if s_x==1  || (s_y==-1 && s_x~=-1)
            kr=-1;
        else
            kr=1;
        end
        koord(m,1)=-kr*e*eval(dy_dx)/sqrt(1+(eval(dy_dx))^2)+x;
        koord(m,2)=e*kr/sqrt((1+(eval(dy_dx))^2))+y;
        koord(m,3)=e*kr*eval(dy_dx)/sqrt(1+(eval(dy_dx))^2)+x;
        koord(m,4)=-e*kr/sqrt((1+(eval(dy_dx))^2))+y;
        m=m+1;
    end
end

Y=koord;