﻿<?xml version="1.0" encoding="utf-8" standalone="yes"?>
<?application progid="SMath Studio Desktop" version="0.97.5346.24640"?>
<regions>
  <settings>
    <identity>
      <id>ffd869cb-8fa9-433f-9af4-3c63373d1c83</id>
      <revision>5</revision>
    </identity>
    <calculation>
      <precision>4</precision>
      <exponentialThreshold>5</exponentialThreshold>
      <fractions>decimal</fractions>
    </calculation>
    <pageModel active="true" printAreas="true" simpleEqualsOnly="false" printBackgroundImages="true">
      <paper id="9" orientation="Portrait" width="827" height="1169" />
      <margins left="39" right="39" top="39" bottom="39" />
      <header alignment="Center" color="#a9a9a9">&amp;[DATE] &amp;[TIME] - &amp;[FILENAME]</header>
      <footer alignment="Center" color="#a9a9a9">&amp;[PAGENUM] / &amp;[COUNT]</footer>
      <backgrounds />
    </pageModel>
    <dependences>
      <assembly name="SMath Studio Desktop" version="0.97.5346.24640" guid="a37cba83-b69c-4c71-9992-55ff666763bd" />
      <assembly name="Text Region" version="1.10.5346.31409" guid="485d28c5-349a-48b6-93be-12a35a1c1e39" />
      <assembly name="Special Functions" version="1.11.5346.31403" guid="2814e667-4e12-48b1-8d51-194e480eabc5" />
      <assembly name="Math Region" version="0.97.5346.24640" guid="02f1ab51-215b-466e-a74d-5d8b1cf85e8d" />
      <assembly name="Nonlinear Solvers" version="1.1.5976.11383" guid="618b7e96-330a-406a-b055-9e577672f0b8" />
    </dependences>
  </settings>
  <region id="0" left="0" top="9" width="103" height="24" color="#000000" bgColor="#ffff80" fontSize="10">
    <text lang="eng">
      <p>References:</p>
    </text>
  </region>
  <region id="1" left="9" top="36" width="614" height="24" color="#000000" bgColor="#ffff80" fontSize="10">
    <text lang="eng">
      <p>http://www2.imm.dtu.dk/pubdb/views/edoc_download.php/3215/pdf/imm3215.pdf</p>
    </text>
  </region>
  <region id="2" left="9" top="63" width="556" height="24" color="#000000" bgColor="#ffff80" fontSize="10">
    <text lang="eng">
      <p>http://www2.imm.dtu.dk/pubdb/views/publication_details.php?id=3215</p>
    </text>
  </region>
  <region id="3" left="9" top="90" width="292" height="24" color="#000000" bgColor="#ffff80" fontSize="10">
    <text lang="eng">
      <p>http://www.brnt.eu/phd/node10.html</p>
    </text>
  </region>
  <region id="4" left="0" top="126" width="961" height="442" color="#000000" bgColor="#ffffff" fontSize="10">
    <math>
      <description active="true" position="Top" lang="eng">
        <p>Newton Method optimization algorithm.</p>
      </description>
      <input>
        <e type="operand">1</e>
        <e type="function" args="1">f</e>
        <e type="operand">x.0</e>
        <e type="operand">ε</e>
        <e type="operand">maxIters</e>
        <e type="operand">StepLength</e>
        <e type="function" args="5">NewtonMethod.ref</e>
        <e type="operand">m</e>
        <e type="operand">x.0</e>
        <e type="function" args="1">f</e>
        <e type="function" preserve="true" args="1">rows</e>
        <e type="function" preserve="true" args="1">eval</e>
        <e type="operator" args="2">:</e>
        <e type="operand">x#</e>
        <e type="function" args="1">F</e>
        <e type="operand">1</e>
        <e type="operand">2</e>
        <e type="operator" args="2">/</e>
        <e type="operand">x#</e>
        <e type="function" args="1">f</e>
        <e type="operand">i</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="bracket">(</e>
        <e type="operand">2</e>
        <e type="operator" args="2">^</e>
        <e type="bracket">(</e>
        <e type="operand">i</e>
        <e type="operand">1</e>
        <e type="operand">m</e>
        <e type="function" preserve="true" args="4">sum</e>
        <e type="operator" args="2">*</e>
        <e type="operator" args="2">:</e>
        <e type="operand">k</e>
        <e type="operand">0</e>
        <e type="operator" args="2">:</e>
        <e type="operand">x</e>
        <e type="operand">x.0</e>
        <e type="operator" args="2">:</e>
        <e type="operand">k</e>
        <e type="operand">maxIters</e>
        <e type="operator" args="2">&lt;</e>
        <e type="operand">h.N</e>
        <e type="operand">x#</e>
        <e type="function" args="1">F</e>
        <e type="operand">x</e>
        <e type="function" preserve="true" args="2">Hessian</e>
        <e type="bracket">(</e>
        <e type="operand">1</e>
        <e type="operator" args="1">-</e>
        <e type="operator" args="2">^</e>
        <e type="operand">x#</e>
        <e type="function" args="1">F</e>
        <e type="operand">x</e>
        <e type="function" preserve="true" args="2">Gradient</e>
        <e type="operator" args="1">-</e>
        <e type="bracket">(</e>
        <e type="operator" args="2">*</e>
        <e type="function" preserve="true" args="1">eval</e>
        <e type="operator" args="2">:</e>
        <e type="operand">h.N</e>
        <e type="operand">2</e>
        <e type="operator" args="2">^</e>
        <e type="function" preserve="true" args="1">sqrt</e>
        <e type="operand">ε</e>
        <e type="operator" args="2">&lt;</e>
        <e type="operand">break</e>
        <e type="operand">StepLength</e>
        <e type="operand" style="string">GSS</e>
        <e type="operator" args="2">≡</e>
        <e type="operand">α</e>
        <e type="operand">x</e>
        <e type="operand">α#</e>
        <e type="operand">h.N</e>
        <e type="operator" args="2">*</e>
        <e type="operator" args="2">+</e>
        <e type="function" args="1">F</e>
        <e type="operand">0</e>
        <e type="operand">10</e>
        <e type="operand">10</e>
        <e type="operand">2</e>
        <e type="operator" args="1">-</e>
        <e type="operator" args="2">^</e>
        <e type="function" preserve="true" args="4">GoldenSectionSearch.min</e>
        <e type="function" preserve="true" args="1">eval</e>
        <e type="operator" args="2">:</e>
        <e type="operand">α</e>
        <e type="operand">1</e>
        <e type="operator" args="2">:</e>
        <e type="function" preserve="true" args="3">if</e>
        <e type="operand">x</e>
        <e type="operand">x</e>
        <e type="operand">α</e>
        <e type="operand">h.N</e>
        <e type="operator" args="2">*</e>
        <e type="operator" args="2">+</e>
        <e type="function" preserve="true" args="1">eval</e>
        <e type="operator" args="2">:</e>
        <e type="operand">k</e>
        <e type="operand">k</e>
        <e type="operand">1</e>
        <e type="operator" args="2">+</e>
        <e type="operator" args="2">:</e>
        <e type="operand">3</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="5">line</e>
        <e type="function" preserve="true" args="3">if</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">line</e>
        <e type="function" preserve="true" args="2">while</e>
        <e type="operand">x</e>
        <e type="operand">k</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="3">mat</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
        <e type="operand">6</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="8">line</e>
        <e type="operator" args="2">:</e>
      </input>
    </math>
  </region>
  <region id="5" left="0" top="612" width="335" height="56" color="#000000" bgColor="#ffff80" fontSize="10">
    <math>
      <input>
        <e type="operand">x</e>
        <e type="function" args="1">Rosenbrock</e>
        <e type="operand">1</e>
        <e type="operand">x</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="operator" args="2">-</e>
        <e type="bracket">(</e>
        <e type="operand">2</e>
        <e type="operator" args="2">^</e>
        <e type="operand">10</e>
        <e type="operand">x</e>
        <e type="operand">2</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="operand">x</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="2">el</e>
        <e type="bracket">(</e>
        <e type="operand">2</e>
        <e type="operator" args="2">^</e>
        <e type="operator" args="2">-</e>
        <e type="bracket">(</e>
        <e type="operand">2</e>
        <e type="operator" args="2">^</e>
        <e type="operator" args="2">*</e>
        <e type="operator" args="2">+</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="3">mat</e>
        <e type="operator" args="2">:</e>
      </input>
    </math>
  </region>
  <region id="6" left="360" top="621" width="167" height="46" color="#000000" bgColor="#80ff80" fontSize="10">
    <math>
      <input>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
        <e type="function" args="1">Rosenbrock</e>
      </input>
      <result action="numeric">
        <e type="operand">0</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="3">mat</e>
      </result>
    </math>
  </region>
  <region id="7" left="18" top="684" width="518" height="66" color="#000000" bgColor="#ffffff" fontSize="10">
    <math decimalPlaces="5">
      <input>
        <e type="operand">x</e>
        <e type="function" args="1">Rosenbrock</e>
        <e type="operand">0.5</e>
        <e type="operator" args="1">-</e>
        <e type="operand">1.5</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
        <e type="operand">0.01</e>
        <e type="operand">50</e>
        <e type="operand">0</e>
        <e type="function" args="5">NewtonMethod.ref</e>
      </input>
      <result action="numeric">
        <e type="operand">0.98893</e>
        <e type="operand">0.97665</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
        <e type="operand">25</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="3">mat</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
      </result>
    </math>
  </region>
  <region id="8" left="18" top="765" width="551" height="66" color="#000000" bgColor="#ffffff" fontSize="10">
    <math decimalPlaces="5">
      <input>
        <e type="operand">x</e>
        <e type="function" args="1">Rosenbrock</e>
        <e type="operand">0.5</e>
        <e type="operator" args="1">-</e>
        <e type="operand">1.5</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
        <e type="operand">0.01</e>
        <e type="operand">50</e>
        <e type="operand" style="string">GSS</e>
        <e type="function" args="5">NewtonMethod.ref</e>
      </input>
      <result action="numeric">
        <e type="operand">0.99746</e>
        <e type="operand">0.99622</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
        <e type="operand">5</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="3">mat</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
      </result>
    </math>
  </region>
  <region id="9" left="0" top="846" width="299" height="44" color="#000000" bgColor="#ffff80" fontSize="10">
    <math>
      <input>
        <e type="operand">x</e>
        <e type="operand">y</e>
        <e type="function" args="2">Rosenbrock</e>
        <e type="operand">1</e>
        <e type="operand">x</e>
        <e type="operator" args="2">-</e>
        <e type="bracket">(</e>
        <e type="operand">2</e>
        <e type="operator" args="2">^</e>
        <e type="operand">10</e>
        <e type="operand">y</e>
        <e type="operand">x</e>
        <e type="operand">2</e>
        <e type="operator" args="2">^</e>
        <e type="operator" args="2">-</e>
        <e type="bracket">(</e>
        <e type="operand">2</e>
        <e type="operator" args="2">^</e>
        <e type="operator" args="2">*</e>
        <e type="operator" args="2">+</e>
        <e type="operator" args="2">:</e>
      </input>
    </math>
  </region>
  <region id="10" left="360" top="864" width="171" height="26" color="#000000" bgColor="#80ff80" fontSize="10">
    <math>
      <input>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" args="2">Rosenbrock</e>
      </input>
      <result action="numeric">
        <e type="operand">0</e>
      </result>
    </math>
  </region>
  <region id="11" left="18" top="918" width="521" height="66" color="#000000" bgColor="#ffffff" fontSize="10">
    <math>
      <input>
        <e type="operand">x</e>
        <e type="operand">y</e>
        <e type="function" args="2">Rosenbrock</e>
        <e type="operand">0.5</e>
        <e type="operator" args="1">-</e>
        <e type="operand">1.5</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">sys</e>
        <e type="operand">0.01</e>
        <e type="operand">50</e>
        <e type="operand">1</e>
        <e type="operand">0</e>
        <e type="operand">0</e>
        <e type="function" preserve="true" args="7">NewtonMethod</e>
      </input>
      <result action="numeric">
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">sys</e>
        <e type="operand">5</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="3">mat</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
      </result>
    </math>
  </region>
  <region id="12" left="18" top="999" width="596" height="66" color="#000000" bgColor="#ffffff" fontSize="10">
    <math decimalPlaces="5">
      <input>
        <e type="operand">x</e>
        <e type="operand">y</e>
        <e type="function" args="2">Rosenbrock</e>
        <e type="operand">0.5</e>
        <e type="operator" args="1">-</e>
        <e type="operand">1.5</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">sys</e>
        <e type="operand">0.01</e>
        <e type="operand">50</e>
        <e type="operand">1</e>
        <e type="operand">0</e>
        <e type="operand">0</e>
        <e type="function" preserve="true" args="7">NewtonMethod.GSS</e>
      </input>
      <result action="numeric">
        <e type="operand">0.99776</e>
        <e type="operand">0.99661</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">sys</e>
        <e type="operand">5</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="3">mat</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
      </result>
    </math>
  </region>
  <region id="13" left="18" top="1107" width="624" height="66" color="#000000" bgColor="#ffffff" fontSize="10">
    <math>
      <input>
        <e type="operand">x</e>
        <e type="operand">y</e>
        <e type="function" args="2">Rosenbrock</e>
        <e type="operand">0.5</e>
        <e type="operator" args="1">-</e>
        <e type="operand">1.5</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">sys</e>
        <e type="operand">10</e>
        <e type="operand">5</e>
        <e type="operator" args="1">-</e>
        <e type="operator" args="2">^</e>
        <e type="operand">0.01</e>
        <e type="operand">50</e>
        <e type="operand">1</e>
        <e type="operand">0</e>
        <e type="operand">0</e>
        <e type="function" preserve="true" args="8">NewtonMethod.CD</e>
      </input>
      <result action="numeric">
        <e type="operand">0.998</e>
        <e type="operand">0.996</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">sys</e>
        <e type="operand">6</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="3">mat</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
      </result>
    </math>
  </region>
  <region id="14" left="18" top="1188" width="674" height="66" color="#000000" bgColor="#ffffff" fontSize="10">
    <math decimalPlaces="5">
      <input>
        <e type="operand">x</e>
        <e type="operand">y</e>
        <e type="function" args="2">Rosenbrock</e>
        <e type="operand">0.5</e>
        <e type="operator" args="1">-</e>
        <e type="operand">1.5</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">sys</e>
        <e type="operand">10</e>
        <e type="operand">5</e>
        <e type="operator" args="1">-</e>
        <e type="operator" args="2">^</e>
        <e type="operand">0.01</e>
        <e type="operand">50</e>
        <e type="operand">1</e>
        <e type="operand">0</e>
        <e type="operand">0</e>
        <e type="function" preserve="true" args="8">NewtonMethod.GSS;CD</e>
      </input>
      <result action="numeric">
        <e type="operand">0.99823</e>
        <e type="operand">0.99627</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">sys</e>
        <e type="operand">50</e>
        <e type="operand">1</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="3">mat</e>
        <e type="operand">2</e>
        <e type="operand">1</e>
        <e type="function" preserve="true" args="4">mat</e>
      </result>
    </math>
  </region>
</regions>