<html><body><div style="color:#000; background-color:#fff; font-family:bookman old style, new york, times, serif;font-size:12pt"><div>Dear Scilab Developers,</div><div><br></div><div>I think that Scilab (version 5.4.0 - alpha - 1) for Mac OS X Lion 10.7.x has bug because it shows error message for linpro function (see below).</div><div><br></div><div>I really want to run this code on my Mac OS X Lion. It works properly in Scilab 5.3 (for Windows).</div><div><br></div><div>In case you wish the code file, I could send to you.</div><div><br></div><div>Thank you in advance.<br></div><div><br></div><div>Best Regards,</div><div>Dr. Reinaldo.</div><div><br></div><div><br></div> <div style="font-size: 12pt; font-family: 'bookman old style', 'new york', times, serif; "> <div style="font-size: 12pt; font-family: 'times new roman', 'new york', times, serif; "> <div dir="ltr"> <font size="2" face="Arial"> ----- Mensagem encaminhada -----<br> <b><span
style="font-weight:bold;">De:</span></b> Scilab Support <scilab.support@scilab.org><br> <b><span style="font-weight: bold;">Para:</span></b> Reinaldo <rei.listas@yahoo.com> <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Quarta-feira, 13 de Junho de 2012 11:13<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [Ticket#2012053010000018] Optimization on Scilab 5.4.0 (for Mac OS X Lion)<br> </font> </div> <br><meta http-equiv="x-dns-prefetch-control" content="off"><div id="yiv829437783"><div>Dear Reinaldo,<br>
<br>
You should ask your question on the Scilab mailing list:<br>
<a rel="nofollow" target="_blank" href="http://www.scilab.org/communities/developer_zone/tools/mailing_list">http://www.scilab.org/communities/developer_zone/tools/mailing_list</a><br>
<br>
Regards,<br>
Scilab Team<br>
--<br>
===================================================<br>
<div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-style: normal; font-variant: normal; font-weight: normal; font-size: 12px; line-height: normal; font-family: Helvetica; ">
The Scilab Consortium</div>
<div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-style: normal; font-variant: normal; font-weight: normal; font-size: 12px; line-height: normal; font-family: Helvetica; ">
Digiteo</div>
<div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-style: normal; font-variant: normal; font-weight: normal; font-size: 12px; line-height: normal; font-family: Helvetica; ">
Domaine de Voluceau</div>
<div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-style: normal; font-variant: normal; font-weight: normal; font-size: 12px; line-height: normal; font-family: Helvetica; ">
Rocquencourt - B.P. 105</div>
<div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-style: normal; font-variant: normal; font-weight: normal; font-size: 12px; line-height: normal; font-family: Helvetica; ">
78153 Le Chesnay Cedex France</div>
===================================================<br>
<br>
<br>
<br>
01.06.2012 08:53 - Reinaldo a écrit:<br>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt;">
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-size: 12pt; font-family: 'bookman old style', 'new york', times, serif; ">
<div>
<span>Hi Scilab Support,</span></div>
<div>
</div>
<div>
<span>I typed at Scilab console:</span></div>
<div>
</div>
<div>
<div>
<span><span>-> atomsInstall("quapro")</span></span></div>
<div>
</div>
<div>
<span><span>And re-started again Scilab. The result is shown as follows:</span></span></div>
<div>
</div>
<div>
</div>
<div>
</div>
<div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> scilab-5.4.0-alpha-1</font></span></span></div>
<div>
</div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> Consórcio Scilab (DIGITEO)</font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> Copyright (c) 1989-2011 (INRIA)</font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> Direitos reservados (c) 1989-2007 (ENPC)</font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> ___________________________________________ </font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff">Loading the inicialization:</font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> carregando o ambiente inicial</font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff">-->atomsInstall("quapro")</font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> ans =</font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></span></div>
<div>
<span><span><font class="yiv829437783Apple-style-span" color="#0000ff">!quapro 1.1-2 allusers SCI/contrib/quapro/1.1-2 I !</font></span></span></div>
<div>
</div>
<div>
</div>
<div>
</div>
<div>
<span><span>So, seems that quapro is installed without any problem (or warning) !</span></span></div>
</div>
<div>
</div>
<div>
<span>The first time I compile that example10, appears the message like:</span></div>
<div>
</div>
<div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->Aeq = [</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 11 53 5 5 29 1 0 0 0 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 3 6 5 1 34 0 1 0 0 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 1 0 0 0 0 0 0 1 0 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 1 0 0 0 0 0 0 1 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 1 0 0 0 0 0 0 1 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 0 1 0 0 0 0 0 0 1 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 0 0 1 0 0 0 0 0 0 1</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> ]; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->beq = [40 20 1 1 1 1 1]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->c = [-13 -16 -16 -14 -39 0 0 0 0 0 0 0]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->[n,p]=size(Aeq);</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->cc = [zeros(p,1);1];</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->AAeq = [Aeq,beq-Aeq*ones(p,1)];</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->bbeq = beq;</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->z0 = ones(p+1,1);</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->zopt=karmarkar(AAeq,bbeq,cc,z0,0,0.99)</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> zopt =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4651227 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.3447924 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4790839 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4846474 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.3723201 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.9937153 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.9969277 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5348773 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.6552076 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5209161 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5153526 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.6276799 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 1.00D-162 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->x0=zopt(1:p) </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> x0 =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4651227 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.3447924 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4790839 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4846474 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.3723201 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.9937153 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.9969277 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5348773 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.6552076 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5209161 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5153526 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.6276799 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->xopt=karmarkar(Aeq,beq,c,x0,1.e-10,0.999) </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> xopt =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 1.0017139 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.2009097 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.9994057 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 1.0011172 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.2881802 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 1.034D-13 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 3.512D-13 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 5.532D-16 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.7990812 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 8.466D-15 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 8.225D-14 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.7119165 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->f = c'*xopt </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> f =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> - 57.481994 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->p = [-13 -16 -16 -14 -39]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->C = [</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 11 53 5 5 29</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 3 6 5 1 34</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 1 0 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 1 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 1 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 0 1 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 0 0 1</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> ]; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->b = [40 20 1 1 1 1 1]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->ci=[0 0 0 0 0]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->cs=[%inf %inf %inf %inf %inf]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->[x1,lagr,f]=linpro(p,C,b,ci,cs) </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b> !--error 4 </b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b>Invalid variable: linpro</b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->[clean(x1) clean(xopt(1:size(p,1)))]</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b> !--error 4 </b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b>Invalid variable: x1</b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b> </b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->a=0;for i=1:size(p,1); a=a+x1(i);end;a </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b> !--error 4 </b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b>Invalid variable: x1</b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b> </b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->b=0;for i=1:size(p,1); b=b+xopt(i);end;b </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> b =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 3.4913266 </font></span></div>
<div>
<span> </span></div>
</div>
<div>
</div>
<div>
<span>Then, when I re-start the Scilab, the result changes as following described:</span></div>
<div>
</div>
<div>
</div>
<div>
</div>
<div>
<div>
<span> <font class="yiv829437783Apple-style-span" color="#0000ff"> ___________________________________________ </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> scilab-5.4.0-alpha-1</font></span></div>
<div>
</div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> Consórcio Scilab (DIGITEO)</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> Copyright (c) 1989-2011 (INRIA)</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> Direitos reservados (c) 1989-2007 (ENPC)</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> ___________________________________________ </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">Execução de iniciação:</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> carregando o ambiente inicial</font></span></div>
<div>
</div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">Start Quapro toolbox</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">Load macros</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">Load gateways</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">Load help</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">Load demos</font></span></div>
<div>
</div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->Aeq = [</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 11 53 5 5 29 1 0 0 0 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 3 6 5 1 34 0 1 0 0 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 1 0 0 0 0 0 0 1 0 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 1 0 0 0 0 0 0 1 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 1 0 0 0 0 0 0 1 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 0 1 0 0 0 0 0 0 1 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 0 0 1 0 0 0 0 0 0 1</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> ]; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->beq = [40 20 1 1 1 1 1]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->c = [-13 -16 -16 -14 -39 0 0 0 0 0 0 0]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->[n,p]=size(Aeq);</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->cc = [zeros(p,1);1];</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->AAeq = [Aeq,beq-Aeq*ones(p,1)];</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->bbeq = beq;</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->z0 = ones(p+1,1);</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->zopt=karmarkar(AAeq,bbeq,cc,z0,0,0.99)</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> zopt =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4651227 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.3447924 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4790839 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4846474 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.3723201 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.9937153 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.9969277 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5348773 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.6552076 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5209161 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5153526 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.6276799 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 1.00D-162 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->x0=zopt(1:p) </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> x0 =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4651227 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.3447924 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4790839 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.4846474 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.3723201 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.9937153 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.9969277 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5348773 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.6552076 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5209161 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.5153526 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.6276799 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->xopt=karmarkar(Aeq,beq,c,x0,1.e-10,0.999) </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> xopt =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 1.0017139 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.2009097 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.9994057 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 1.0011172 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.2881802 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 1.034D-13 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 3.512D-13 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 5.532D-16 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.7990812 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 8.466D-15 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 8.225D-14 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0.7119165 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->f = c'*xopt </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> f =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> - 57.481994 </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->p = [-13 -16 -16 -14 -39]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->C = [</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 11 53 5 5 29</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 3 6 5 1 34</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 1 0 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 1 0 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 1 0 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 0 1 0</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> 0 0 0 0 1</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> ]; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->b = [40 20 1 1 1 1 1]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->ci=[0 0 0 0 0]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->cs=[%inf %inf %inf %inf %inf]'; </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->[x1,lagr,f]=linpro(p,C,b,ci,cs) </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b>Warning:</b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b>/Applications/scilab-5.4.0-alpha-1.app/Contents/MacOS/share/scilab/contrib/quapro/1.1-2/macros/quapro.bin</b></font></span></div>
<div>
</div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> f =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 1. </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> lagr =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0. </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> x1 =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 0. </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><span class="yiv829437783Apple-style-span"><font class="yiv829437783Apple-style-span" color="#ff0000"><b>Warning: stack problem..., solved</b></font></span></span></div>
<div>
</div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">--></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->[clean(x1) clean(xopt(1:size(p,1)))]</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font><font class="yiv829437783Apple-style-span" color="#ff0000"><b> !--error 5 </b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b>Incoherence column / row dimension.</b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->a=0;for i=1:size(p,1); a=a+x1(i);end;a </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font><font class="yiv829437783Apple-style-span" color="#ff0000"><b> !--error 21 </b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#ff0000"><b>Invalid index.</b></font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff">-->b=0;for i=1:size(p,1); b=b+xopt(i);end;b </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> b =</font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> </font></span></div>
<div>
<span><font class="yiv829437783Apple-style-span" color="#0000ff"> 3.4913266 </font></span></div>
<div>
<span> </span></div>
</div>
<div>
</div>
<div>
<span>So, it's this ! :-(((</span></div>
<div>
</div>
<div>
<span>Seems that "quapro" is not so compatible to Scilab 5.4.0-alpha1 for Mac OS X Lion 10.7.4 (XCode 4.2) :-(</span></div>
<div>
</div>
<div>
<span>However, it did not show any error message when it was installed !</span></div>
<div>
</div>
<div>
<span>Thank you in advance !</span></div>
<div>
</div>
<div>
<span>All best,</span></div>
<div>
<span>Reinaldo.</span></div>
<div>
</div>
<div>
<span>PS: The same file run properly using Scilab 5.3.3 for Windows.</span></div>
<div>
</div>
</div>
<div>
</div>
<div style="font-size: 12pt; font-family: times, serif; ">
<div style="font-size: 12pt; font-family: times, serif; ">
<div dir="ltr">
<hr size="1">
<font face="Arial" size="2"><b><span style="font-weight:bold;">De:</span></b> Scilab Support <scilab.support@scilab.org><br>
<b><span style="font-weight:bold;">Para:</span></b> rei.listas@yahoo.com<br>
<b><span style="font-weight:bold;">Enviadas:</span></b> Quarta-feira, 30 de Maio de 2012 4:01<br>
<b><span style="font-weight:bold;">Assunto:</span></b> Re: [Ticket#2012053010000018] Optimization on Scilab 5.4.0 (for Mac OS X Lion)<br>
</font></div>
</div>
</div>
</div>
</div>
<br>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt;">
<div style="color: rgb(0, 0, 0); background-color: rgb(255, 255, 255); font-size: 12pt; font-family: 'bookman old style', 'new york', times, serif; ">
<div style="font-size: 12pt; font-family: times, serif; ">
<div style="font-size: 12pt; font-family: times, serif; ">
<div id="yiv829437783">
<div>
Dear Reinaldo,<br>
<br>
Linpro is now part of an ATOMS module: <a rel="nofollow" target="_blank" href="http://atoms.scilab.org/toolboxes/quapro" title="http://atoms.scilab.org/toolboxes/quapro">http://atoms.scilab.org/toolboxes/quapro</a><br>
<br>
Scilab Team<br>
--<br>
===================================================<br>
<div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-style: normal; font-variant: normal; font-weight: normal; font-size: 12px; line-height: normal; font-family: Helvetica; ">
The Scilab Consortium</div>
<div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-style: normal; font-variant: normal; font-weight: normal; font-size: 12px; line-height: normal; font-family: Helvetica; ">
Digiteo</div>
<div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-style: normal; font-variant: normal; font-weight: normal; font-size: 12px; line-height: normal; font-family: Helvetica; ">
Domaine de Voluceau</div>
<div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-style: normal; font-variant: normal; font-weight: normal; font-size: 12px; line-height: normal; font-family: Helvetica; ">
Rocquencourt - B.P. 105</div>
<div style="margin-top: 0px; margin-right: 0px; margin-bottom: 0px; margin-left: 0px; font-style: normal; font-variant: normal; font-weight: normal; font-size: 12px; line-height: normal; font-family: Helvetica; ">
78153 Le Chesnay Cedex France</div>
===================================================<br>
<br>
<br>
<b><span style="font-size:10.0pt;">De :</span></b> Reinaldo [mailto:<a rel="nofollow" ymailto="mailto:rei.listas@yahoo.com" target="_blank" href="mailto:rei.listas@yahoo.com">rei.listas@yahoo.com</a>]<br>
<div style="border:none;border-left:solid blue 1.5pt;padding:0cm 0cm 0cm 4.0pt;">
<div class="yiv829437783gmail_quote">
<div lang="FR">
<div>
<div>
<div style="border:none;border-top:solid #b5c4df 1.0pt;padding:3.0pt 0cm 0cm 0cm;">
<span style="font-size:10.0pt;"><b>Envoyé :</b> vendredi 25 mai 2012 14:06<br>
<b>À :</b> SABIRON William<br>
<b>Objet :</b> Re: Optimization on Scilab 5.4.0 (for Mac OS X Lion)</span></div>
</div>
<br>
<div>
<div>
<span style="">Hi William,</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">I sent a mail to Michaël, but he didn't work on Scilab. That''s a pity ! :-(</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">He's good guy and helps me on my doubts on optimization function.</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">I would like to know if we could discuss about this issue freely.</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">If positive, please answer me why "linpro function" does not work on Scilab 5.4.</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">I wrote a program and it ran well on Scilab previous version (i.e., 5.3). However, using Scilab 5.4, </span></div>
<div>
<span style="">the same program has "errors", as follows:</span></div>
<div>
<span style=""> </span></div>
<div>
<div>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// Example 10 - page 80 - "Operations Research: applications and algorithms" (Winston)</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// PROJECT SELECTION PROBLEM: find the optimum value of following LP model.</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// Maximize 13*x1 + 16*x2 + 16*x3 + 14*x4 + 39*x5</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// such as</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// 11*x1 + 53*x2 + 5*x3 + 5*x4 + 29*x5 <= 40</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// 3*x1 + 6*x2 + 5*x3 + x4 + 34*x5 <= 20</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x1 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x2 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x3 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x4 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x5 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// xi >= 0 (i = 1,...,5)</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// Converting to min LP problem</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// Minimize -13*x1 - 16*x2 - 16*x3 - 14*x4 - 39*x5</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// such as</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// 11*x1 + 53*x2 + 5*x3 + 5*x4 + 29*x5 <= 40</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// 3*x1 + 6*x2 + 5*x3 + x4 + 34*x5 <= 20</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x1 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x2 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x3 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x4 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x5 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// xi >= 0 (i = 1,...,5)</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// First resolution of LP problem: (THIS RESOLUTION WORKED !)</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// You can try the karmarkar function built in Scilab :</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// Use slack variables:</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// Minimize -13*x1 - 16*x2 - 16*x3 - 14*x4 - 39*x5</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// such as</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// 11*x1 + 53*x2 + 5*x3 + 5*x4 + 29*x5 + e1 = 40</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// 3*x1 + 6*x2 + 5*x3 + x4 + 34*x5 + e2 = 20</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x1 + e3 = 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x2 + e4 = 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x3 + e5 = 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x4 + e6 = 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x5 + e7 = 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// xi, ej >= 0 (i = 1,...,5) (j = 1,..., 7)</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// The following script solves the problem. Here, the initial guess x0 is given. </span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">Aeq </span><span style="color:#5c5c5c;">=</span><span style=""> </span><span style="color:#4a55db;">[</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">11</span><span style=""> </span><span style="color:rosybrown;">53</span><span style=""> </span><span style="color:rosybrown;">5</span><span style=""> </span><span style="color:rosybrown;">5</span><span style=""> </span><span style="color:rosybrown;">29</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">3</span><span style=""> </span><span style="color:rosybrown;">6</span><span style=""> </span><span style="color:rosybrown;">5</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">34</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:#4a55db;">]</span><span style="">; </span><i><span style="color:#01a801;">// coefficients of constraint matrix</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">beq </span><span style="color:#5c5c5c;">=</span><span style=""> </span><span style="color:#4a55db;">[</span><span style="color:rosybrown;">40</span><span style=""> </span><span style="color:rosybrown;">20</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">1</span><span style="color:#4a55db;">]</span><span style="color:#5c5c5c;">'</span><span style="">; </span><i><span style="color:#01a801;">// vector of independent terms</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">c </span><span style="color:#5c5c5c;">=</span><span style=""> </span><span style="color:#4a55db;">[</span><span style="color:#5c5c5c;">-</span><span style="color:rosybrown;">13</span><span style=""> </span><span style="color:#5c5c5c;">-</span><span style="color:rosybrown;">16</span><span style=""> </span><span style="color:#5c5c5c;">-</span><span style="color:rosybrown;">16</span><span style=""> </span><span style="color:#5c5c5c;">-</span><span style="color:rosybrown;">14</span><span style=""> </span><span style="color:#5c5c5c;">-</span><span style="color:rosybrown;">39</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style="">
</span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style="color:#4a55db;">]</span><span style="color:#5c5c5c;">'</span><span style="">; </span><i><span style="color:#01a801;">// coefficients of objective function</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="color:#4a55db;">[</span><span style="">n,p</span><span style="color:#4a55db;">]</span><span style="color:#5c5c5c;">=</span><span style="color:#32b9b9;">size</span><span style="color:#4a55db;">(</span><span style="">Aeq</span><span style="color:#4a55db;">)</span><span style="">;</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// Now, assume that the initial guess x0 is unknown.</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// To find a feasible point, we solve the following problem : </span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">cc </span><span style="color:#5c5c5c;">=</span><span style=""> </span><span style="color:#4a55db;">[</span><span style="color:#32b9b9;">zeros</span><span style="color:#4a55db;">(</span><span style="">p,</span><span style="color:rosybrown;">1</span><span style="color:#4a55db;">)</span><span style="">;</span><span style="color:rosybrown;">1</span><span style="color:#4a55db;">]</span><span style="">;</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">AAeq </span><span style="color:#5c5c5c;">=</span><span style=""> </span><span style="color:#4a55db;">[</span><span style="">Aeq,beq</span><span style="color:#5c5c5c;">-</span><span style="">Aeq</span><span style="color:#5c5c5c;">*</span><span style="color:#32b9b9;">ones</span><span style="color:#4a55db;">(</span><span style="">p,</span><span style="color:rosybrown;">1</span><span style="color:#4a55db;">)]</span><span style="">;</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">bbeq </span><span style="color:#5c5c5c;">=</span><span style=""> beq;</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">z0 </span><span style="color:#5c5c5c;">=</span><span style=""> </span><span style="color:#32b9b9;">ones</span><span style="color:#4a55db;">(</span><span style="">p</span><span style="color:#5c5c5c;">+</span><span style="color:rosybrown;">1</span><span style="">,</span><span style="color:rosybrown;">1</span><span style="color:#4a55db;">)</span><span style="">;</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">zopt</span><span style="color:#5c5c5c;">=</span><u><span style="color:#ae5cb0;">karmarkar</span></u><span style="color:#4a55db;">(</span><span style="">AAeq,bbeq,cc,z0,</span><span style="color:rosybrown;">0</span><span style="">,</span><span style="color:rosybrown;">0.99</span><span style="color:#4a55db;">)</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">x0</span><span style="color:#5c5c5c;">=</span><span style="">zopt</span><span style="color:#4a55db;">(</span><span style="color:rosybrown;">1</span><span style="color:#ffaa00;">:</span><span style="">p</span><span style="color:#4a55db;">)</span><span style=""> </span><i><span style="color:#01a801;">// This gives a very small value of zopt(p+1)</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// The previous script allows to produces the initial guess: x0</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// We now plug the initial guess x0 into the original problem and get: </span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">xopt</span><span style="color:#5c5c5c;">=</span><u><span style="color:#ae5cb0;">karmarkar</span></u><span style="color:#4a55db;">(</span><span style="">Aeq,beq,c,x0,</span><span style="color:rosybrown;">1.e-10</span><span style="">,</span><span style="color:rosybrown;">0.999</span><span style="color:#4a55db;">)</span><span style=""> </span><i><span style="color:#01a801;">// optimization result</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// where the components xopt(1:4) are the original variables and xopt(5:) are slack variables.</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">f </span><span style="color:#5c5c5c;">=</span><span style=""> c</span><span style="color:#5c5c5c;">'*</span><span style="">xopt </span><i><span style="color:#01a801;">// result of objetive function</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// Second resolution of LP problem: (THIS RESOLUTION WORKED !)</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// You can also try the "quapro" module and the linpro function.</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// Minimize -13*x1 - 16*x2 - 16*x3 - 14*x4 - 39*x5</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// such as</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// 11*x1 + 53*x2 + 5*x3 + 5*x4 + 29*x5 <= 40</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// 3*x1 + 6*x2 + 5*x3 + x4 + 34*x5 <= 20</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x1 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x2 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x3 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x4 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// x5 <= 1</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// xi >= 0 (i = 1,...,5)</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">p </span><span style="color:#5c5c5c;">=</span><span style=""> </span><span style="color:#4a55db;">[</span><span style="color:#5c5c5c;">-</span><span style="color:rosybrown;">13</span><span style=""> </span><span style="color:#5c5c5c;">-</span><span style="color:rosybrown;">16</span><span style=""> </span><span style="color:#5c5c5c;">-</span><span style="color:rosybrown;">16</span><span style=""> </span><span style="color:#5c5c5c;">-</span><span style="color:rosybrown;">14</span><span style=""> </span><span style="color:#5c5c5c;">-</span><span style="color:rosybrown;">39</span><span style="color:#4a55db;">]</span><span style="color:#5c5c5c;">'</span><span style="">; </span><i><span style="color:#01a801;">// coefficients of objective function</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">C </span><span style="color:#5c5c5c;">=</span><span style=""> </span><span style="color:#4a55db;">[</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">11</span><span style=""> </span><span style="color:rosybrown;">53</span><span style=""> </span><span style="color:rosybrown;">5</span><span style=""> </span><span style="color:rosybrown;">5</span><span style=""> </span><span style="color:rosybrown;">29</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">3</span><span style=""> </span><span style="color:rosybrown;">6</span><span style=""> </span><span style="color:rosybrown;">5</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">34</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">0</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">1</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span><span style="color:#4a55db;">]</span><span style="">; </span><i><span style="color:#01a801;">// coefficients of constraint matrix</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">b </span><span style="color:#5c5c5c;">=</span><span style=""> </span><span style="color:#4a55db;">[</span><span style="color:rosybrown;">40</span><span style=""> </span><span style="color:rosybrown;">20</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">1</span><span style=""> </span><span style="color:rosybrown;">1</span><span style="color:#4a55db;">]</span><span style="color:#5c5c5c;">'</span><span style="">; </span><i><span style="color:#01a801;">// vector of independent terms</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">ci</span><span style="color:#5c5c5c;">=</span><span style="color:#4a55db;">[</span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style=""> </span><span style="color:rosybrown;">0</span><span style="color:#4a55db;">]</span><span style="color:#5c5c5c;">'</span><span style="">; </span><i><span style="color:#01a801;">// initial values (boundary condition)</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">cs</span><span style="color:#5c5c5c;">=</span><span style="color:#4a55db;">[</span><span style="color:orchid;">%inf</span><span style=""> </span><span style="color:orchid;">%inf</span><span style=""> </span><span style="color:orchid;">%inf</span><span style=""> </span><span style="color:orchid;">%inf</span><span style=""> </span><span style="color:orchid;">%inf</span><span style="color:#4a55db;">]</span><span style="color:#5c5c5c;">'</span><span style="">; </span><i><span style="color:#01a801;">// final values</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="color:#4a55db;">[</span><span style="">x1,lagr,f</span><span style="color:#4a55db;">]</span><span style="color:#5c5c5c;">=</span><span style="">linpro</span><span style="color:#4a55db;">(</span><span style="">p,C,b,ci,cs</span><span style="color:#4a55db;">)</span><span style=""> </span><i><span style="color:#01a801;">// x1 represents the optimum solution</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// PS: </span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// [x,lagr,f]=linpro(p,C,b,ci,cs [,x0]) when minimize p'*x under the constraints C*x <= b</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><i><span style="color:#01a801;">// [x,lagr,f]=linpro(p,-C,-b,ci,cs [,x0]) when minimize p'*x under the constraints C*x >= b</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style=""> </span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="color:#4a55db;">[</span><span style="color:#32b9b9;">clean</span><span style="color:#4a55db;">(</span><span style="">x1</span><span style="color:#4a55db;">)</span><span style=""> </span><span style="color:#32b9b9;">clean</span><span style="color:#4a55db;">(</span><span style="">xopt</span><span style="color:#4a55db;">(</span><span style="color:rosybrown;">1</span><span style="color:#ffaa00;">:</span><span style="color:#32b9b9;">size</span><span style="color:#4a55db;">(</span><span style="">p,</span><span style="color:rosybrown;">1</span><span style="color:#4a55db;">)))]</span></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">a</span><span style="color:#5c5c5c;">=</span><span style="color:rosybrown;">0</span><span style="">;</span><span style="color:#a020f0;">for</span><span style=""> i</span><span style="color:#5c5c5c;">=</span><span style="color:rosybrown;">1</span><span style="color:#ffaa00;">:</span><span style="color:#32b9b9;">size</span><span style="color:#4a55db;">(</span><span style="">p,</span><span style="color:rosybrown;">1</span><span style="color:#4a55db;">)</span><span style="">; a</span><span style="color:#5c5c5c;">=</span><span style="">a</span><span style="color:#5c5c5c;">+</span><span style="">x1</span><span style="color:#4a55db;">(</span><span style="">i</span><span style="color:#4a55db;">)</span><span style="">;</span><span style="color:#a020f0;">end</span><span style="">;a </span><i><span style="color:#01a801;">//sum of optimal values using
linpro</span></i></pre>
<pre style="margin-right:0cm;margin-bottom:12.0pt;margin-left:0cm;background:white;"><span style="">b</span><span style="color:#5c5c5c;">=</span><span style="color:rosybrown;">0</span><span style="">;</span><span style="color:#a020f0;">for</span><span style=""> i</span><span style="color:#5c5c5c;">=</span><span style="color:rosybrown;">1</span><span style="color:#ffaa00;">:</span><span style="color:#32b9b9;">size</span><span style="color:#4a55db;">(</span><span style="">p,</span><span style="color:rosybrown;">1</span><span style="color:#4a55db;">)</span><span style="">; b</span><span style="color:#5c5c5c;">=</span><span style="">b</span><span style="color:#5c5c5c;">+</span><span style="">xopt</span><span style="color:#4a55db;">(</span><span style="">i</span><span style="color:#4a55db;">)</span><span style="">;</span><span style="color:#a020f0;">end</span><span style="">;b </span><i><span style="color:#01a801;">//sum of optimal values using
Karmarkar</span></i></pre>
</div>
<div>
<span style=""> </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">However, seems that some functions are modified in Scilab 5.4, aren't they ?</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">Before, that program runs correctly ! Now, in Scilab 5.4, the result is:</span></div>
<div>
<span style=""> </span></div>
<div>
<div>
<span style="">-></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->Aeq = [</span></div>
<div>
<span style="">--> 11 53 5 5 29 1 0 0 0 0 0 0</span></div>
<div>
<span style="">--> 3 6 5 1 34 0 1 0 0 0 0 0</span></div>
<div>
<span style="">--> 1 0 0 0 0 0 0 1 0 0 0 0</span></div>
<div>
<span style="">--> 0 1 0 0 0 0 0 0 1 0 0 0</span></div>
<div>
<span style="">--> 0 0 1 0 0 0 0 0 0 1 0 0</span></div>
<div>
<span style="">--> 0 0 0 1 0 0 0 0 0 0 1 0</span></div>
<div>
<span style="">--> 0 0 0 0 1 0 0 0 0 0 0 1</span></div>
<div>
<span style="">--> ]; </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->beq = [40 20 1 1 1 1 1]'; </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->c = [-13 -16 -16 -14 -39 0 0 0 0 0 0 0]'; </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->[n,p]=size(Aeq);</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->cc = [zeros(p,1);1];</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->AAeq = [Aeq,beq-Aeq*ones(p,1)];</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->bbeq = beq;</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->z0 = ones(p+1,1);</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->zopt=karmarkar(AAeq,bbeq,cc,z0,0,0.99)</span></div>
<div>
<span style=""> zopt =</span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> 0.4651227 </span></div>
<div>
<span style=""> 0.3447924 </span></div>
<div>
<span style=""> 0.4790839 </span></div>
<div>
<span style=""> 0.4846474 </span></div>
<div>
<span style=""> 0.3723201 </span></div>
<div>
<span style=""> 0.9937153 </span></div>
<div>
<span style=""> 0.9969277 </span></div>
<div>
<span style=""> 0.5348773 </span></div>
<div>
<span style=""> 0.6552076 </span></div>
<div>
<span style=""> 0.5209161 </span></div>
<div>
<span style=""> 0.5153526 </span></div>
<div>
<span style=""> 0.6276799 </span></div>
<div>
<span style=""> 1.00D-162 </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->x0=zopt(1:p) </span></div>
<div>
<span style=""> x0 =</span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> 0.4651227 </span></div>
<div>
<span style=""> 0.3447924 </span></div>
<div>
<span style=""> 0.4790839 </span></div>
<div>
<span style=""> 0.4846474 </span></div>
<div>
<span style=""> 0.3723201 </span></div>
<div>
<span style=""> 0.9937153 </span></div>
<div>
<span style=""> 0.9969277 </span></div>
<div>
<span style=""> 0.5348773 </span></div>
<div>
<span style=""> 0.6552076 </span></div>
<div>
<span style=""> 0.5209161 </span></div>
<div>
<span style=""> 0.5153526 </span></div>
<div>
<span style=""> 0.6276799 </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->xopt=karmarkar(Aeq,beq,c,x0,1.e-10,0.999) </span></div>
<div>
<span style=""> xopt =</span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> 1.0017139 </span></div>
<div>
<span style=""> 0.2009097 </span></div>
<div>
<span style=""> 0.9994057 </span></div>
<div>
<span style=""> 1.0011172 </span></div>
<div>
<span style=""> 0.2881802 </span></div>
<div>
<span style=""> 1.034D-13 </span></div>
<div>
<span style=""> 3.512D-13 </span></div>
<div>
<span style=""> 5.532D-16 </span></div>
<div>
<span style=""> 0.7990812 </span></div>
<div>
<span style=""> 8.466D-15 </span></div>
<div>
<span style=""> 8.225D-14 </span></div>
<div>
<span style=""> 0.7119165 </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->f = c'*xopt </span></div>
<div>
<span style=""> f =</span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> - 57.481994 </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->p = [-13 -16 -16 -14 -39]'; </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->C = [</span></div>
<div>
<span style="">--> 11 53 5 5 29</span></div>
<div>
<span style="">--> 3 6 5 1 34</span></div>
<div>
<span style="">--> 1 0 0 0 0</span></div>
<div>
<span style="">--> 0 1 0 0 0</span></div>
<div>
<span style="">--> 0 0 1 0 0</span></div>
<div>
<span style="">--> 0 0 0 1 0</span></div>
<div>
<span style="">--> 0 0 0 0 1</span></div>
<div>
<span style="">--> ]; </span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->b = [40 20 1 1 1 1 1]'; </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->ci=[0 0 0 0 0]'; </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->cs=[%inf %inf %inf %inf %inf]'; </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->[x1,lagr,f]=linpro(p,C,b,ci,cs) </span></div>
<div>
<span style=""> !--error 4 </span></div>
<div>
<span style="">Variável indefinida: linpro</span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">--></span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->[clean(x1) clean(xopt(1:size(p,1)))]</span></div>
<div>
<span style=""> !--error 4 </span></div>
<div>
<span style="">Variável indefinida: x1</span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->a=0;for i=1:size(p,1); a=a+x1(i);end;a </span></div>
<div>
<span style=""> !--error 4 </span></div>
<div>
<span style="">Variável indefinida: x1</span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">-->b=0;for i=1:size(p,1); b=b+xopt(i);end;b </span></div>
<div>
<span style=""> b =</span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> 3.4913266 </span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">Do you know why this program has errors ?</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">Please update me about the optimization functions (if there are modified in Scilab 5.4).</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">Thank you so much, Willian !!!</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">All best.</span></div>
<div>
<span style="">Reinaldo Golmia Dante.</span></div>
<div>
<span style=""> </span></div>
<div>
<span style="">PS: this program is attached in this mail in case you wish run it !</span></div>
</div>
</div>
<div>
<span style=""> </span></div>
<div>
<span style=""> </span></div>
<div>
<span style=""> </span></div>
<div>
<div>
<div>
<hr align="center" size="1" width="100%">
<br>
<b><span style="font-size:10.0pt;">De:</span></b><span style="font-size:10.0pt;"> "<a rel="nofollow" ymailto="mailto:michael.baudin-noreply@scilab.org" target="_blank" href="mailto:michael.baudin-noreply@scilab.org">michael.baudin-noreply@scilab.org</a>" <<a rel="nofollow" ymailto="mailto:michael.baudin-noreply@scilab.org" target="_blank" href="mailto:michael.baudin-noreply@scilab.org">michael.baudin-noreply@scilab.org</a>><br>
<b>Para:</b> <a rel="nofollow" ymailto="mailto:rei.listas@yahoo.com" target="_blank" href="mailto:rei.listas@yahoo.com">rei.listas@yahoo.com</a><br>
<b>Enviadas:</b> Sexta-feira, 25 de Maio de 2012 8:33<br>
<b>Assunto:</b> Re: Optimization on Scilab 5.4.0 (for Mac OS X Lion)</span></div>
<br>
<span style="">From: Michaël BAUDIN <<a rel="nofollow" ymailto="mailto:michael.baudin@scilab.org" target="_blank" href="mailto:michael.baudin@scilab.org">michael.baudin@scilab.org</a>><br>
<br>
Bonjour,<br>
<br>
Je ne fais plus partie des effectifs de la Fondation Scilab Digiteo.<br>
Merci de contacter <a rel="nofollow" ymailto="mailto:scilab.support@scilab.org" target="_blank" href="mailto:scilab.support@scilab.org">scilab.support@scilab.org</a> pour toute question relative à Scilab,<br>
William SABIRON (<a rel="nofollow" ymailto="mailto:william.sabiron@fcs-digiteotrianglephysique.fr" target="_blank" href="mailto:william.sabiron@fcs-digiteotrianglephysique.fr">william.sabiron@fcs-digiteotrianglephysique.fr</a>) sinon.<br>
Merci de ne plus utiliser cette adresse mail.<br>
<br>
Cordialement,<br>
<br>
Michaël BAUDIN<br>
<br>
-------- Original Message --------</span><br>
<br>
<br>
<br>
<br>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
<br>
</div></div><meta http-equiv="x-dns-prefetch-control" content="on"><br><br> </div> </div> </div></body></html>