<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=iso-8859-1">
<meta name="Generator" content="Microsoft Word 14 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Tahoma;
panose-1:2 11 6 4 3 5 4 4 2 4;}
@font-face
{font-family:Consolas;
panose-1:2 11 6 9 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;}
p
{mso-style-priority:99;
mso-margin-top-alt:auto;
margin-right:0cm;
mso-margin-bottom-alt:auto;
margin-left:0cm;
font-size:12.0pt;
font-family:"Times New Roman","serif";}
pre
{mso-style-priority:99;
mso-style-link:"Préformaté HTML Car";
margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
p.MsoAcetate, li.MsoAcetate, div.MsoAcetate
{mso-style-priority:99;
mso-style-link:"Texte de bulles Car";
margin:0cm;
margin-bottom:.0001pt;
font-size:8.0pt;
font-family:"Tahoma","sans-serif";}
span.PrformatHTMLCar
{mso-style-name:"Préformaté HTML Car";
mso-style-priority:99;
mso-style-link:"Préformaté HTML";
font-family:"Courier New";
mso-fareast-language:FR;}
p.HTMLPreformatted, li.HTMLPreformatted, div.HTMLPreformatted
{mso-style-name:"HTML Preformatted";
mso-style-link:"HTML Preformatted Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:11.0pt;
font-family:"Calibri","sans-serif";}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:Consolas;}
span.EmailStyle22
{mso-style-type:personal;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.EmailStyle23
{mso-style-type:personal;
font-family:"Calibri","sans-serif";
color:black;
font-weight:normal;
font-style:normal;
text-decoration:none none;}
span.TextedebullesCar
{mso-style-name:"Texte de bulles Car";
mso-style-priority:99;
mso-style-link:"Texte de bulles";
font-family:"Tahoma","sans-serif";}
span.EmailStyle27
{mso-style-type:personal-reply;
font-family:"Calibri","sans-serif";
color:#1F497D;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></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]-->
</head>
<body lang="FR" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D">Thanks Rafael for the answer<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D">One explanation of the error is that I used a and b letter instead of aa and bb in your example; your code, if I change aa -> a and bb -> it fails again … why ? (I do not understand)<o:p></o:p></span></p>
<p class="MsoNormal"><o:p></o:p></p>
<p><b><i><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:gray">EXPORT CONTROL :
<br>
</span></i></b><b><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:black">Cet email ne contient pas de données techniques<br>
This email does not contain technical data</span></b><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><a name="_____replyseparator"></a><b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif"">De :</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif""> users [mailto:users-bounces@lists.scilab.org]
<b>De la part de</b> Rafael Guerra<br>
<b>Envoyé :</b> jeudi 26 mai 2016 12:26<br>
<b>À :</b> 'Users mailing list for Scilab'<br>
<b>Objet :</b> [EXTERNAL] Re: [Scilab-users] Vectorization issue<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black">Hi Paul,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black">It seems that there may be some issues with the .*, .^ syntax and also the use of the same variable in different context (b).<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black">Editing the vectorization part of your code as follows:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black"><o:p> </o:p></span></p>
<p class="MsoNormal"><i><span lang="EN-US" style="font-size:12.0pt;font-family:"Monospaced","serif";color:#64AE64">// using vectorization</span></i><span lang="EN-US" style="font-size:12.0pt;font-family:"Monospaced","serif""><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:12.0pt;font-family:"Monospaced","serif";color:black">aa</span><span lang="EN-US" style="font-size:12.0pt;font-family:"Monospaced","serif"">
<span style="color:#5C5C5C">=</span> <span style="color:#4A55DB">[</span><span style="color:rosybrown">0</span><span style="color:#FFAA00">:</span><span style="color:black">n</span><span style="color:#4A55DB">]</span><span style="color:#5C5C5C">'</span><span style="color:black">;</span><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif";color:black">bb</span><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif"">
<span style="color:#5C5C5C">=</span> <span style="color:#32B9B9">ones</span><span style="color:#4A55DB">(</span><span style="color:black">n,</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:black">;</span><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif";color:black">i</span><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif"">
<span style="color:#5C5C5C">=</span> <span style="color:black">aa</span><span style="color:#5C5C5C">.*.</span><span style="color:black">bb;</span><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif";color:black">j</span><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif"">
<span style="color:#5C5C5C">=</span> <span style="color:black">bb</span><span style="color:#5C5C5C">.*.</span><span style="color:black">aa;</span><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif""><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif";color:black">Km2</span><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif"">
<span style="color:#5C5C5C">=</span> <span style="color:#32B9B9">zeros</span><span style="color:#4A55DB">(</span><span style="color:black">n</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">(</span><span style="color:black">n</span><span style="color:#5C5C5C">+</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:black">,</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:black">;</span><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif";color:black">Km2</span><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif"">
<span style="color:#5C5C5C">=</span> <span style="color:orchid">%pi</span><span style="color:#5C5C5C">^</span><span style="color:rosybrown">6</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">((</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:black">i</span><span style="color:#5C5C5C">+</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">2</span>
<span style="color:#5C5C5C">.*</span> <span style="color:black">j</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">2</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">.*</span><span style="color:#4A55DB">((</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:black">i</span><span style="color:#5C5C5C">+</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">2.</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">(</span><span style="color:black">b</span><span style="color:#5C5C5C">/</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:black">d</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">^</span><span style="color:rosybrown">2</span>
<span style="color:#5C5C5C">+</span> <span style="color:black">j</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">2</span><span style="color:#4A55DB">)</span>
<span style="color:black">;</span><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif";color:black">Km2</span><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif"">
<span style="color:#5C5C5C">=</span> <span style="color:#4A55DB">(</span><span style="color:rosybrown">144</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">(</span><span style="color:black">b</span><span style="color:#5C5C5C">/</span><span style="color:black">d</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">4</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">/</span>
<span style="color:#32B9B9">sum</span><span style="color:#4A55DB">(</span><span style="color:black">Km2</span><span style="color:#4A55DB">)</span><span style="color:black">;</span><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif";color:black">Km2</span><span lang="PT" style="font-size:12.0pt;font-family:"Monospaced","serif"">
<span style="color:#5C5C5C">=</span> <span style="color:rosybrown">1.2</span> <span style="color:#5C5C5C">
+</span> <span style="color:#4A55DB">(</span><span style="color:black">nu</span><span style="color:#5C5C5C">/</span><span style="color:#4A55DB">(</span><span style="color:rosybrown">1</span><span style="color:#5C5C5C">+</span><span style="color:black">nu</span><span style="color:#4A55DB">))</span><span style="color:#5C5C5C">*</span><span style="color:black">Km2;</span><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:12.0pt;font-family:"Monospaced","serif";color:black">k2</span><span lang="EN-US" style="font-size:12.0pt;font-family:"Monospaced","serif"">
<span style="color:#5C5C5C">=</span> <span style="color:rosybrown">1</span><span style="color:#5C5C5C">/</span><span style="color:black">Km2</span><o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black">produces:<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black"> k = 0.8333090
<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black"> k2 = 0.8333333 <o:p>
</o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black">A small difference. I did not check the formulas but it might be numerical error only.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black">Regards,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black">Rafael<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:black"><o:p> </o:p></span></p>
<div>
<div style="border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span lang="EN-US">From:</span></b><span lang="EN-US"> users [<a href="mailto:users-bounces@lists.scilab.org">mailto:users-bounces@lists.scilab.org</a>]
<b>On Behalf Of </b>Carrico, Paul<br>
<b>Sent:</b> Thursday, May 26, 2016 11:30 AM<br>
<b>To:</b> International users mailing list for Scilab. (<a href="mailto:users@lists.scilab.org">users@lists.scilab.org</a>) <<a href="mailto:users@lists.scilab.org">users@lists.scilab.org</a>><br>
<b>Subject:</b> [Scilab-users] Vectorization issue<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">Dear</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">I failed in using vectorization in the example immediately bellow; but I do not remember if it’s possible : any advice ?</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">Thanks</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">Paul</span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US">######################################################################################</span><o:p></o:p></p>
<pre><span style="font-family:"Monospaced","serif";color:#32B9B9">mode</span><span style="font-family:"Monospaced","serif";color:#4A55DB">(</span><span style="font-family:"Monospaced","serif";color:rosybrown">0</span><span style="font-family:"Monospaced","serif";color:#4A55DB">)</span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span style="font-family:"Monospaced","serif""> </span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">n</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">20</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">Km</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">0</span><span style="color:black">;</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif""> </span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">b</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">100</span><span style="color:black">;</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">d</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">20</span><span style="color:black">;</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">nu</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">0.3</span><span style="color:black">;</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif""> </span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><i><span lang="EN-US" style="font-family:"Monospaced","serif";color:#64AE64">// using loops</span></i><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:#A020F0">for</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:black">i</span> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">0</span> <span style="color:#FFAA00">:</span> <span style="color:black">n</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#A020F0">for</span> <span style="color:black">j</span> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">1</span> <span style="color:#FFAA00">:</span> <span style="color:black">n</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:black">Km</span> <span style="color:#5C5C5C">=</span> <span style="color:black">Km</span> <span style="color:#5C5C5C">+</span> <span style="color:#4A55DB">(</span><span style="color:rosybrown">144</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">(</span><span style="color:black">b</span><span style="color:#5C5C5C">/</span><span style="color:black">d</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">^</span><span style="color:rosybrown">4</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">/</span><span style="color:#4A55DB">(</span> <span style="color:orchid">%pi</span><span style="color:#5C5C5C">^</span><span style="color:rosybrown">6</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">(</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:black">i</span><span style="color:#5C5C5C">+</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">^</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:black">j</span><span style="color:#5C5C5C">^</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">((</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:black">i</span><span style="color:#5C5C5C">+</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">^</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">(</span><span style="color:black">b</span><span style="color:#5C5C5C">/</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:black">d</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">^</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">+</span><span style="color:black">j</span><span style="color:#5C5C5C">^</span><span style="color:rosybrown">2</span><span style="color:#4A55DB">)</span> <span style="color:#4A55DB">)</span><span style="color:black">;</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#A020F0">end</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:#A020F0">end</span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif""> </span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">Km</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">1.2</span> <span style="color:#5C5C5C">+</span> <span style="color:#4A55DB">(</span><span style="color:black">nu</span><span style="color:#5C5C5C">/</span><span style="color:#4A55DB">(</span><span style="color:rosybrown">1</span><span style="color:#5C5C5C">+</span><span style="color:black">nu</span><span style="color:#4A55DB">))</span><span style="color:#5C5C5C">*</span><span style="color:black">Km</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">k</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">1</span><span style="color:#5C5C5C">/</span><span style="color:black">Km</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif""> </span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><i><span lang="EN-US" style="font-family:"Monospaced","serif";color:#64AE64">// using vectorization</span></i><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">a</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:#4A55DB">[</span><span style="color:rosybrown">0</span><span style="color:#FFAA00">:</span><span style="color:black">n</span><span style="color:#4A55DB">]</span><span style="color:#5C5C5C">'</span><span style="color:black">;</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">b</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:#32B9B9">ones</span><span style="color:#4A55DB">(</span><span style="color:black">n,</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:black">;</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">i</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:black">a</span><span style="color:#5C5C5C">.*.</span><span style="color:black">b;</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">j</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:black">b</span><span style="color:#5C5C5C">.*.</span><span style="color:black">a;</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif""> </span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">Km2</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:#32B9B9">zeros</span><span style="color:#4A55DB">(</span><span style="color:black">n</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">(</span><span style="color:black">n</span><span style="color:#5C5C5C">+</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:black">,</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:black">;</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><i><span lang="EN-US" style="font-family:"Monospaced","serif";color:#64AE64">//toto = i. *j</span></i><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">Km2</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:#4A55DB">(</span> <span style="color:orchid">%pi</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">6</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">((</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:black">i</span><span style="color:#5C5C5C">+</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">2.</span> <span style="color:#5C5C5C">*</span> <span style="color:black">j</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">2</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">((</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*.</span> <span style="color:black">i</span><span style="color:#5C5C5C">+</span><span style="color:rosybrown">1</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">(</span><span style="color:black">b</span><span style="color:#5C5C5C">/</span><span style="color:rosybrown">2</span><span style="color:#5C5C5C">*</span><span style="color:black">d</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">2</span> <span style="color:#5C5C5C">+</span> <span style="color:black">j</span><span style="color:#5C5C5C">^</span><span style="color:rosybrown">2</span><span style="color:#4A55DB">)</span> <span style="color:#4A55DB">)</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">Km2</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:#4A55DB">(</span><span style="color:rosybrown">144</span><span style="color:#5C5C5C">*</span><span style="color:#4A55DB">(</span><span style="color:black">b</span><span style="color:#5C5C5C">/</span><span style="color:black">d</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">.^</span><span style="color:rosybrown">4</span><span style="color:#4A55DB">)</span><span style="color:#5C5C5C">/</span> <span style="color:black">Km2</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span lang="EN-US" style="font-family:"Monospaced","serif";color:black">Km2</span><span lang="EN-US" style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">1.2</span> <span style="color:#5C5C5C">+</span> <span style="color:#4A55DB">(</span><span style="color:black">nu</span><span style="color:#5C5C5C">/</span><span style="color:#4A55DB">(</span><span style="color:rosybrown">1</span><span style="color:#5C5C5C">+</span><span style="color:black">nu</span><span style="color:#4A55DB">))</span><span style="color:#5C5C5C">*</span><span style="color:black">Km2</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<pre><span style="font-family:"Monospaced","serif";color:black">k2</span><span style="font-family:"Monospaced","serif""> <span style="color:#5C5C5C">=</span> <span style="color:rosybrown">1</span><span style="color:#5C5C5C">/</span><span style="color:black">Km2</span></span><span style="font-size:10.0pt;font-family:"Courier New""><o:p></o:p></span></pre>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal" style="mso-margin-top-alt:auto;mso-margin-bottom-alt:auto"><span lang="EN-US" style="font-size:7.5pt;font-family:"Times New Roman","serif";color:red"> </span><o:p></o:p></p>
<p><b><i><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:gray">EXPORT CONTROL :
<br>
</span></i></b><b><span style="font-size:10.0pt;font-family:"Arial","sans-serif";color:black">Cet email ne contient pas de données techniques<br>
This email does not contain technical data</span></b><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</body>
</html>