sigmaS=0.01
sigmabeta=0.1
theta=60*pi/200
S=254
omega=200/pi
Kb=[power(sigmaS,2) 0 ; 0 power(sigmabeta,2)]
G=[sin(theta) S*cos(theta)/omega; cos(theta) -S*sin(theta)/omega]
KEN=G*Kb*transpose(G)
corr=KEN(1,2)/sqrt(KEN(1,1)*KEN(2,2))
corr=-0.0757/sqrt(0.0551*0.1042)
