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.
-----