spm function

spm calculates the dynamics of a Schaefer or Fox model