cos(a-b)=cosacosb+sinasinb(1)
cos(a+b)=cosacosb-sinasinb(2)
(1)-(2)
cos(a-b)-cos(a+b)=2sinasinb
sinasinb=[cos(a-b)-cos(a+b)]/2