#OK to post homework

#Nicholas DiMarzio, 11/1/21, Assignment 16

#
#Problem 1

Orbk(2, x, x[1]*(5/3 - x[2]), [0, 1], 1000, 1020);
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

Orbk(2, x, x[1]*(5/3 - x[2]), [0.5, 1], 1000, 1020);
   [0.6666666667, 0.6666666667, 0.6666666667, 0.6666666667, 

     0.6666666667, 0.6666666667, 0.6666666667, 0.6666666667, 

     0.6666666667, 0.6666666667, 0.6666666667, 0.6666666667, 

     0.6666666667, 0.6666666667, 0.6666666667, 0.6666666667, 

     0.6666666667, 0.6666666667, 0.6666666667, 0.6666666667, 

     0.6666666667]


Orbk(2, x, x[1]*(2 - x[2]), [0, 1], 1000, 1020);
[0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0]

Orbk(2, x, x[1]*(2 - x[2]), [1, 1], 1000, 1020);
[1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1]


#Problem 4
#
plot(Dis1(z*(3 - z)*(5 - z), z, 0.01, 0.01, 20));
plot(Dis1(z*(3 - z)*(5 - z), z, 0.01, 3.01, 20));
plot(Dis1(z*(3 - z)*(5 - z), z, 0.01, 5.01, 20));
plot(Dis1(x^2*(3 - x)*(5 - x)*(7 - x), x, 0.01, 0.01, 20));
plot(Dis1(x^2*(3 - x)*(5 - x)*(7 - x), x, 0.01, 3.01, 20));
plot(Dis1(x^2*(3 - x)*(5 - x)*(7 - x), x, 0.01, 7.01, 20));