Moi j'aime bien... II

acidezabs a écrit:
[ironie]

j'adore résoudre des équations différentielles à 5h50 du matin

[/ironie]

// Programe solve equa(3) ; sujet 1 - question 4.1
// Détermination des coordonnées de l'état stationnaire pour nu=0.04

// d(alpha)/d(temps)=0 ==> nu-sigma*PHI=0 ==> PHI = nu/sigma
// d(gamma)/d(temps)=0 ==> gamma = (q*nu)/(h*k)

clear;

//définition de la fonction
deff('y=fct(alpha)','y=nu-sigma*(alpha*(1+alpha)/((L0/(1+g)^2)+(1+alpha)^2))');

//paramètres
L0=10^6;
q=100;
nu=0.04;
h=10;
k=0.4;
sigma=1.2;
g=(q*nu)/(h*k); // ==> g=1

// alpha0 est le point permettant d'initialiser l'algorithme de recherche
alpha0=rand(1);
sol=fsolve(alpha0,fct);

//affichage des cordonnées : alpha et gamma
printf('alpha = %f, gamma=%f',sol,g);

résultat obtenu :
--> alpha = 92.366351, gamma=1.000000

acidezabs a bossé sur:
// Programe solve equa(3) ; sujet 1 - question 4.1
// Détermination des coordonnées de l'état stationnaire pour nu=0.04

// d(alpha)/d(temps)=0 ==> nu-sigma*PHI=0 ==> PHI = nu/sigma
// d(gamma)/d(temps)=0 ==> gamma = (q*nu)/(h*k)

clear;

//définition de la fonction
deff('y=fct(alpha)','y=nu-sigma*(alpha*(1+alpha)/((L0/(1+g)^2)+(1+alpha)^2))');

//paramètres
L0=10^6;
q=100;
nu=0.04;
h=10;
k=0.4;
sigma=1.2;
g=(q*nu)/(h*k); // ==> g=1

// alpha0 est le point permettant d'initialiser l'algorithme de recherche
alpha0=rand(1);
sol=fsolve(alpha0,fct);

//affichage des cordonnées : alpha et gamma
printf('alpha = %f, gamma=%f',sol,g);

résultat obtenu :
--> alpha = 92.366351, gamma=1.000000

stargate a écrit:
mieux vaut tard que jamais comme ont dit !!

Marouz a écrit:
j'aime bien faire de longues parties de uno...