[mar 12] disks.py: demonstration of loopless programming in numpy