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

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

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

end

