|
From: Prahas D. N. <pra...@gm...> - 2015-03-29 17:07:17
|
Hi Mat-Plotters,
I'm trying to modify the below code so that I can
set the initial conditions to (-1,0,0.5).
The code below randomly sets the initial conditions:
**************
# I changed the equation -- it's not Lorenz.
N_trajectories = 1
def lorentz_deriv((x, y, z), t0, aa=1.1, yy=0.87):
"""Compute the time-derivative of a Lorentz system."""
return [y*(z-1+x*x)+yy*x, x*(3*z+1-x*x)+yy*y, -2*z*(aa+x*y)]
# Choose random starting points, uniformly distributed from -15 to 15
np.random.seed(1)
*# Here's the statement which assigns the initial conditions:*
x0 = -15 + 30 * np.random.random((N_trajectories, 3))
******************
I tried simply doing this:
x0 = (-1,0,0.5)
but I get this error:
ValueError: need more than 1 value to unpack
What am I missing? What is the correct way to make the assignment?
Thanks!
--Prahas
|