Essencial function for working with dynamic system

3D_plot: cobweb() = return cobweb plot orbit_diagram() = return orbit plot Lyap() = return Lyapunov exponent plot

Bisection_method: bisection_2_order() = return root of second order polynomial bisection_3_order() = return root of third order polynomial eigvalues_bisection_method() = return eigvalues of 3x3 matrix

Dynamic: Solution_discrete_system() = solve the dynamic plot_all_X() = plot all variables