|
From: Prahas D. N. <pra...@gm...> - 2015-03-30 15:45:30
|
Hi, For those of you following the trials and tribulations of moi, I hacked the solution. The assignment is: x0 = [ [-1.0,0.0,0.5] ] I printed the orig x0. Printed mine. Noticed mine was missing a set of brackets. Tried it. Success! --Prahas On Sun, Mar 29, 2015 at 10:07 AM, Prahas David Nafissian < pra...@gm...> wrote: > > 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 > > > > > |