Le forum de XCAS

Publié : **jeu. nov. 13, 2014 9:08 pm**

Bonjour, dans le but d'étudier le modèle proie prédateur en TS, j'ai créer deux procédures qui calculent le nombres de lapins et de renards.

Code : Tout sélectionner

lapin(n):={
  local k, u, v, a, b, al, bl, L;
  u:=100;
  v:=100;
  a:=.001;
  b:=.045;
  al:=.0002;
  bl:=.025;
  L:=u;
  pour k de 1 jusque n faire
    u:=-a*u*v+(b+1)*u;
    v:=al*u*v+(1-bl)*v;
    L:=L,u;
  fpour;
  retourne L;
}:;

et

Code : Tout sélectionner

renard(n):={
  local k, u, v, a, b, al, bl, L;
  u:=100;
  v:=100;
  a:=.001;
  b:=.045;
  al:=.0002;
  bl:=.025;
  L:=v;
  pour k de 1 jusque n faire
    u:=-a*u*v+(b+1)*u;
    v:=al*u*v+(1-bl)*v;
    L:=L,v;
  fpour;
  retourne L;
}:;

Par contre je en sais pas quelle commande utiliser pour avoir les deux évolutions sur le même graphique mis à part après avoir ouvert une fenêtre graphique (nuage_points).

Je ne sais pas non plus représenter les renards en fonction des lapins à partir des deux listes crées.

Merci d'avance.

Bruno

Publié : **jeu. nov. 13, 2014 11:27 pm**

Bonjour, pour les 2 courbes sur le même graphique, j'ai trouvé dans la documentation :

Code : Tout sélectionner

display(nuage_points(lapin(500)), [rouge]) ; display(nuage_points(renard(500)), [vert])

Par contre, pour faire apparaitre le cycle, je sèche...

Bruno

Publié : **ven. nov. 14, 2014 6:00 am**

salut,
on peut creer une matrice avant d'appliquer nuage_points
nuage_points(tran([[lapin(500)],[renard(500)]])

Publié : **ven. nov. 14, 2014 12:47 pm**

Bonjour, merci pour la réponse.

je m'en étais tiré avec :

Code : Tout sélectionner

cycle(n):={
  local R,L,k;
  R:=renard(n);
  L:=lapin(n);
  pour k de 0 jusque n faire
    point((R[k], L[k]));
  fpour;
  retourne -1;
}:;

Merci encore

Bruno

Le forum de XCAS

Proie Prédateur discret - Affichage des graphqiues.

Proie Prédateur discret - Affichage des graphqiues.

Re: Proie Prédateur discret - Affichage des graphqiues.

Re: Proie Prédateur discret - Affichage des graphqiues.

Re: Proie Prédateur discret - Affichage des graphqiues.