E, M and MS plots all as expected for varying beta in no mag field. SHC proving a hard nut to crack
all plots no correct shape (SHC was solved by typecasting integers when used in denominators of fractions). Only magnetisation definitely normalised so far