<html>
<head>
<style><!--
.hmmessage P
{
margin:0px;
padding:0px
}
body.hmmessage
{
font-size: 10pt;
font-family:Tahoma
}
--></style>
</head>
<body class='hmmessage'><div dir='ltr'>
<font class="Apple-style-span" face="Tahoma" size="2">Hi, </font><div style="font-family: Tahoma; font-size: 10pt; "><br></div><div style="font-family: Tahoma; font-size: 10pt; ">I am using scilab to solve a differential equation and I have an error ("error 27") I don't understand.</div><div style="font-family: Tahoma; font-size: 10pt; "><br></div><div style="font-family: Tahoma; font-size: 10pt; ">When K=1 pr K>1 it's okay but as soon as K<1 there is the error.</div><div style="font-family: Tahoma; font-size: 10pt; "><br></div><div style="font-family: Tahoma; font-size: 10pt; ">This is my code : </div><div style="font-family: Tahoma; font-size: 10pt; "><br></div><div><pre style="font-family: Monospaced; font-style: normal; font-size: 14px; "><span style="color:rgb(176,24,19);">function</span> <span style="color:rgb(74,85,219);">[</span><span style="color:rgb(131,67,16);font-weight:bold;">Xprime</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(0,0,0);text-decoration:underline;">VanDerPol</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">t</span><span style="color:rgb(0,0,0);">, </span><span style="color:rgb(131,67,16);font-weight:bold;">X</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(131,67,16);font-weight:bold;">Xprime</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(131,67,16);font-weight:bold;">X</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(131,67,16);font-weight:bold;">Xprime</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span> <span style="color:rgb(92,92,92);">=</span> <span style="color:rgb(0,0,0);">g</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(0,0,0);">R</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(50,185,185);">sin</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">X</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(218,112,214);">%pi</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(131,67,16);font-weight:bold;">X</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">mu</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">K</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">R</span><span style="color:rgb(92,92,92);">^</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(131,67,16);font-weight:bold;">X</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">^</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">+</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">R</span><span style="color:rgb(92,92,92);">^</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">mu</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(218,112,214);">%pi</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(92,92,92);">-</span><span style="color:rgb(131,67,16);font-weight:bold;">X</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(176,24,19);">endfunction</span>
<span style="color:rgb(0,0,0);">t</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(174,92,176);text-decoration:underline;">linspace</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">0</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">4</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">3000</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span> <span style="color:rgb(0,0,0);">X0</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(74,85,219);">[</span><span style="color:rgb(188,143,143);">85</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(218,112,214);">%pi</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(188,143,143);">180</span><span style="color:rgb(0,0,0);">;</span><span style="color:rgb(188,143,143);">0</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(74,85,219);">[</span><span style="color:rgb(0,0,0);">u</span><span style="color:rgb(74,85,219);">]</span><span style="color:rgb(92,92,92);">=</span><span style="color:rgb(50,185,185);">ode</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">X0</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">0</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">t</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);text-decoration:underline;">VanDerPol</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(0,0,0);">;</span>
<span style="color:rgb(174,92,176);text-decoration:underline;">plot</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">t</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">u</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">1</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(188,143,143);">180</span><span style="color:rgb(92,92,92);">/</span><span style="color:rgb(218,112,214);">%pi</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(188,143,143);">r</span><span style="color:rgb(188,143,143);">'</span><span style="color:rgb(74,85,219);">)</span>
<span style="color:rgb(174,92,176);text-decoration:underline;">plot</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(0,0,0);">t</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(0,0,0);">R</span><span style="color:rgb(92,92,92);">*</span><span style="color:rgb(0,0,0);">u</span><span style="color:rgb(74,85,219);">(</span><span style="color:rgb(188,143,143);">2</span><span style="color:rgb(0,0,0);">,</span><span style="color:rgb(255,170,0);">:</span><span style="color:rgb(74,85,219);">)</span><span style="color:rgb(74,85,219);">)</span></pre><pre style="font-family: Monospaced; font-style: normal; font-size: 14px; "><span style="color:rgb(74,85,219);"><br></span></pre><pre style="font-family: Monospaced; font-style: normal; font-size: 14px; "><span style="color:rgb(74,85,219);"><br></span></pre><pre style="font-family: Monospaced; font-style: normal; font-size: 14px; "><span style="color:rgb(74,85,219);">And there is the Scilab message when K<1 : </span></pre><pre style="font-family: Monospaced; font-style: normal; font-size: 14px; "><span style="color:rgb(74,85,219);"><br></span></pre><pre><font class="Apple-style-span" color="#4a55db" face="Monospaced"><span class="Apple-style-span" style="font-size: 14px;">-->exec('C:\Documents and Settings\ademollerat\My Documents\accelerometre\Ecoulement\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\Ecoulement\VanDerPol-2.sci', -1)</span></font></pre><pre><font class="Apple-style-span" color="#4a55db" face="Monospaced"><span class="Apple-style-span" style="font-size: 14px;"><br></span></font></pre><pre><font class="Apple-style-span" color="#4a55db" face="Monospaced"><span class="Apple-style-span" style="font-size: 14px;"><br></span></font></pre><pre><font class="Apple-style-span" color="#4a55db" face="Monospaced"><span class="Apple-style-span" style="font-size: 14px;">Maybe I have made a mistake or maybe it is a bug.</span></font></pre><pre><font class="Apple-style-span" color="#4a55db" face="Monospaced"><span class="Apple-style-span" style="font-size: 14px;"><br></span></font></pre><pre><font class="Apple-style-span" color="#4a55db" face="Monospaced"><span class="Apple-style-span" style="font-size: 14px;">Thank you for helping me.</span></font></pre><pre><font class="Apple-style-span" color="#4a55db" face="Monospaced"><span class="Apple-style-span" style="font-size: 14px;"><br></span></font></pre><pre><font class="Apple-style-span" color="#4a55db" face="Monospaced"><span class="Apple-style-span" style="font-size: 14px;">Best regards</span></font></pre></div> </div></body>
</html>