“…PlotJoined ->True,PlotRange->AU, AxesLabel->{ "time (s)"," Q total (W/mA3)"}]; vv[i-,j,dtJ:= 11;If[Mod[numtimesteps,5]=0 II numtimesteps >=60, Print[numtimesteps," '',time," ",To [ 11," ",xv[ 13," ",phase[ 13," 'Iq[ l,time,qtot]," ",0rfC1],'~ ",qs [qtot,time]," ",qchem [ PlotJoined ->True, PlotRange->All, AxesLabel->{ "time (s)","average temp (K)" }]…”