r code to run the webworld model of calderelli et al 1998