Problème de division par zéro avec Scilab

[invite3751aec8]

Bonjour,

j'ai un problème avec mon programme Scilab.

je résous une équation différentielle et j'ai remarqué qua quand je passais un paramètre (K) en dessous de 1 Scilab me signalait une erreur ("erreur 27") à cause d'une soi disant division par zero. Pourtant ce K est au numérateur...


Voici le programme Scilab :

function [Xprime] = VanDerPol(t,X)

Xprime(1) = X(2)
Xprime(2) = g/R*(1-sin(X(1)))/(%pi/2-X(1))-(mu*(K)*R^2*((X(2))^2)/2)+/(R^2*mu*(%pi/2-X(1)))

endfunction

t=linspace(0,4,3000); X0=[85*%pi/180;0];
[u]=ode(X0,0,t,VanDerPol);

plot(t,u(1,*180/%pi,'r')
plot(t,R*u(2,)


Et voilà ce que la console Scilab me répond :

-->exec('C:\Documents and Settings\ademollerat\My Documents\accelerometre\Ecoule ment\VanDerPol-2.sci', -1)
Attention: redéfinition de la fonction: VanDerPol . Utilisez funcprot(0) pour éviter ce message

!--error 27
Division par zéro ...

at line 7 of function VanDerPol called by :
[u]=ode(X0,0,t,VanDerPol);
at line 35 of exec file called by :
exec('C:\Documents and Settings\ademollerat\My Documents\accelerometre\Ecoule ment\VanDerPol-2.sci', -1)


Est ce qu'il y a une erreur dans mon programme, sachant que cela fonctionne quand K>1 ?
Est ce que c'est un bug de Scilab ?
Quelqu'un a t-il déjà eu à faire à ce genre de problèmes précédemment ?

D'avance merci pour vos conseils et remarques.
-----

[invite3751aec8]

en fait c'est que quand K=1 qu'il n'y a pas d'erreurs ...

[maclag]

Peut-être que tu devrais essayer de décomposer la formule avec affichage des résultats intermédiaires pour voir un peu ce qui se passe?

S'il n'y a pas d'erreur avec les éléments décomposés, on les regroupe un à la fois jusqu'à ce que ça coince.