read SMAZ: DIAGpaper(1,(1-x)^2+(1-y)^2+(1-z)^2-2,a,[x,y,z],n,{a},F,A); print(`This took`, time(), `seconds. `): quit: