Inputs:
sigma
beta
rho
x0
y0
z0
epsilon