<html>
<head>
<meta content="text/html; charset=ISO-8859-1"
http-equiv="Content-Type">
</head>
<body bgcolor="#FFFFFF" text="#000000">
<div class="moz-cite-prefix">Le 13/11/2012 08:47, Paul Carrico a
écrit :<br>
</div>
<blockquote
cite="mid:000301cdc173$2e80df90$8b829eb0$@carrico@free.fr"
type="cite">
<meta http-equiv="Content-Type" content="text/html;
charset=ISO-8859-1">
<meta name="Generator" content="Microsoft Word 12 (filtered
medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Wingdings;
panose-1:5 0 0 0 0 0 0 0 0 0;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Monospaced;
panose-1:0 0 0 0 0 0 0 0 0 0;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
pre
{mso-style-priority:99;
mso-style-link:"Préformaté HTML Car";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
{mso-style-priority:34;
margin-top:0cm;
margin-right:0cm;
margin-bottom:0cm;
margin-left:36.0pt;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
span.EmailStyle17
{mso-style-type:personal-compose;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.PrformatHTMLCar
{mso-style-name:"Préformaté HTML Car";
mso-style-priority:99;
mso-style-link:"Préformaté HTML";
font-family:"Courier New";}
.MsoChpDefault
{mso-style-type:export-only;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
/* List Definitions */
@list l0
{mso-list-id:1275937926;
mso-list-type:hybrid;
mso-list-template-ids:1302116462 -248580400 67895299 67895301 67895297 67895299 67895301 67895297 67895299 67895301;}
@list l0:level1
{mso-level-start-at:0;
mso-level-number-format:bullet;
mso-level-text:-;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;
font-family:"Calibri","sans-serif";
mso-fareast-font-family:Calibri;}
ol
{margin-bottom:0cm;}
ul
{margin-bottom:0cm;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div class="WordSection1"><tt>.../...</tt><tt><span
style="font-size: 9pt;" lang="EN-US"><o:p> <br>
</o:p></span></tt>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
black;" lang="EN-US">Epsilon</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">=</span> <span
style="color:rosybrown">10</span><span
style="color:#5C5C5C">^</span><span
style="color:#4A55DB">(</span><span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">8</span><span
style="color:#4A55DB">)</span> <span
style="color:#5C5C5C">*</span><span
style="color:#4A55DB">[</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rosybrown;" lang="EN-US">4.441D-08</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:rosybrown">0.</span> <span
style="color:rosybrown">0.</span> <span
style="color:rosybrown">0.</span> <span
style="color:rosybrown">0.</span> <span
style="color:rosybrown">0.</span> <span
style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(92, 92, 92);" lang="EN-US">-</span></tt><tt><span
style="font-size: 9pt; color: rosybrown;" lang="EN-US">1022.</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">6278.</span> <span
style="color:rosybrown">6914.</span> <span
style="color:rosybrown">210.6</span> <span
style="color:rosybrown">783100.</span> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">33910.</span> <span
style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(92, 92, 92);" lang="EN-US">-</span></tt><tt><span
style="font-size: 9pt; color: rosybrown;" lang="EN-US">2316.</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">24190.</span> <span
style="color:rosybrown">25930.</span> <span
style="color:rosybrown">894.6</span> <span
style="color:rosybrown">1545000.</span> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">67320.</span> <span
style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(92, 92, 92);" lang="EN-US">-</span></tt><tt><span
style="font-size: 9pt; color: rosybrown;" lang="EN-US">3608.</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">52660.</span> <span
style="color:rosybrown">55610.</span> <span
style="color:rosybrown">2004.</span> <span
style="color:rosybrown">2284000.</span> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">99750.</span> <span
style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(92, 92, 92);" lang="EN-US">-</span></tt><tt><span
style="font-size: 9pt; color: rosybrown;" lang="EN-US">4895.</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">90810.</span> <span
style="color:rosybrown">94970.</span> <span
style="color:rosybrown">3508.</span> <span
style="color:rosybrown">3002000.</span> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">131300.</span> <span
style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(92, 92, 92);" lang="EN-US">-</span></tt><tt><span
style="font-size: 9pt; color: rosybrown;" lang="EN-US">6171.</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">137700.</span> <span
style="color:rosybrown">143100.</span> <span
style="color:rosybrown">5365.</span> <span
style="color:rosybrown">3699000.</span> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">161900.</span> <span
style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(92, 92, 92);" lang="EN-US">-</span></tt><tt><span
style="font-size: 9pt; color: rosybrown;" lang="EN-US">7429.</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">192400.</span> <span
style="color:rosybrown">199000.</span> <span
style="color:rosybrown">7533.</span> <span
style="color:rosybrown">4373000.</span> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">191600.</span> <span
style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(92, 92, 92);" lang="EN-US">-</span></tt><tt><span
style="font-size: 9pt; color: rosybrown;" lang="EN-US">8675.</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">254200.</span> <span
style="color:rosybrown">262000.</span> <span
style="color:rosybrown">9997.</span> <span
style="color:rosybrown">5027000.</span> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">220400.</span> <span
style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(92, 92, 92);" lang="EN-US">-</span></tt><tt><span
style="font-size: 9pt; color: rosybrown;" lang="EN-US">8675.</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">254200.</span> <span
style="color:rosybrown">262000.</span> <span
style="color:rosybrown">9997.</span> <span
style="color:rosybrown">5027000.</span> <span
style="color:#5C5C5C">-</span><span
style="color:rosybrown">220400.</span> <o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(74, 85, 219);" lang="EN-US">]</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"><o:p> </o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(74, 85, 219);" lang="EN-US">[</span></tt><tt><span
style="font-size: 9pt; color: black;" lang="EN-US">p,c</span></tt><tt><span
style="font-size: 9pt; color: rgb(74, 85, 219);"
lang="EN-US">]</span></tt><tt><span style="font-size:
9pt;" lang="EN-US"> <span style="color:#5C5C5C">=</span>
<span style="color:#32B9B9">size</span><span
style="color:#4A55DB">(</span><span style="color:black">Epsilon</span><span
style="color:#4A55DB">)</span><span style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><i><span style="font-size: 9pt; color:
rgb(100, 174, 100);" lang="EN-US">// number of time
steps</span></i></tt><tt><span style="font-size: 9pt;"
lang="EN-US"><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
black;" lang="EN-US">step_</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">=</span> <span
style="color:rosybrown">2</span><span
style="color:black">;</span> <i><span
style="color:#64AE64">//</span></i><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"><o:p> </o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
black;" lang="EN-US">nb_loops</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">=</span> <span
style="color:rosybrown">0</span><span
style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><u><span style="font-size: 9pt; color:
rgb(174, 92, 176);" lang="EN-US">tic</span></u></tt><tt><span
style="font-size: 9pt; color: rgb(74, 85, 219);"
lang="EN-US">()</span></tt><tt><span style="font-size:
9pt;" lang="EN-US"><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(160, 32, 240);" lang="EN-US">for</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:black">i</span> <span
style="color:#5C5C5C">=</span> <span
style="color:rosybrown">1</span> <span
style="color:#FFAA00">:</span> <span
style="color:black">p</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> <span style="color:#A020F0">for</span> <span
style="color:black">theta</span> <span
style="color:#5C5C5C">=</span> <span
style="color:rosybrown">0</span> <span
style="color:#FFAA00">:</span> <span
style="color:black">step_</span> <span
style="color:#FFAA00">:</span> <span
style="color:rosybrown">180</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> <span style="color:#A020F0">for</span>
<span style="color:black">fi</span> <span
style="color:#5C5C5C">=</span> <span
style="color:rosybrown">0</span> <span
style="color:#FFAA00">:</span> <span
style="color:black">step_</span> <span
style="color:#FFAA00">:</span> <span
style="color:rosybrown">180</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> </span></tt><tt><span
style="font-size: 9pt; color: black;">Epsilon_matrix</span></tt><tt><span
style="font-size: 9pt;"> <span style="color:#5C5C5C">=</span>
<u><span style="color:#AE5CB0">fct_vector2matrix</span></u><span
style="color:#4A55DB">(</span><span style="color:black">Epsilon</span><span
style="color:#4A55DB">(</span><span style="color:black">i,</span><span
style="color:#FFAA00">:</span><span
style="color:#4A55DB">))</span><span style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;">
<span style="color:black">Epsilon_prime_matrix</span> <span
style="color:#5C5C5C">=</span> <u><span
style="color:#AE5CB0">fct_transformation_M</span></u><span
style="color:#4A55DB">(</span><span style="color:black">Epsilon_matrix,theta,fi</span><span
style="color:#4A55DB">)</span><span style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;">
</span></tt><tt><span style="font-size: 9pt; color: black;"
lang="EN-US">Epsilon_prime</span></tt><tt><span
style="font-size: 9pt; color: rgb(74, 85, 219);"
lang="EN-US">(</span></tt><tt><span style="font-size: 9pt;
color: black;" lang="EN-US">i,</span></tt><tt><span
style="font-size: 9pt; color: rgb(255, 170, 0);"
lang="EN-US">:</span></tt><tt><span style="font-size: 9pt;
color: rgb(74, 85, 219);" lang="EN-US">)</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">=</span> <u><span
style="color:#AE5CB0">fct_matrix2vector</span></u><span
style="color:#4A55DB">(</span><span style="color:black">Epsilon_prime_matrix</span><span
style="color:#4A55DB">)</span><span style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> <span style="color:black">nb_loops</span>
<span style="color:#5C5C5C">=</span> <span
style="color:black">nb_loops</span> <span
style="color:#5C5C5C">+</span> <span
style="color:rosybrown">1</span><span
style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> <span style="color:#A020F0">end</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> <span style="color:#A020F0">end</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(160, 32, 240);" lang="EN-US">end</span></tt></p>
</div>
</blockquote>
Since <tt><span style="font-size: 9pt;"><span style="color:black"><tt><span
style="font-size: 9pt; color: black;">Epsilon_matrix </span></tt></span><span
style="color:#4A55DB"></span><span style="color:black"></span><span
style="color:#FFAA00"></span><span style="color:#4A55DB">does
not depend neither on theta nor on fi, </span></span></tt>it
can be <br>
set only once when opening the for { i} loop:<br>
<br>
<div class="WordSection1">
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(160, 32, 240);" lang="EN-US">for</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:black">i</span> <span style="color:#5C5C5C">=</span>
<span style="color:rosybrown">1</span> <span
style="color:#FFAA00">:</span> <span style="color:black">p</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> </span></tt><tt><span style="font-size:
9pt; color: black;">Epsilon_matrix</span></tt><tt><span
style="font-size: 9pt;"> <span style="color:#5C5C5C">=</span>
<u><span style="color:#AE5CB0">fct_vector2matrix</span></u><span
style="color:#4A55DB">(</span><span style="color:black">Epsilon</span><span
style="color:#4A55DB">(</span><span style="color:black">i,</span><span
style="color:#FFAA00">:</span><span style="color:#4A55DB">))</span><span
style="color:black">;</span><o:p></o:p></span></tt></p>
<tt><span style="font-size: 9pt;"></span></tt>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> <span style="color:#A020F0">for</span> <span
style="color:black">theta</span> <span
style="color:#5C5C5C">=</span> <span
style="color:rosybrown">0</span> <span
style="color:#FFAA00">:</span> <span style="color:black">step_</span>
<span style="color:#FFAA00">:</span> <span
style="color:rosybrown">180</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> <span style="color:#A020F0">for</span>
<span style="color:black">fi</span> <span
style="color:#5C5C5C">=</span> <span
style="color:rosybrown">0</span> <span
style="color:#FFAA00">:</span> <span style="color:black">step_</span>
<span style="color:#FFAA00">:</span> <span
style="color:rosybrown">180</span><o:p></o:p></span></tt></p>
<tt><span style="font-size: 9pt;"> <span
style="color:black">Epsilon_prime_matrix</span> <span
style="color:#5C5C5C">=</span> <u><span
style="color:#AE5CB0">fct_transformation_M</span></u><span
style="color:#4A55DB">(</span><span style="color:black">Epsilon_matrix,theta,fi</span><span
style="color:#4A55DB">)</span><span style="color:black">;</span><o:p></o:p></span></tt>
<p class="MsoNormal"><tt><span style="font-size: 9pt;">
</span></tt><tt><span style="font-size: 9pt; color: black;"
lang="EN-US">Epsilon_prime</span></tt><tt><span
style="font-size: 9pt; color: rgb(74, 85, 219);"
lang="EN-US">(</span></tt><tt><span style="font-size: 9pt;
color: black;" lang="EN-US">i,</span></tt><tt><span
style="font-size: 9pt; color: rgb(255, 170, 0);"
lang="EN-US">:</span></tt><tt><span style="font-size: 9pt;
color: rgb(74, 85, 219);" lang="EN-US">)</span></tt><tt><span
style="font-size: 9pt;" lang="EN-US"> <span
style="color:#5C5C5C">=</span> <u><span
style="color:#AE5CB0">fct_matrix2vector</span></u><span
style="color:#4A55DB">(</span><span style="color:black">Epsilon_prime_matrix</span><span
style="color:#4A55DB">)</span><span style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> <span style="color:black">nb_loops</span>
<span style="color:#5C5C5C">=</span> <span
style="color:black">nb_loops</span> <span
style="color:#5C5C5C">+</span> <span
style="color:rosybrown">1</span><span style="color:black">;</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> <span style="color:#A020F0">end</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt;"
lang="EN-US"> <span style="color:#A020F0">end</span><o:p></o:p></span></tt></p>
<p class="MsoNormal"><tt><span style="font-size: 9pt; color:
rgb(160, 32, 240);" lang="EN-US">end</span></tt></p>
</div>
<br>
<br>
Regards<br>
Samuel<br>
<br>
</body>
</html>