% TeXgraph version 1.93
% Fenetre Xmin Xmax Ymin Ymax Xscale Yscale
100#0.013#0.0215#-8.91610081502822#-0.1695761651398#2235.29411764705#2.1722913683485##
% Marges gauche droite haut bas cadre gestion_couleur comptgraph
101#0.5#0.5#0.5#0.5#1#1#8##
% Affectation des Variables theta et phi
18##[theta:=0.5236, phi:=1.0472, OriginalCoord(1)]##
% Déclaration des Variables Globales
15#A#1.1##
15#omega#2*pi/300##
15#gamma#0.002083##
15#m#0.0334##
15#n#100##
15#x0#0.25##
15#y0#1##
15#N#50 {points / intervalle}##
% Déclaration des Macros
16#f#[$t:=%1, $x:=%2, $y:=%3,
{ equadif du type: x'(t) + i*y'(t) }
x*(A*y-1-m*sin(omega*t)) +i*gamma*(1-y*(1+x))
]
##
% Déclaration des Eléments graphiques
% Tracer (Utilisateur)
18##[AutoReCalc:=0,tMin:=-5,tMax:=5,LabelSize:=1]##
14#Tracer#[$T:=2*pi/omega,tMin:=0, tMax:=n*T, NbPoints:=N*n+1,

$S:=for m from 0.013 to 0.022 step 0.0005
  do
   stock:=Get(EquaDif( f(t,x,y), 0, x0+i*y0),0),
   Del(stock, 1, N*50),
   k:=0,
   for $z in stock do
      Inc(k,1),
      if k=1 then m+i*ln(Re(z)) elif k=N then k:=0 fi
   od
  od,
R:=rectangle(S),
Fenetre(R[1],R[2]), size(20,0),
Point(S)
]
#-1##