Projet scilab...

[invite42f97068]

Bonjour à tous, je vous écris car jusqu'à maintenant vous m'avez tous bien aidé et je vous en remercie. Alors je vous donne tout d'abord le code de tout mon programme:

clear
clf()

a1=4;//reproduction
b1=1;//mortalité
c1=2;//effort de pêche
a2=3;b2=1;c2=3;
k1=1;//migration de zone 1 vers 2
k2=2;//migration de zone 2 vers 1
//h=0.01;
N=1000;
dx=0.10;
h=dx/2

function z=f1(x,y)
 z=(a1-b1-c1-k1)*x+k2*y;
endfunction

function z=f2(x,y)
 z=(a2-b2-c2-k2)*y+k1*x;
endfunction



p1(1)=5;p2(1)=5;
for k=1:N
    p1(k+1)=p1(k)+h*f1(p1(k),p2(k));
    p2(k+1)=p2(k)+h*f2(p1(k),p2(k));
end;


xset("window",0)
plot2d(p1,p2)

u1=zeros(N+1,10); u2=u1;

for i=1:10
  u1(1,i)=p1(1)/10;
  u2(1,i)=p2(1)/10;
end;
for k =1:N
  u1(k,1)=a1;
  u2(k,1)=a2;
  p1(k+1)= p1(k)*dx;
  p2(k+1)= p2(k)*dx;
for i=2:10
    u1(k+1,i)= (u1(k,i)-(h/dx)*(u1(k,i)- u1(k,i-1)))/(1+h*(-b1-c1-k1));
    u2(k+1,i)= (u2(k,i)-(h/dx)*(u2(k,i)- u2(k,i-1)))/(1+h*(-b2-c2-k2));
    p1(k+1) = p1(k+1)+ u1(k+1,i)*dx;
    p2(k+1) = p2(k+1)+ u2(k+1,i)*dx;
end;
end;

xset("window",1)
plot2d(p1,p2)

Alors ça peut vous sembler bizarre mais notre prof n'a pas réussi à nous dire là ou ça bloque. Notre projet s'inspire de "Proie Prédateur".
Pour vous guider de la même manière qu'on l'a été, je vous dit que pour la premier graphique on est censé obtenir une spirale or comme vous le constaterez ce n'est pas du tout ce qu'on obtient et notre prof à faillit avoir une rupture d'anévrisme....Son aide nous embrouille plus qu'autre chose mais malheureusement nous avons besoin d'aide donc j'espère qu'ici je trouverai ce qui me mande.

Dans les deux graphique, pour le premier j'utilise euler pour obtenir la courbe de population. et dans le deuxième la même chose en fonction de la taille de la population.

Merci pour votre aide.
-----

[invite551c2897]

Bonjour.
Peut-être l'erreur :

function z=f2(x,y)
 z=(a2-b2-c2-k2)*x+k1*y;
endfunction