for i=[1:100] do
    kat=w*t(i);
    T1=[cos(kat), 0, -sin(kat); 0, 1, 0; sin(kat), 0, cos(kat);];
    Wspolrzedne1=[xPS(i); yPS(i); zPS(i)];
    Wspolrzedne1=T1*Wspolrzedne1;
    xPS(i)=Wspolrzedne1(1);
    yPS(i)=Wspolrzedne1(2);
    zPS(i)=Wspolrzedne1(3);

    T=[cos(kat), -sin(kat), 0; sin(kat), cos(kat), 0; 0, 0, 1];
    Wspolrzedne=[xPS(i); yPS(i); zPS(i)];
    Wspolrzedne=T*Wspolrzedne;
    xPS(i)=Wspolrzedne(1);
    yPS(i)=Wspolrzedne(2);
    zPS(i)=Wspolrzedne(3);
  end