function rhs = GlobalEquiriblium(X,edd, wx,wy,wz,Cd)
rx=X(1);
ry=X(2);
rz =X(3);
kx=X(4);
kz=X(5);
rhs(1)= -(kz^2-kx^2)-3*Cd*edd*((kx^2-kz^2)-(2*kx^2+kz^2)*f(kz/kx,kz/kx)/2)/(rx*ry*rz);
rhs(2)= rx*ry*rz*kx*kx*kz-1;
rhs(3)= rx^2-(2*kx^2+kz^2)/3 ...
    -Cd*edd*(f(kz/kx, kz/kx)- f(rx*wz/(rz*wx),ry*wz/(rz*wy)) ...
    +rx*wz/(rz*wx)*f1(rx*wz/(rz*wx),ry*wz/(rz*wy)))/(rz*ry*rx);
rhs(4)=ry^2-(2*kx^2+kz^2)/3 ...
    -Cd*edd*(f(kz/kx, kz/kx)- f(rx*wz/(rz*wx),ry*wz/(rz*wy)) ...
    +ry*wz/(rz*wy)*f2(rx*wz/(rz*wx),ry*wz/(rz*wy)))/(rz*ry*rx);
rhs(5)=rz^2-(2*kx^2+kz^2)/3 ...
    -Cd*edd*(f(kz/kx, kz/kx)- f(rx*wz/(rz*wx),ry*wz/(rz*wy)) ...
    -rx*wz/(rz*wx)*f1(rx*wz/(rz*wx),ry*wz/(rz*wy)) ...
    -ry*wz/(rz*wy)*f2(rx*wz/(rz*wx),ry*wz/(rz*wy)))/(rz*ry*rx);