gistfile1.py

from math import sin, cos, pi


l = list()
for alpha in np.linspace(0,90,90):
    alpha = alpha*pi/180         # plane angle in radians
    R = np.matrix([ [cos(alpha) , 0, sin(alpha)], 
                    [0          , 1,          0],
                    [-sin(alpha), 0, cos(alpha)] ])

    rotated = np.matrix(data[['X','Y','Z']])*R
    l.append(rotated.std(axis=0))

standard_devs = pd.DataFrame(np.vstack(l), columns=['std(X)','std(Y)','std(Z)'])    
standard_devs.plot()