<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:Wingdings;
panose-1:5 0 0 0 0 0 0 0 0 0;}
@font-face
{font-family:Wingdings;
panose-1:5 0 0 0 0 0 0 0 0 0;}
@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: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";
color:black;
mso-fareast-language:EN-US;}
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";
color:black;}
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";
color:black;}
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";
color:black;
mso-fareast-language:EN-US;}
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;}
span.EmailStyle21
{mso-style-type:personal;
font-family:"Calibri","sans-serif";
color:windowtext;}
span.EmailStyle22
{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;}
/* List Definitions */
@list l0
{mso-list-id:666054383;
mso-list-type:hybrid;
mso-list-template-ids:2131676516 1599924016 67895299 67895301 67895297 67895299 67895301 67895297 67895299 67895301;}
@list l0:level1
{mso-level-start-at:13;
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;
mso-bidi-font-family:"Times New Roman";}
@list l0:level2
{mso-level-number-format:bullet;
mso-level-text:o;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;
font-family:"Courier New";}
@list l0:level3
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;
font-family:Wingdings;}
@list l0:level4
{mso-level-number-format:bullet;
mso-level-text:\F0B7;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;
font-family:Symbol;}
@list l0:level5
{mso-level-number-format:bullet;
mso-level-text:o;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;
font-family:"Courier New";}
@list l0:level6
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;
font-family:Wingdings;}
@list l0:level7
{mso-level-number-format:bullet;
mso-level-text:\F0B7;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;
font-family:Symbol;}
@list l0:level8
{mso-level-number-format:bullet;
mso-level-text:o;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;
font-family:"Courier New";}
@list l0:level9
{mso-level-number-format:bullet;
mso-level-text:\F0A7;
mso-level-tab-stop:none;
mso-level-number-position:left;
text-indent:-18.0pt;
font-family:Wingdings;}
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]-->
</head>
<body bgcolor="white" lang="FR" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span style="color:#1F497D">Hi Stéphane<o:p></o:p></span></p>
<p class="MsoNormal"><span style="color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="color:#1F497D">I used the sqrt function in the idea of the “Sum of the Square Error” (while there’s not sum here) .. but I understand your remarks<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">Thanks<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">Paul<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"">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";color:windowtext;mso-fareast-language:FR">De :</span></b><span style="font-size:10.0pt;font-family:"Tahoma","sans-serif";color:windowtext;mso-fareast-language:FR">
users [mailto:users-bounces@lists.scilab.org] <b>De la part de</b> Stéphane Mottelet<br>
<b>Envoyé :</b> lundi 16 janvier 2017 10:45<br>
<b>À :</b> Users mailing list for Scilab<br>
<b>Objet :</b> [EXTERNAL] Re: [Scilab-users] Optim & NelderMead use [closed]<o:p></o:p></span></p>
</div>
</div>
<p class="MsoNormal"><o:p> </o:p></p>
<p>Hi Paul,<o:p></o:p></p>
<p>your cost function<o:p></o:p></p>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">f</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">sqrt</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">((</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">val_lin </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">-</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> val_rac</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">^</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">2</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<p>hasn't changed, since sqrt(x^2)=abs(x). What I meant before is replacing<o:p></o:p></p>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">f</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">abs</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">val_lin </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">-</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> val_rac</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<p>by<o:p></o:p></p>
<p><span style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">f</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">
</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">
</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">val_lin
</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">-</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> val_rac</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">^</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">2</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;
</span><o:p></o:p></p>
<p class="MsoNormal" style="margin-bottom:12.0pt">in order to make it differentiable. When using a non-differentiable cost function together with numderivative, it seems logical that tweaking the step size could artificially help convergence.<br>
<br>
S.<br>
<br>
<o:p></o:p></p>
<div>
<p class="MsoNormal">Le 16/01/2017 à 10:30, Carrico, Paul a écrit :<o:p></o:p></p>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<p class="MsoNormal"><span lang="EN-US">Hi all</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">After performing tests (and modifying the target function as it should have been done first), I can better understand how to use ‘optim’ and ‘Neldermead’ procedures.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">For my needs the mean flags are :</span><o:p></o:p></p>
<pre style="margin-left:36.0pt;text-indent:-18.0pt;mso-list:l0 level1 lfo2"><![if !supportLists]><span style="font-family:"Calibri","sans-serif""><span style="mso-list:Ignore">-<span style="font:7.0pt "Times New Roman""> </span></span></span><![endif]><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";mso-fareast-language:EN-US">Step h in numderivative </span><span lang="EN-US" style="font-size:11.0pt;font-family:Wingdings;mso-fareast-language:EN-US">à</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";mso-fareast-language:EN-US"> usefull reading as “EE 221 Numerical Computing" Scott Hudson</span><o:p></o:p></pre>
<p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo2"><![if !supportLists]><span style="mso-list:Ignore">-<span style="font:7.0pt "Times New Roman"">
</span></span><![endif]><span lang="EN-US">The threshold epsg in optim (%eps is the default value – such high accuracy is not necessary for my application – furthermore using a value such as 1e-5 leads to err=1 that is correct for checking)</span><o:p></o:p></p>
<pre style="margin-left:36.0pt;text-indent:-18.0pt;mso-list:l0 level1 lfo2"><![if !supportLists]><span style="font-family:"Calibri","sans-serif""><span style="mso-list:Ignore">-<span style="font:7.0pt "Times New Roman""> </span></span></span><![endif]><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";mso-fareast-language:EN-US">Ditto for Nelder-Mead and ‘-tolfunrelative’ & ‘-tolxrelative’</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";mso-fareast-language:EN-US"> </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";mso-fareast-language:EN-US">Now it works fine :-)</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";mso-fareast-language:EN-US"> </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";mso-fareast-language:EN-US">Thanks all for the support</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";mso-fareast-language:EN-US"> </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Calibri","sans-serif";mso-fareast-language:EN-US">Paul</span><o:p></o:p></pre>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US">#####################################################################</span><o:p></o:p></p>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">mode</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">0</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">clear</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">global</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">count_use</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">count_use </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">0</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// ****</span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#B01813">function</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">f</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">lineaire</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">x</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">a2</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">b2</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">f</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">a2</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">*</span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">x</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">+</span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">b2</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#B01813">endfunction</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// ****</span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#B01813">function</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">g</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">racine</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">x</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">a1</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">b1</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">g</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">sqrt</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">a1</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">*</span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">x</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">+</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">b1</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#B01813">endfunction</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// ****</span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#B01813">function</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">f</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">target</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">x</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">a1</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">b1</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">a2</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">b2</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">val_lin </span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span style="font-size:11.0pt;font-family:"Monospaced","serif""> <u>lineaire</u></span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><b><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">x</span></b><span style="font-size:11.0pt;font-family:"Monospaced","serif"">,</span><b><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">a2</span></b><span style="font-size:11.0pt;font-family:"Monospaced","serif"">,</span><b><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">b2</span></b><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif""> val_rac </span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span style="font-size:11.0pt;font-family:"Monospaced","serif""> <u>racine</u></span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><b><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">x</span></b><span style="font-size:11.0pt;font-family:"Monospaced","serif"">,</span><b><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">a1</span></b><span style="font-size:11.0pt;font-family:"Monospaced","serif"">,</span><b><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">b1</span></b><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">f</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">sqrt</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">((</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">val_lin </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">-</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> val_rac</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">^</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">2</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">global</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">count_use</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> count_use </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> count_use </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">+</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#B01813">endfunction</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// Cost function: </span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#B01813">function</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">[</span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">f</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">g</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">ind</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">]</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">cost</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">x</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">ind</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">a1</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">b1</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">a2</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">b2</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">f</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> <u>target</u></span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">x</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// g = numderivative(target, x.',order = 4); </span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">g</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AE5CB0">numderivative</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">target</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#834310">x</span></b><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.'</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">,</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1e-3</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, order </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">4</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// 1E-3 => see EE 221 "Numerical Computing" Scott Hudson </span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// Study of the influence of h on the number of target function calculation & the fopt accuracy:</span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// (epsg = %eps here)</span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// h = 1.e-1 => number = 220 & fopt = 2.242026e-05</span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// h = 1.e-2 => number = 195 & fopt = 2.267564e-07</span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// h = 1.e-3 => number = 170 & fopt = 2.189495e-09 ++++++</span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// h = 1.e-4 => number = 190 & fopt = 1.941203e-11</span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// h = 1.e-5 => number = 215 & fopt = 2.131628e-13</span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// h = 1.e-6 => number = 235 & fopt = 0.</span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// h = 1.e-7 => number = 255 & fopt = 7.105427e-15</span></i><o:p></o:p></pre>
<pre><i><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// h = 1.e-8 => number = 275 & fopt = 0.</span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#B01813">endfunction</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// *************************</span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// optimisation with optim</span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">initial_parameters </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">[</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">10</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">]</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">lower_bounds </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">[</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">0</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">]</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">upper_bounds </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">[</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1000</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">]</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">nocf </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1000</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// number max of call of f</span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">niter </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1000</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// number max of iterations</span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a1 </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">30</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">b1 </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">2.5</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a2 </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">b2 </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">2</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">epsg </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1e-5</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// gradient norm threshold (%eps by defaut) --> lead to err = 1 !!!</span></i><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">//epsg = %eps; // lead to Err = 13</span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">epsf </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">0</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">//threshold controlling decreasing of f (epsf = 0 by defaut)</span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">costf </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">list</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">cost</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, a1, b1, a2, b2</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">[</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">fopt, xopt, gopt, work, iters, evals, err</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">]</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">optim</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">costf,</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">'b'</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">,lower_bounds,upper_bounds,initial_parameters,</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">'qn'</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">,</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">'ar'</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">,nocf,niter,epsg,epsf</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">printf</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">"Optimized value : %g\n"</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">,xopt</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">printf</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">"min cost function value (should be as closed as possible to 0) ; %e\n"</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">,fopt</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">printf</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">'Number of calculations = %d !!!\n'</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">,count_use</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><i><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// Curves definition</span></i><o:p></o:p></pre>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif"">x </span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><u><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AE5CB0">linspace</span></u><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">0</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">,</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">50</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">,</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1000</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">'</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif"">plot_raci </span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span style="font-size:11.0pt;font-family:"Monospaced","serif""> <u>racine</u></span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">x,a1,b1</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif"">plot_lin </span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span style="font-size:11.0pt;font-family:"Monospaced","serif""> <u>lineaire</u></span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">x,a2,b2</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><u><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AE5CB0">scf</span></u><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">drawlater</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">()</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">xgrid</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">3</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">f </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AE5CB0">gcf</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">()</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">//f </span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">f</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">figure_size</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">[</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1000</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">, </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1000</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">]</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">f</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">background</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">color</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">255</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">,</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">255</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">,</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">255</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AE5CB0">gca</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">()</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">//a </span></i><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">font_size</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">2</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">x_label</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">text</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">"X axis"</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> ; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">x_location</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">"bottom"</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">x_label</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">font_angle</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">0</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">x_label</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">font_size</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">4</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">y_label</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">text</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">"Y axis"</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">y_location</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">"left"</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">y_label</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">font_angle</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=-</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">90</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">Y_label</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">font_size</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">4</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">title</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">text</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">"Title"</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">title</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">font_size</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">5</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">a</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">line_style</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><i><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#64AE64">// Curves plot</span></i><o:p></o:p></pre>
<pre><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AE5CB0">plot</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">x,plot_lin</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">e1 </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AE5CB0">gce</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">()</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">p1 </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> e1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">children</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">p1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">thickness</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">p1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">line_style</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">p1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">foreground</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">3</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><o:p></o:p></pre>
<pre><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AE5CB0">plot</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">(</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">x,plot_raci</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">)</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">e2 </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AE5CB0">gce</span></u><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">()</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">p2 </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> e2</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">children</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">; </span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">p2</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">thickness</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">p2</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">line_style</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">1</span><span lang="EN-US" style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif"">p2</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">.</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#AAAAAA">foreground</span><span style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#5C5C5C">=</span><span style="font-size:11.0pt;font-family:"Monospaced","serif""> </span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:rosybrown">2</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<pre><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#32B9B9">drawnow</span><span style="font-size:11.0pt;font-family:"Monospaced","serif";color:#4A55DB">()</span><span style="font-size:11.0pt;font-family:"Monospaced","serif"">;</span><o:p></o:p></pre>
<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"">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>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman","serif";mso-fareast-language:FR"><br>
<br>
<br>
<o:p></o:p></span></p>
<pre>_______________________________________________<o:p></o:p></pre>
<pre>users mailing list<o:p></o:p></pre>
<pre><a href="mailto:users@lists.scilab.org">users@lists.scilab.org</a><o:p></o:p></pre>
<pre><a href="https://urldefense.proofpoint.com/v2/url?u=http-3A__lists.scilab.org_mailman_listinfo_users&d=DwMG-g&c=0hKVUfnuoBozYN8UvxPA-w&r=4TCz--8bXfJhZZvIxJAemAJyz7Vfx78XvgYu3LN7eLo&m=ukUEk5nhdzIuOGpMpafIK_3yNqraOC0llI0bD1BOJ-w&s=PA1LLL81Y0BFOpHSSlQteCC8Mkvdc6GJPtsGf2XA1k8&e=">http://lists.scilab.org/mailman/listinfo/users</a><o:p></o:p></pre>
</blockquote>
<p class="MsoNormal"><span style="font-size:12.0pt;font-family:"Times New Roman","serif";mso-fareast-language:FR"><o:p> </o:p></span></p>
</div>
</body>
</html>