HI I CAN HELP. Give me the details
HERE IS A LITTLE SAMPLE OF PREVIOUS WORK
void equation8(float theta,float gamma, float omega_g,double c[3][3])
{
c[0][0]=cos(gamma)*cos(omega_g) - sin(gamma)*sin(omega_g)*sin(theta); c[0][1]=-cos(theta)*sin(omega_g); c[0][2]=cos(omega_g)*sin(gamma) + cos(gamma)*sin(omega_g)*sin(theta);
c[1][0]=cos(gamma)*sin(omega_g) + cos(omega_g)*sin(gamma)*sin(theta); c[1][1]=cos(omega_g)*cos(theta); c[1][2]=sin(gamma)*sin(omega_g) - cos(gamma)*cos(omega_g)*sin(theta);
c[2][0]= -cos(theta)*sin(gamma); c[2][1]=sin(theta); c[2][2]=cos(gamma)*cos(theta);
}
void equation10(float theta,float gamma, float omega,double c[3][3])
{
c[0][0]=cos(gamma)*cos(omega) - sin(gamma)*sin(omega)*sin(theta); c[0][1]=-cos(theta)*sin(omega); c[0][2]=cos(omega)*sin(gamma) + cos(gamma)*sin(omega)*sin(theta);
c[1][0]=cos(gamma)*sin(omega) + cos(omega)*sin(gamma)*sin(theta); c[1][1]=cos(omega)*cos(theta); c[1][2]=sin(gamma)*sin(omega) - cos(gamma)*cos(omega)*sin(theta);
c[2][0]= -cos(theta)*sin(gamma); c[2][1]=sin(theta); c[2][2]=cos(gamma)*cos(theta);
}
void equation11(float theta,float theta_r,float gamma,float gamma_r, float omega, float omega_r,double Wb[3][3])
{
Wb[0][0]=cos(gamma*theta_r)-(sin(gamma)*cos(theta*omega_r));
Wb[1][0]=gamma_r+sin(theta*omega_r);
Wb[2][0]=sin(gamma*theta)+(cos(gamma)*cos(theta*omega_r));
}