<< nsq.m m=100; (* Number to represent *) r=10; (* Number of squares *) Print["m=",m," Number to represent"]; Print["r=",r," Number of squares"]; n=Ceiling[ Sqrt[m] ]; din = Date[]; coeff = Coefficient[ sq[n,q,r], q, m ]; dout = Date[]; Print[ " Time in: ", din ]; Print[ "Representations: ", coeff ]; Print[ " Time out: ", dout ]; tin = din[[6]] + 60 (din[[5]] + 60 * (din[[4]] + 24 * din[[3]])); tout = dout[[6]] + 60 (dout[[5]] + 60 * (dout[[4]] + 24 * dout[[3]])); time = tout - tin; Print[ " Time: ", time ]; (* Okay if both times are in the same month *)