share/dynamics/rtest_dynamics.mac

   1 (kill(all),0);
   2 0$
   3
   4 last(rk(-n,n,1,[t,0,0.3,0.1]));
   5 [0.3, .7408184220011778]$
   6
   7 last(rk(t-x^2,x,1,[t,0,8,0.1]));
   8 [8.0, 2.796249427928215]$
   9
  10 length(rk(t-x^2,x,1,[t,0,8,0.1]));
  11 81$
  12
  13 last(rk(t-x^2,x,1,[t,2,3,0.1]));
  14 [3.0, 1.618275167116153]$
  15
  16 (evolution(cos(y), 2, 11),0);
  17 0;
  18
  19 (staircase(cos(y), 1, 11, [y, 0, 1.2]),0);
  20 0;
  21
  22 (orbits(x^2+a, 0, 50, 100, [a, -2, 0.25], [style, dots]),0);
  23 0;
  24
  25 (
  26 f: 0.6*x*(1+2*x)+0.8*y*(x-1)-y^2-0.9,
  27 g: 0.1*x*(1-6*x+4*y)+0.1*y*(1+9*y)-0.4,
  28 evolution2d([f,g], [x,y], [-0.5,0], 5000, [style,dots]),0
  29 );
  30 0;
  31
  32 (chaosgame([[0, 0], [1, 0], [0.5, sqrt(3)/2]], [0.1, 0.1], 1/2,
  33                  3000, [style, dots]),0);
  34 0;
  35
  36 (
  37 a1: matrix([0.85,0.04],[-0.04,0.85]),
  38 a2: matrix([0.2,-0.26],[0.23,0.22]),
  39 a3: matrix([-0.15,0.28],[0.26,0.24]),
  40 a4: matrix([0,0],[0,0.16]),
  41 p1: [0,1.6],
  42 p2: [0,1.6],
  43 p3: [0,0.44],
  44 p4: [0,0],
  45 w: [85,92,99,100],
  46 ifs(w, [a1,a2,a3,a4], [p1,p2,p3,p4], [5,0], 5000, [style,dots]),0
  47 );
  48 0;
  49
  50 set_plot_option([plot_format,gnuplot]);
  51
  52 (evolution(cos(y), 2, 11),0);
  53 0;
  54
  55 (staircase(cos(y), 1, 11, [y, 0, 1.2]),0);
  56 0;
  57
  58 (orbits(x^2+a, 0, 50, 100, [a, -2, 0.25], [style, dots]),0);
  59 0;
  60
  61 (
  62 f: 0.6*x*(1+2*x)+0.8*y*(x-1)-y^2-0.9,
  63 g: 0.1*x*(1-6*x+4*y)+0.1*y*(1+9*y)-0.4,
  64 evolution2d([f,g], [x,y], [-0.5,0], 5000, [style,dots]),0
  65 );
  66 0;
  67
  68 (chaosgame([[0, 0], [1, 0], [0.5, sqrt(3)/2]], [0.1, 0.1], 1/2,
  69                  3000, [style, dots]),0);
  70 0;
  71
  72 (
  73 a1: matrix([0.85,0.04],[-0.04,0.85]),
  74 a2: matrix([0.2,-0.26],[0.23,0.22]),
  75 a3: matrix([-0.15,0.28],[0.26,0.24]),
  76 a4: matrix([0,0],[0,0.16]),
  77 p1: [0,1.6],
  78 p2: [0,1.6],
  79 p3: [0,0.44],
  80 p4: [0,0],
  81 w: [85,92,99,100],
  82 ifs(w, [a1,a2,a3,a4], [p1,p2,p3,p4], [5,0], 5000, [style,dots]),0
  83 );
  84 0;
  85
  86 set_plot_option([plot_format,gnuplot_pipes]);
  87
  88 (reset(numer),0);
  89 0$