The physical data used for simulation can be found at : https://drive.google.com/drive/folders/194ralzQNBT4iwupcq72sTHUX19cXKTNB?usp=sharing
For more information on the physics and some sample results, see: https://confluence.mcmasterneudose.ca/x/BAChAg
#private contains matlab functions needed for the ode solvers to work, it needs to be there for the simulations to run properly
#ECEF_basic_propagation.m is the main code that propagates the attitude, it also plots the pointing error, angular velocities and all modelled torques on the satellite
#eqset.m Sets up the 10 coupled ODES needed to propagate the attitude, angular velocity and hysteresis model (flatley). This utilizes many of the other functions
#ode113quat and ode420 are modified versions of ode113 and ode45 which include a normalization step for the quaternions
#CalcT caluclates all torques on the satellite
#calcBeta caluclates the pointing error of the satellite
#PlotT calculates all the torques for plotting