Problème de division par zéro avec Scilab
Répondre à la discussion
Affichage des résultats 1 à 3 sur 3

Problème de division par zéro avec Scilab



  1. #1
    invite3751aec8

    Problème de division par zéro avec Scilab


    ------

    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.

    -----

  2. #2
    invite3751aec8

    Re : Problème de division par zéro avec Scilab

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

  3. #3
    invitef433700c

    Re : Problème de division par zéro avec Scilab

    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.

Discussions similaires

  1. Probleme avec fscanfMat sur scilab
    Par invite3a97487a dans le forum Programmation et langages, Algorithmique
    Réponses: 0
    Dernier message: 27/03/2011, 11h17
  2. Problème avec Scilab
    Par invite2f949375 dans le forum Logiciel - Software - Open Source
    Réponses: 2
    Dernier message: 21/10/2008, 23h18
  3. Problème de résolution d'équat-diff avec Scilab
    Par invitec12155c1 dans le forum Mathématiques du supérieur
    Réponses: 0
    Dernier message: 26/09/2008, 10h13