<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>