write me a method to calculate reaction wheel torque authority for a pyramid configuration of four reaction wheels in matlab

main.m
m = ...; % mass of satellite [kg]
r = ...; % distance from the center of mass [m]
w = ...; % angular velocity [rad/s]
118 chars
4 lines

main.m
Ixx = ...; % moment of inertia around x-axis [kg.m^2]
Iyy = ...; % moment of inertia around y-axis [kg.m^2]
Izz = ...; % moment of inertia around z-axis [kg.m^2]

I = [ Ixx   0    0  ;
        0  Iyy   0  ;
        0    0  Izz ];
230 chars
8 lines

main.m
h = I * w; % angular momentum [kg.m^2/s]
41 chars
2 lines

main.m
J = [ r/2   r/2   0    0   ;
        0     0   r/2  r/2 ;
       -r/4   r/4  r/4 -r/4 ];

T = J * h; % torque authority [N.m]
126 chars
6 lines