Skip to content

Instantly share code, notes, and snippets.

Created May 5, 2017 17:36
Show Gist options
  • Star 0 You must be signed in to star a gist
  • Fork 0 You must be signed in to fork a gist
  • Save anonymous/434025eb71084b77b8f6ed53bab0f17b to your computer and use it in GitHub Desktop.
Save anonymous/434025eb71084b77b8f6ed53bab0f17b to your computer and use it in GitHub Desktop.
Matrix3D M_2=Area_radius_Vector(Minutes_C -KIntercalary*(Days_C+1)-350*(IntercalaryYear-1),Hours_C,Days_C,Astronom_K);
Matrix3D M_0=(Rotation_axis(Minutes_C,Hours_C,Astronom_K))*inclinationPlanet_23_angle()*M_2;
Matrix3D M_1=M_rotation_Z(90-W_Fi)*M_0;
Matrix3D M_Inv=new Matrix3D(1,0,0,0,0,1,0,0,0,0,1,0,0,0,(RadiusEarth/K_RadiusEarth),1)*M_rotation_Z(90-W_Fi)*M_0;
M_Inv.Invert();
/// <summary>
/// вектор центр Сонця центр Землі
/// </summary>
Vector3D V_Solar=new Vector3D(M_2.OffsetX/100,M_2.OffsetY/100,M_2.OffsetZ/100);
/// <summary>
///
/// </summary>
Vector3D V_Earth=MultiplyMatrixVector(M_Inv,V_Solar);
/// <summary>
///
/// </summary>
Fi_Solar_Earthx=V_Earth.X;
Fi_Solar_Earthy=V_Earth.Y;
Fi_Solar_Earthz=V_Earth.Z;
GeometrFigure1.Transform=new MatrixTransform3D(M_0);
GeometrFigure2.Transform=new MatrixTransform3D(M_0);
GeometrFigure3.Transform=new MatrixTransform3D((M_1));
GeometrFigure4.Transform=new MatrixTransform3D((Rotation_axis(Minutes_C,Hours_C,Astronom_K))*M_2);
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment