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

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

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

end

