function [A] = rotatez(angle)
%rotate around x axis
%A=rotatex("angle in radians")

co=cos(angle);
si=sin(angle);

A=[co,-si,0;si,co,0;0,0,1];

end

