Write library for effective area calculation about matlab-system-model

[Beredos2]

Here's a verbalized algorithm for calculating the effective area of one CubeSat face relative to any celestial body.

1. Find the satellite's position relative to the Earth-centered-inertial (ECI) frame: this is either provided by ADCS, or must be derived from orbital elements.
2. Find the position of the celestial body relative to the ECI frame.
3. Subtract the celestial body position from the spacecraft position, which gives the spacecraft-body vector.
4. Normalize the spacecraft-body vector (such that its magnitude is 1)
5. Find the area vector relative to the spacecraft-fixed frame of reference (the area vector is normal to the surface of interest and outward pointing)
6. Calculate the angle between the area vector and the spacecraft-body vector. If angle is greater than pi (180 degrees), the effective area is zero (that face of the spacecraft is in shadow). If the angle is less than 180, the effective area is equal to the true area times the point product of the two vectors.