x1=45i; x2=36i; q=cos(x2*pi/180-x1*pi/180)^2; c=cos(x1*pi/180)^2*cos(x2*pi/180)^2+sin(x1*pi/180)^2*sin(x2*pi/180)^2