update of modeling script and addition of numerical integration script

modeling script - better visual, correction of scaling
numerical integration script - compare of relative errors of integration methods with Stefan - Boltzman law, temperature set for 5778 K