<html><body><div style="color:#000; background-color:#fff; font-family:bookman old style, new york, times, serif;font-size:12pt"><div><span>Hi Adrien,</span></div><div><span><br></span></div><div><span>Thank you for your answer and suggestion. I install it. Then I run it again and the results is shown, as follows:</span></div><div><span><br></span></div><div><span><div>-->Aeq = [</div><div>--> 11 53 5 5 29 1 0 0 0 0 0 0</div><div>--> 3 6 5 1 34 0 1 0 0 0 0 0</div><div>--> 1 0 0 0 0 0 0 1 0 0 0 0</div><div>--> 0 1 0 0 0 0 0 0 1 0 0 0</div><div>--> 0 0 1 0 0 0 0 0 0 1 0 0</div><div>--> 0 0 0 1 0 0 0 0 0 0 1 0</div><div>--> 0 0 0 0 1 0 0 0 0 0 0 1</div><div>--> ];
</div><div> </div><div>--></div><div> </div><div>-->beq = [40 20 1 1 1 1 1]'; </div><div> </div><div>-->c = [-13 -16 -16 -14 -39 0 0 0 0 0 0 0]'; </div><div> </div><div>-->[n,p]=size(Aeq);</div><div> </div><div>--></div><div> </div><div>--></div><div> </div><div>--></div><div> </div><div>-->cc = [zeros(p,1);1];</div><div> </div><div>-->AAeq = [Aeq,beq-Aeq*ones(p,1)];</div><div> </div><div>-->bbeq = beq;</div><div> </div><div>-->z0 = ones(p+1,1);</div><div> </div><div>-->zopt=karmarkar(AAeq,bbeq,cc,z0,0,0.99)</div><div> zopt =</div><div> </div><div> 0.4651227 </div><div> 0.3447924 </div><div> 0.4790839 </div><div> 0.4846474 </div><div> 0.3723201 </div><div> 0.9937153 </div><div>
0.9969277 </div><div> 0.5348773 </div><div> 0.6552076 </div><div> 0.5209161 </div><div> 0.5153526 </div><div> 0.6276799 </div><div> 1.00D-162 </div><div> </div><div>-->x0=zopt(1:p) </div><div> x0 =</div><div> </div><div> 0.4651227 </div><div> 0.3447924 </div><div> 0.4790839 </div><div> 0.4846474 </div><div> 0.3723201 </div><div> 0.9937153 </div><div> 0.9969277 </div><div> 0.5348773 </div><div> 0.6552076 </div><div> 0.5209161 </div><div> 0.5153526 </div><div> 0.6276799
</div><div> </div><div>--></div><div> </div><div>--></div><div> </div><div>-->xopt=karmarkar(Aeq,beq,c,x0,1.e-10,0.999) </div><div> xopt =</div><div> </div><div> 1.0017139 </div><div> 0.2009097 </div><div> 0.9994057 </div><div> 1.0011172 </div><div> 0.2881802 </div><div> 1.034D-13 </div><div> 3.512D-13 </div><div> 5.532D-16 </div><div> 0.7990812 </div><div> 8.466D-15 </div><div> 8.225D-14 </div><div> 0.7119165 </div><div> </div><div>--></div><div> </div><div>-->f = c'*xopt </div><div> f =</div><div> </div><div> - 57.481994 </div><div> </div><div>--></div><div> </div><div>-->p = [-13 -16 -16 -14
-39]'; </div><div> </div><div>--></div><div> </div><div>-->C = [</div><div>--> 11 53 5 5 29</div><div>--> 3 6 5 1 34</div><div>--> 1 0 0 0 0</div><div>--> 0 1 0 0 0</div><div>--> 0 0 1 0 0</div><div>--> 0 0 0 1 0</div><div>--> 0 0 0 0 1</div><div>--> ]; </div><div> </div><div>--></div><div> </div><div>--> </div><div> </div><div>-->b = [40 20 1 1 1 1 1]'; </div><div> </div><div>--></div><div> </div><div>-->ci=[0 0 0 0 0]'; </div><div> </div><div>--></div><div> </div><div>-->cs=[%inf %inf %inf %inf
%inf]'; </div><div> </div><div>--></div><div> </div><div>-->[x1,lagr,f]=linpro(p,C,b,ci,cs) </div><div>Warning:</div><div>/Applications/scilab-5.4.0-alpha-1.app/Contents/MacOS/share/scilab/contrib/quapro/1.1-2/macros/quapro.bin</div><div><br></div><div> f =</div><div> </div><div> 1. </div><div> lagr =</div><div> </div><div> 0. </div><div> x1 =</div><div> </div><div> 0. </div><div> </div><div>Warning: stack problem..., solved</div><div><br></div><div>--></div><div> </div><div>--></div><div> </div><div>--></div><div> </div><div>--></div><div> </div><div>--></div><div> </div><div>--></div><div> </div><div>-->[clean(x1) clean(xopt(1:size(p,1)))]</div><div>
!--error 5 </div><div>Incoerente rows/columns dimension.</div><div> </div><div> </div><div>-->a=0;for i=1:size(p,1); a=a+x1(i);end;a </div><div> !--error 21 </div><div>Invalid index.</div><div> </div><div> </div><div>-->b=0;for i=1:size(p,1); b=b+xopt(i);end;b </div><div> b =</div><div> </div><div> 3.4913266 </div><div> </div><div>In case you want, I could send to you the code in your mail.</div><div><br></div><div>All best,</div><div>Reinaldo.</div></span></div><div><br></div><div>PS: Interesting that it worked properly on previous version of Scilab.</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"> <hr size="1"> <b><span style="font-weight:bold;">De:</span></b> Adrien Vogt-Schilb <vogt@centre-cired.fr><br> <b><span style="font-weight: bold;">Para:</span></b> users@lists.scilab.org <br> <b><span style="font-weight: bold;">Enviadas:</span></b> Sexta-feira, 25 de Maio de 2012 9:29<br> <b><span style="font-weight: bold;">Assunto:</span></b> Re: [scilab-Users] Problem on optimization function (Scilab 5.4)<br> </font> </div> <br><meta http-equiv="x-dns-prefetch-control" content="off"><div id="yiv542005231">
<div>
On 25/05/2012 14:19, Reinaldo wrote:
<blockquote type="cite">
<pre style="font-family: Monospaced; "><span>-->[x1,lagr,f]=linpro(p,C,b,ci,cs)
!--error 4
Variável indefinida: linpro</span></pre>
</blockquote>
Hi<br>
<br>
you have to install the quapro module on your new scilab<br>
you can do so from the scilab console<br>
<br>
atomsInstall('quapro') <br>
<br>
<br>
<a rel="nofollow" class="yiv542005231moz-txt-link-freetext" target="_blank" href="http://atoms.scilab.org/toolboxes/quapro">http://atoms.scilab.org/toolboxes/quapro</a><br>
<br>
Hope this helps<br>
</div>
</div><meta http-equiv="x-dns-prefetch-control" content="on"><br><br> </div> </div> </div></body></html>