#a,i#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1;
D = vector(RDF,3000);
D[0]=2;
SumLogD=log(abs(D[0]))
for n in range(1,2001):
D[n]=0.5-(3/((D[n-1])^2));
D[n]=(D[n]^2)^(1/2);
SumLogD=SumLogD+log((D[n]));
L=(1/2000)*(SumLogD);
L;
0.896550885635 0.896550885635 |
#a,ii#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1;
max = P[0];
for n in range(1,2001):
if P[n]>max:
max=P[n];
min=P[0];
for n in range(1,2001):
if P[n]<min:
min=P[n];
binnumber=100;
rangevalue=max-min;
binsize=rangevalue/binnumber;
def f(x):
return (3/x)+(x/2)-1
B=vector(RDF,2000)
B[0]=1
for i in range (1,2000):
B[i]=f(B[i-1])
d=vector(RDF,1000)
for j in range(0,100):
d[j]=0
for j in range(0,2000):
d[floor((B[j])/binsize)] +=1
g=bar_chart(d)
show(g)
|
|
#b,i#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.5;
D = vector(RDF,3000);
D[0] = 2;
SumLogD=log(abs(D[0]))
for n in range(1,2001):
D[n]=0.5-(3/((D[n-1])^2));
D[n]=(D[n]^2)^(1/2);
SumLogD=SumLogD+log((D[n]));
L=(1/2000)*(SumLogD);
L;
0.896550885635 0.896550885635 |
#b,ii#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.5;
max = P[0];
for n in range(1,2001):
if P[n]>max:
max=P[n];
min=P[0];
for n in range(1,2001):
if P[n]<min:
min=P[n];
binnumber=100;
rangevalue=max-min;
binsize=rangevalue/binnumber;
def f(x):
return (3/x)+(x/2)-1.5
B=vector(RDF,2000)
B[0]=1
for i in range (1,2000):
B[i]=f(B[i-1])
d=vector(RDF,1000)
for j in range(0,100):
d[j]=0
for i in range(0,2000):
d[floor((B[j])/binsize)] +=1
g=bar_chart(d)
show(g)
|
|
#c,i#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.85;
D = vector(RDF,3000);
D[0] = 2;
SumLogD=log(abs(D[0]))
for n in range(1,2001):
D[n]=0.5-(3/((D[n-1])^2));
D[n]=(D[n]^2)^(1/2);
SumLogD=SumLogD+log((D[n]));
L=(1/2000)*(SumLogD);
L;
0.896550885635 0.896550885635 |
#c,ii#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.85;
max = P[0];
for n in range(1,2001):
if P[n]>max:
max=P[n];
min=P[0];
for n in range(1,2001):
if P[n]<min:
min=P[n];
binnumber=100;
rangevalue=max-min;
binsize=rangevalue/binnumber;
def f(x):
return (3/x)+(x/2)-1.85
B=vector(RDF,2000)
B[0]=1
for i in range (1,2000):
B[i]=f(B[i-1])
d=vector(RDF,1000)
for j in range(0,100):
d[j]=0
for j in range(0,2000):
d[floor((B[j])/binsize)] +=1
g=bar_chart(d)
show(g)
|
|
#d,i#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.92;
D = vector(RDF,3000);
D[0] = 2;
SumLogD=log(abs(D[0]))
for n in range(1,2001):
D[n]=0.5-(3/((D[n-1])^2));
D[n]=(D[n]^2)^(1/2);
SumLogD=SumLogD+log((D[n]));
L=(1/2000)*(SumLogD);
L;
0.896550885635 0.896550885635 |
#d,ii#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.92;
max = P[0];
for n in range(1,2001):
if P[n]>max:
max=P[n];
min=P[0];
for n in range(1,2001):
if P[n]<min:
min=P[n];
binnumber=100;
rangevalue=max-min;
binsize=rangevalue/binnumber;
def f(x):
return (3/x)+(x/2)-1.92
B=vector(RDF,2000)
B[0]=1
for i in range (1,2000):
B[i]=f(B[i-1])
d=vector(RDF,1000)
for j in range(0,100):
d[j]=0
for j in range(0,2000):
d[floor((B[j])/binsize)] +=1
g=bar_chart(d)
show(g)
|
|
show(g)#e,i#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.999;
D = vector(RDF,3000);
D[0] = 2;
SumLogD=log(abs(D[0]))
for n in range(1,2001):
D[n]=0.5-(3/((D[n-1])^2));
D[n]=(D[n]^2)^(1/2);
SumLogD=SumLogD+log((D[n]));
L=(1/2000)*(SumLogD);
L;
0.896550885635 0.896550885635 |
#e,ii#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.999;
max = P[0];
for n in range(1,2001):
if P[n]>max:
max=P[n];
min=P[0];
for n in range(1,2001):
if P[n]<min:
min=P[n];
binnumber=100;
rangevalue=max-min;
binsize=rangevalue/binnumber;
def f(x):
return (3/x)+(x/2)-1.999
B=vector(RDF,2000)
B[0]=1
for i in range (1,2000):
B[i]=f(B[i-1])
d=vector(RDF,1000)
for j in range(0,100):
d[j]=0
for j in range(0,2000):
d[floor((B[j])/binsize)] +=1
g=bar_chart(d)
show(g)
|
|
#f,i,c = 1.925#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.925;
D = vector(RDF,3000);
D[0] = 2;
SumLogD=log(abs(D[0]))
for n in range(1,2001):
D[n]=0.5-(3/((D[n-1])^2));
D[n]=(D[n]^2)^(1/2);
SumLogD=SumLogD+log((D[n]));
L=(1/2000)*(SumLogD);
L;
0.896550885635 0.896550885635 |
#f, i, c = 1.915#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.915;
D = vector(RDF,3000);
D[0] = 2;
SumLogD=log(abs(D[0]))
for n in range(1,2001):
D[n]=0.5-(3/((D[n-1])^2));
D[n]=(D[n]^2)^(1/2);
SumLogD=SumLogD+log((D[n]));
L=(1/2000)*(SumLogD);
L;
0.896550885635 0.896550885635 |
#f,ii, c = 1.915#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.915;
max = P[0];
for n in range(1,2001):
if P[n]>max:
max=P[n];
min=P[0];
for n in range(1,2001):
if P[n]<min:
min=P[n];
binnumber=100;
rangevalue=max-min;
binsize=rangevalue/binnumber;
def f(x):
return (3/x)+(x/2)-1.915
B=vector(RDF,2000)
B[0]=1
for i in range (1,2000):
B[i]=f(B[i-1])
d=vector(RDF,1000)
for j in range(0,100):
d[j]=0
for j in range(0,2000):
d[floor((B[j])/binsize)] +=1
g=bar_chart(d)
show(g)
|
|
#f,ii, c = 1.925#
P = vector(RDF,3000);
P[0]=1;
for n in range(1,2001):
P[n]=(3/P[n-1])+(P[n-1]/2)-1.925;
max = P[0];
for n in range(1,2001):
if P[n]>max:
max=P[n];
min=P[0];
for n in range(1,2001):
if P[n]<min:
min=P[n];
binnumber=100;
rangevalue=max-min;
binsize=rangevalue/binnumber;
def f(x):
return (3/x)+(x/2)-1.925
B=vector(RDF,2000)
B[0]=1
for i in range (1,2000):
B[i]=f(B[i-1])
d=vector(RDF,1000)
for j in range(0,100):
d[j]=0
for j in range(0,2000):
d[floor((B[j])/binsize)] +=1
g=bar_chart(d)
show(g)
|
|
#When we changed the value of c slightly, the Lyapunov exponent remains constant. The frequency distribution changed as the values of c change.#
|
|
#g: As c increases, the Lyapunov exonent remains constant while the frequency distribution increases.
|
|
|
|
|
|
|
|