﻿<?xml version="1.0" encoding="utf-8" standalone="yes"?>
<?application progid="SMath Studio Desktop" version="1.0.8253.4763"?>
<worksheet xmlns="http://smath.info/schemas/worksheet/1.0">
  <settings ppi="96">
    <identity>
      <id>e2f724d7-2e0d-4c4c-9159-c091642fd4cc</id>
      <revision>32</revision>
    </identity>
    <metadata lang="rus">
      <author>nvmnushkin</author>
    </metadata>
    <calculation>
      <precision>4</precision>
      <exponentialThreshold>5</exponentialThreshold>
      <trailingZeros>false</trailingZeros>
      <significantDigitsMode>false</significantDigitsMode>
      <roundingMode>0</roundingMode>
      <approximateEqualAccuracy>0</approximateEqualAccuracy>
      <fractions>decimal</fractions>
    </calculation>
    <pageModel active="true" viewMode="0" printGrid="false" printAreas="true" simpleEqualsOnly="false" printBackgroundImages="true">
      <paper id="9" orientation="Portrait" width="827" height="1169" />
      <margins left="39" right="39" top="49" bottom="49" />
      <header alignment="Center" color="#a9a9a9">  </header>
      <footer alignment="Center" color="#a9a9a9">&amp;[PAGENUM] / &amp;[COUNT]</footer>
      <backgrounds />
    </pageModel>
    <dependencies>
      <assembly name="SMath Studio Desktop" version="1.0.8253.4763" guid="a37cba83-b69c-4c71-9992-55ff666763bd" />
      <assembly name="MathRegion" version="1.11.8253.4763" guid="02f1ab51-215b-466e-a74d-5d8b1cf85e8d" />
      <assembly name="SpecialFunctions" version="1.12.8253.4763" guid="2814e667-4e12-48b1-8d51-194e480eabc5" />
      <assembly name="TextRegion" version="1.11.8253.4763" guid="485d28c5-349a-48b6-93be-12a35a1c1e39" />
      <assembly name="AreaRegion" version="1.11.8253.4763" guid="4974b228-4974-44cf-8274-bf2936b4a766" />
      <assembly name="Statistical Tools" version="1.0.6824.32561" guid="76f6e9f5-ac86-4ffb-9e4d-284c8daecaf9" />
      <assembly name="X-Y Plot Region (JXCharts)" version="0.2.8209.38826" guid="c12231ec-4873-43c1-a7d0-a167ebd17066" />
      <assembly name="PlotRegion" version="1.11.8253.4763" guid="c451c2b5-798b-4f08-b9ec-b90963d1ddaa" />
    </dependencies>
  </settings>
  <regions type="content">
    <region top="9" color="#000000">
      <area collapsed="false" />
      <region left="18" top="45" width="134" height="23" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p style="font-weight: bold;">Comments about #</p>
          </content>
        </text>
      </region>
      <region left="18" top="90" width="727" height="38" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>SMatth doesn't have true local variables, so it's usual put some decoration for avoid that issue. Maple have the sentence local ( ) and Mathematica have Block[ ] and Module[ ] for provide that functionality. For example</p>
          </content>
        </text>
      </region>
      <region left="72" top="135" width="116" height="26" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">x</e>
            <e type="function" args="1">f</e>
            <e type="operand">k</e>
            <e type="operand">x</e>
            <e type="operator" args="2">*</e>
            <e type="operand">1</e>
            <e type="operator" args="2">+</e>
            <e type="operator" args="2">:</e>
          </input>
        </math>
      </region>
      <region left="72" top="171" width="282" height="84" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">1</e>
            <e type="function" args="1">f</e>
            <e type="operand">x</e>
            <e type="operand">n</e>
            <e type="function" args="3">MySum</e>
            <e type="operand">s</e>
            <e type="operand">0</e>
            <e type="operator" args="2">:</e>
            <e type="operand">k</e>
            <e type="operand">1</e>
            <e type="operand">n</e>
            <e type="function" args="2">range</e>
            <e type="operand">s</e>
            <e type="operand">s</e>
            <e type="operand">x</e>
            <e type="function" args="1">f</e>
            <e type="operator" args="2">+</e>
            <e type="operator" args="2">:</e>
            <e type="function" args="3">for</e>
            <e type="operand">s</e>
            <e type="operand">3</e>
            <e type="operand">1</e>
            <e type="function" args="5">line</e>
            <e type="operator" args="2">:</e>
          </input>
        </math>
      </region>
      <region left="396" top="171" width="336" height="84" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">1</e>
            <e type="function" args="1">f#</e>
            <e type="operand">x#</e>
            <e type="operand">n#</e>
            <e type="function" args="3">MySum#</e>
            <e type="operand">s#</e>
            <e type="operand">0</e>
            <e type="operator" args="2">:</e>
            <e type="operand">k#</e>
            <e type="operand">1</e>
            <e type="operand">n#</e>
            <e type="function" args="2">range</e>
            <e type="operand">s#</e>
            <e type="operand">s#</e>
            <e type="operand">x#</e>
            <e type="function" args="1">f#</e>
            <e type="operator" args="2">+</e>
            <e type="operator" args="2">:</e>
            <e type="function" args="3">for</e>
            <e type="operand">s#</e>
            <e type="operand">3</e>
            <e type="operand">1</e>
            <e type="function" args="5">line</e>
            <e type="operator" args="2">:</e>
          </input>
        </math>
      </region>
      <region left="72" top="279" width="175" height="28" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">x</e>
            <e type="function" args="1">f</e>
            <e type="operand">2</e>
            <e type="operand">3</e>
            <e type="function" args="3">MySum</e>
          </input>
          <result action="symbolic">
            <e type="operand">15</e>
          </result>
        </math>
      </region>
      <region left="396" top="279" width="247" height="28" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">x</e>
            <e type="function" args="1">f</e>
            <e type="operand">2</e>
            <e type="operand">3</e>
            <e type="function" args="3">MySum#</e>
          </input>
          <result action="symbolic">
            <e type="operand">3</e>
            <e type="operand">1</e>
            <e type="operand">2</e>
            <e type="operand">k</e>
            <e type="operator" args="2">*</e>
            <e type="operator" args="2">+</e>
            <e type="bracket">(</e>
            <e type="operator" args="2">*</e>
          </result>
        </math>
      </region>
      <region left="18" top="324" width="157" height="23" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p style="font-weight: bold;">Comments about feval</p>
          </content>
        </text>
      </region>
      <region left="18" top="360" width="717" height="38" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>SMath have the equivalent of matlab function handler @ indicating the "signature" of the function, that is the number of paraters that it have. That enable to define f, f(x), f(x,y) as different things. But have a lot of issues handling 3 cases</p>
          </content>
        </text>
      </region>
      <region left="72" top="414" width="110" height="23" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>Functions with if</p>
          </content>
        </text>
      </region>
      <region left="72" top="441" width="264" height="23" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>Functions with functions defined in plugins</p>
          </content>
        </text>
      </region>
      <region left="72" top="468" width="179" height="23" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>Functions that call functions</p>
          </content>
        </text>
      </region>
      <region left="18" top="504" width="713" height="23" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>You can use Unknowns(f) for check if procedures like solve, rotts, diff, int and others goes to work or not. For example </p>
          </content>
        </text>
      </region>
      <region left="72" top="549" width="139" height="84" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">x</e>
            <e type="function" args="1">f</e>
            <e type="operand">x</e>
            <e type="operand">0</e>
            <e type="operator" args="2">&gt;</e>
            <e type="operand">2</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="operand">3</e>
            <e type="operator" args="2">-</e>
            <e type="operand">3</e>
            <e type="operand">x</e>
            <e type="operator" args="2">*</e>
            <e type="function" args="3">if</e>
            <e type="operator" args="2">:</e>
          </input>
        </math>
      </region>
      <region left="243" top="549" width="173" height="26" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">x</e>
            <e type="function" args="1">g</e>
            <e type="operand">5</e>
            <e type="operand">x</e>
            <e type="operand">4</e>
            <e type="function" args="2">Bessel</e>
            <e type="operator" args="2">*</e>
            <e type="operator" args="2">:</e>
          </input>
        </math>
      </region>
      <region left="450" top="558" width="240" height="158" color="#000000" fontSize="10">
        <plot type="2d" render="lines" scale_x="1.96773979691896" scale_y="8.95430243255239" scale_z="4.21802500966716" rotate_x="0" rotate_y="0" rotate_z="0" transpose_x="-63" transpose_y="-12" transpose_z="0">
          <input>
            <e type="operand">x</e>
            <e type="function" args="1">f</e>
            <e type="operand">x</e>
            <e type="function" args="1">g</e>
            <e type="operand">2</e>
            <e type="operand">1</e>
            <e type="function" args="4">sys</e>
          </input>
        </plot>
      </region>
      <region left="72" top="648" width="149" height="28" color="#000000" fontSize="10">
        <math error="2">
          <input>
            <e type="operand">x</e>
            <e type="function" args="1">f</e>
            <e type="function" args="1">Unknowns</e>
          </input>
          <result action="symbolic">
            <e type="operand">#</e>
          </result>
        </math>
      </region>
      <region left="72" top="693" width="166" height="28" color="#000000" fontSize="10">
        <math error="2">
          <input>
            <e type="operand">x</e>
            <e type="function" args="1">g</e>
            <e type="function" args="1">Unknowns</e>
          </input>
          <result action="symbolic">
            <e type="operand">x</e>
            <e type="operand">1</e>
            <e type="operand">1</e>
            <e type="function" args="3">mat</e>
          </result>
        </math>
      </region>
      <region left="72" top="738" width="242" height="24" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">lastError</e>
          </input>
          <result action="numeric">
            <e type="operand" style="string">x - not defined.</e>
          </result>
        </math>
      </region>
      <region left="18" top="783" width="721" height="38" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>Using this feval you can handle them, or more or less something like a workaround. Also, feval is like feval's old matlab, before the introduction of functions handlers.</p>
          </content>
        </text>
      </region>
      <region left="81" top="837" width="586" height="33" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">#fs#</e>
            <e type="operand">#x#</e>
            <e type="function" args="2">feval</e>
            <e type="operand">#fs#</e>
            <e type="function" args="1">num2str</e>
            <e type="operand" style="string">(</e>
            <e type="operand">#x#</e>
            <e type="function" args="1">num2str</e>
            <e type="operand" style="string">)</e>
            <e type="function" args="4">concat</e>
            <e type="function" args="1">str2num</e>
            <e type="operand">1</e>
            <e type="operand">1</e>
            <e type="function" args="3">line</e>
            <e type="operator" args="2">:</e>
          </input>
        </math>
      </region>
      <region left="81" top="882" width="121" height="26" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">f</e>
            <e type="operand">2</e>
            <e type="function" args="2">feval</e>
          </input>
          <result action="numeric">
            <e type="operand">5</e>
          </result>
        </math>
      </region>
      <region left="315" top="882" width="175" height="26" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand" style="string">g</e>
            <e type="operand">2</e>
            <e type="function" args="2">feval</e>
          </input>
          <result action="numeric">
            <e type="operand">1.8206</e>
          </result>
        </math>
      </region>
      <region left="18" top="918" width="147" height="23" color="#000000" fontSize="10" isBreakable="false">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p style="font-weight: bold;">Comments about Bis</p>
          </content>
        </text>
      </region>
      <region left="18" top="954" width="709" height="38" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>In SMath if you write a procedure and in its body you change the value of one of its parameters, then you change that variable globally. Actually, that's a great SMath feature. For example</p>
          </content>
        </text>
      </region>
      <region left="63" top="1017" width="128" height="45" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">x</e>
            <e type="operand">k</e>
            <e type="function" args="2">f</e>
            <e type="operand">k</e>
            <e type="operand">2</e>
            <e type="operator" args="2">:</e>
            <e type="operand">k</e>
            <e type="operand">x</e>
            <e type="operator" args="2">*</e>
            <e type="operand">2</e>
            <e type="operand">1</e>
            <e type="function" args="4">line</e>
            <e type="operator" args="2">:</e>
          </input>
        </math>
      </region>
      <region left="252" top="1017" width="46" height="24" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">k</e>
            <e type="operand">3</e>
            <e type="operator" args="2">:</e>
          </input>
        </math>
      </region>
      <region left="315" top="1017" width="89" height="26" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">1</e>
            <e type="operand">k</e>
            <e type="function" args="2">f</e>
          </input>
          <result action="numeric">
            <e type="operand">2</e>
          </result>
        </math>
      </region>
      <region left="432" top="1017" width="44" height="24" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">k</e>
          </input>
          <result action="numeric">
            <e type="operand">2</e>
          </result>
        </math>
      </region>
      <region left="18" top="1080" width="733" height="53" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>Following Bisection takes as arguments the string of an univariate function name, the bounds of some real interval where the function have a sign change and two epsilons: one for  check if the function value is enough small, εy, and another for check if  the interval where you are working is enough small.</p>
          </content>
        </text>
      </region>
      <region left="18" top="1152" width="729" height="53" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>About εy, notice that in numerical procedures the assertion a = b have not too much sense, it is just a casuality. Better is ask if abs(a-b) &lt; ε for some ε value. In the examples I use εy = 0, this is, check if for casuallity f(x) = 0, and try to find and interval if not.</p>
          </content>
        </text>
      </region>
      <region left="18" top="1215" width="447" height="23" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>About εx notice that its value determine the number of iterations, given by</p>
          </content>
        </text>
      </region>
      <region left="486" top="1215" width="121" height="43" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">N</e>
            <e type="operand">b</e>
            <e type="operand">a</e>
            <e type="operator" args="2">-</e>
            <e type="operand">ε</e>
            <e type="operator" args="2">/</e>
            <e type="operand">2</e>
            <e type="function" args="2">log</e>
            <e type="operator" args="2">≡</e>
          </input>
        </math>
      </region>
      <region left="18" top="1269" width="424" height="23" color="#000000" fontSize="10">
        <text lang="eng" fontFamily="Arial" fontSize="10">
          <content>
            <p>Finally, the returned value c indicates that f have a zero in the interval</p>
          </content>
        </text>
      </region>
      <region left="486" top="1269" width="129" height="24" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">c</e>
            <e type="operand">c</e>
            <e type="operand">εx</e>
            <e type="operator" args="2">-</e>
            <e type="operand">c</e>
            <e type="operand">εx</e>
            <e type="operator" args="2">+</e>
            <e type="function" args="3">lele</e>
          </input>
        </math>
      </region>
      <region left="54" top="1314" width="645" height="631" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">f#</e>
            <e type="operand">ao#</e>
            <e type="operand">bo#</e>
            <e type="operand">εy#</e>
            <e type="operand">εx#</e>
            <e type="function" args="5">Bis</e>
            <e type="operand" style="string">Avoid parameter modification</e>
            <e type="operand">a#</e>
            <e type="operand">b#</e>
            <e type="operand">1</e>
            <e type="operand">2</e>
            <e type="function" args="4">mat</e>
            <e type="operand">ao#</e>
            <e type="operand">bo#</e>
            <e type="operand">1</e>
            <e type="operand">2</e>
            <e type="function" args="4">mat</e>
            <e type="operator" args="2">:</e>
            <e type="operand" style="string">Initialize</e>
            <e type="operand">ya#</e>
            <e type="operand">yb#</e>
            <e type="operand">1</e>
            <e type="operand">2</e>
            <e type="function" args="4">mat</e>
            <e type="operand">f#</e>
            <e type="operand">a#</e>
            <e type="function" args="2">feval</e>
            <e type="operand">f#</e>
            <e type="operand">b#</e>
            <e type="function" args="2">feval</e>
            <e type="operand">1</e>
            <e type="operand">2</e>
            <e type="function" args="4">mat</e>
            <e type="operator" args="2">:</e>
            <e type="operand">ya#</e>
            <e type="function" args="1">sign</e>
            <e type="operand">yb#</e>
            <e type="function" args="1">sign</e>
            <e type="operator" args="2">≡</e>
            <e type="operand" style="string">Check interval</e>
            <e type="operand" style="string">Bad interval</e>
            <e type="function" args="1">error</e>
            <e type="operand">2</e>
            <e type="operand">1</e>
            <e type="function" args="4">line</e>
            <e type="operand">ya#</e>
            <e type="function" args="1">abs</e>
            <e type="operand">εy#</e>
            <e type="operand">ya#</e>
            <e type="function" args="1">UnitsOf</e>
            <e type="operator" args="2">*</e>
            <e type="operator" args="2">≤</e>
            <e type="operand" style="string">Check if one bound the solution</e>
            <e type="operand">c#</e>
            <e type="operand">a#</e>
            <e type="operator" args="2">:</e>
            <e type="operand">2</e>
            <e type="operand">1</e>
            <e type="function" args="4">line</e>
            <e type="operand">yb#</e>
            <e type="function" args="1">abs</e>
            <e type="operand">εy#</e>
            <e type="operand">yb#</e>
            <e type="function" args="1">UnitsOf</e>
            <e type="operator" args="2">*</e>
            <e type="operator" args="2">≤</e>
            <e type="operand" style="string">Check the other bound</e>
            <e type="operand">c#</e>
            <e type="operand">b#</e>
            <e type="operator" args="2">:</e>
            <e type="operand">2</e>
            <e type="operand">1</e>
            <e type="function" args="4">line</e>
            <e type="operand" style="string">Number of iterations</e>
            <e type="operand">N#</e>
            <e type="operand">1</e>
            <e type="operand">b#</e>
            <e type="operand">a#</e>
            <e type="operator" args="2">-</e>
            <e type="operand">b#</e>
            <e type="operand">a#</e>
            <e type="operator" args="2">-</e>
            <e type="function" args="1">UnitsOf</e>
            <e type="operator" args="2">/</e>
            <e type="function" args="1">abs</e>
            <e type="function" args="1">ln</e>
            <e type="operand">εx#</e>
            <e type="function" args="1">ln</e>
            <e type="operator" args="2">-</e>
            <e type="operand">2</e>
            <e type="function" args="1">ln</e>
            <e type="operator" args="2">/</e>
            <e type="operand">0</e>
            <e type="function" args="2">round</e>
            <e type="operator" args="2">+</e>
            <e type="operator" args="2">:</e>
            <e type="operand">iter#</e>
            <e type="operand">1</e>
            <e type="operand">N#</e>
            <e type="function" args="2">range</e>
            <e type="operand" style="string">Find the middle point</e>
            <e type="operand">c#</e>
            <e type="operand">0.5</e>
            <e type="operand">a#</e>
            <e type="operand">b#</e>
            <e type="operator" args="2">+</e>
            <e type="bracket">(</e>
            <e type="operator" args="2">*</e>
            <e type="function" args="1">eval</e>
            <e type="operator" args="2">:</e>
            <e type="operand">yc#</e>
            <e type="operand">f#</e>
            <e type="operand">c#</e>
            <e type="function" args="2">feval</e>
            <e type="operator" args="2">:</e>
            <e type="operand">1</e>
            <e type="operand">2</e>
            <e type="function" args="4">mat</e>
            <e type="operand">yc#</e>
            <e type="function" args="1">abs</e>
            <e type="operand">εy#</e>
            <e type="operand">yc#</e>
            <e type="function" args="1">UnitsOf</e>
            <e type="operator" args="2">*</e>
            <e type="operator" args="2">≤</e>
            <e type="operand" style="string">The function value is enough small</e>
            <e type="operand">break</e>
            <e type="operand">2</e>
            <e type="operand">1</e>
            <e type="function" args="4">line</e>
            <e type="operand">yb#</e>
            <e type="function" args="1">sign</e>
            <e type="operand">yc#</e>
            <e type="function" args="1">sign</e>
            <e type="operator" args="2">*</e>
            <e type="operand">0</e>
            <e type="operator" args="2">&gt;</e>
            <e type="operand" style="string">The solution is in the right interval</e>
            <e type="operand">b#</e>
            <e type="operand">yb#</e>
            <e type="operand">1</e>
            <e type="operand">2</e>
            <e type="function" args="4">mat</e>
            <e type="operand">c#</e>
            <e type="operand">yc#</e>
            <e type="operand">1</e>
            <e type="operand">2</e>
            <e type="function" args="4">mat</e>
            <e type="operator" args="2">:</e>
            <e type="operand">2</e>
            <e type="operand">1</e>
            <e type="function" args="4">line</e>
            <e type="operand" style="string">The solution is in the left interval</e>
            <e type="operand">a#</e>
            <e type="operand">ya#</e>
            <e type="operand">1</e>
            <e type="operand">2</e>
            <e type="function" args="4">mat</e>
            <e type="operand">c#</e>
            <e type="operand">yc#</e>
            <e type="operand">1</e>
            <e type="operand">2</e>
            <e type="function" args="4">mat</e>
            <e type="operator" args="2">:</e>
            <e type="operand">2</e>
            <e type="operand">1</e>
            <e type="function" args="4">line</e>
            <e type="function" args="3">if</e>
            <e type="function" args="3">if</e>
            <e type="operand">3</e>
            <e type="operand">1</e>
            <e type="function" args="5">line</e>
            <e type="function" args="3">for</e>
            <e type="operand">3</e>
            <e type="operand">1</e>
            <e type="function" args="5">line</e>
            <e type="function" args="7">if</e>
            <e type="operand">c#</e>
            <e type="operand">6</e>
            <e type="operand">1</e>
            <e type="function" args="8">line</e>
            <e type="operator" args="2">:</e>
          </input>
        </math>
      </region>
      <region left="18" top="1971" width="555" height="50" color="#000000" fontSize="10">
        <math>
          <input>
            <e type="operand">f#</e>
            <e type="operand">a#</e>
            <e type="operand">δx#</e>
            <e type="function" args="3">ndiff.1</e>
            <e type="operand">f#</e>
            <e type="operand">a#</e>
            <e type="operand">δx#</e>
            <e type="operand">a#</e>
            <e type="function" args="1">UnitsOf</e>
            <e type="operator" args="2">*</e>
            <e type="operator" args="2">+</e>
            <e type="function" args="2">feval</e>
            <e type="operand">f#</e>
            <e type="operand">a#</e>
            <e type="function" args="2">feval</e>
            <e type="operator" args="2">-</e>
            <e type="operand">δx#</e>
            <e type="operand">a#</e>
            <e type="function" args="1">UnitsOf</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" args="3">line</e>
            <e type="operator" args="2">:</e>
          </input>
        </math>
      </region>
      <region top="2025" color="#000000">
        <area terminator="true" />
      </region>
    </region>
    <region left="18" top="2052" width="90" height="23" color="#000000" fontSize="10">
      <text lang="eng" fontFamily="Arial" fontSize="10">
        <content>
          <p>Usual values</p>
        </content>
      </text>
    </region>
    <region left="162" top="2052" width="52" height="24" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand" style="unit">εy</e>
          <e type="operand">0</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region left="252" top="2052" width="90" height="33" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand" style="unit">εx</e>
          <e type="operand">10</e>
          <e type="operand">15</e>
          <e type="operator" args="1">-</e>
          <e type="operator" args="2">^</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region left="369" top="2052" width="82" height="33" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand" style="unit">δx</e>
          <e type="operand">10</e>
          <e type="operand">7</e>
          <e type="operator" args="1">-</e>
          <e type="operator" args="2">^</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region left="18" top="2106" width="61" height="23" color="#000000" fontSize="10">
      <text lang="eng" fontFamily="Arial" fontSize="10">
        <content>
          <p>function</p>
        </content>
      </text>
    </region>
    <region left="126" top="2106" width="347" height="29" color="#000000" bgColor="#ffff80" fontSize="10">
      <math decimalPlaces="2" exponentialThreshold="15">
        <input>
          <e type="operand">x</e>
          <e type="function" args="1">F</e>
          <e type="operand">x</e>
          <e type="operand">1</e>
          <e type="operand">x</e>
          <e type="function" args="2">BesselJ</e>
          <e type="operator" args="2">*</e>
          <e type="operand">1.01</e>
          <e type="operand">0</e>
          <e type="operand">x</e>
          <e type="function" args="2">BesselJ</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" args="3">line</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region left="18" top="2151" width="71" height="23" color="#000000" fontSize="10">
      <text lang="eng" fontFamily="Arial" fontSize="10">
        <content>
          <p>derivative</p>
        </content>
      </text>
    </region>
    <region left="117" top="2151" width="193" height="33" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand">x</e>
          <e type="function" args="1">f</e>
          <e type="operand">F</e>
          <e type="operand">x</e>
          <e type="operand" style="unit">δx</e>
          <e type="function" args="3">ndiff.1</e>
          <e type="operand">1</e>
          <e type="operand">1</e>
          <e type="function" args="3">line</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region left="18" top="2187" width="57" height="23" color="#000000" fontSize="10">
      <text lang="eng" fontFamily="Arial" fontSize="10">
        <content>
          <p>interval</p>
        </content>
      </text>
    </region>
    <region left="117" top="2187" width="114" height="27" color="#000000" bgColor="#ffff80" fontSize="10">
      <math>
        <input>
          <e type="operand">a</e>
          <e type="operand">b</e>
          <e type="operand">1</e>
          <e type="operand">2</e>
          <e type="function" args="4">mat</e>
          <e type="operand">0</e>
          <e type="operand">60</e>
          <e type="operand">1</e>
          <e type="operand">2</e>
          <e type="function" args="4">mat</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region left="18" top="2223" width="89" height="23" color="#000000" fontSize="10">
      <text lang="eng" fontFamily="Arial" fontSize="10">
        <content>
          <p>sub intervals</p>
        </content>
      </text>
    </region>
    <region left="117" top="2223" width="62" height="24" color="#000000" bgColor="#ffff80" fontSize="10">
      <math>
        <input>
          <e type="operand">N</e>
          <e type="operand">100</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region left="243" top="2223" width="160" height="41" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand">X</e>
          <e type="operand">a</e>
          <e type="operand">b</e>
          <e type="operand">a</e>
          <e type="operator" args="2">-</e>
          <e type="operand">N</e>
          <e type="operator" args="2">/</e>
          <e type="operand">0</e>
          <e type="operand">N</e>
          <e type="function" args="2">range</e>
          <e type="operator" args="2">*</e>
          <e type="operator" args="2">+</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region left="432" top="2223" width="46" height="24" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand">n</e>
          <e type="operand">0</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region left="486" top="2223" width="46" height="24" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand">m</e>
          <e type="operand">0</e>
          <e type="operator" args="2">:</e>
        </input>
      </math>
    </region>
    <region left="18" top="2268" width="601" height="23" color="#000000" fontSize="10">
      <text lang="eng" fontFamily="Arial" fontSize="10">
        <content>
          <p>Scan each subinterval for a root, and store it in a vector, detecting where there are a change of sign</p>
        </content>
      </text>
    </region>
    <region left="18" top="2304" width="342" height="83" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand">k</e>
          <e type="operand">1</e>
          <e type="operand">N</e>
          <e type="function" args="2">range</e>
          <e type="operand">X</e>
          <e type="operand">k</e>
          <e type="function" args="2">el</e>
          <e type="function" args="1">F</e>
          <e type="function" args="1">sign</e>
          <e type="operand">X</e>
          <e type="operand">k</e>
          <e type="operand">1</e>
          <e type="operator" args="2">+</e>
          <e type="function" args="2">el</e>
          <e type="function" args="1">F</e>
          <e type="function" args="1">sign</e>
          <e type="operator" args="2">*</e>
          <e type="operand">0</e>
          <e type="operator" args="2">≤</e>
          <e type="operand">xo</e>
          <e type="operand">n</e>
          <e type="operand">n</e>
          <e type="operand">1</e>
          <e type="operator" args="2">+</e>
          <e type="operator" args="2">:</e>
          <e type="function" args="2">el</e>
          <e type="operand">F</e>
          <e type="operand">X</e>
          <e type="operand">k</e>
          <e type="function" args="2">el</e>
          <e type="operand">X</e>
          <e type="operand">k</e>
          <e type="operand">1</e>
          <e type="operator" args="2">+</e>
          <e type="function" args="2">el</e>
          <e type="operand" style="unit">εy</e>
          <e type="operand" style="unit">εx</e>
          <e type="function" args="5">Bis</e>
          <e type="operator" args="2">:</e>
          <e type="function" args="2">if</e>
          <e type="function" args="3">for</e>
        </input>
      </math>
    </region>
    <region left="378" top="2304" width="342" height="83" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand">k</e>
          <e type="operand">1</e>
          <e type="operand">N</e>
          <e type="function" args="2">range</e>
          <e type="operand">X</e>
          <e type="operand">k</e>
          <e type="function" args="2">el</e>
          <e type="function" args="1">f</e>
          <e type="function" args="1">sign</e>
          <e type="operand">X</e>
          <e type="operand">k</e>
          <e type="operand">1</e>
          <e type="operator" args="2">+</e>
          <e type="function" args="2">el</e>
          <e type="function" args="1">f</e>
          <e type="function" args="1">sign</e>
          <e type="operator" args="2">*</e>
          <e type="operand">0</e>
          <e type="operator" args="2">≤</e>
          <e type="operand">uo</e>
          <e type="operand">m</e>
          <e type="operand">m</e>
          <e type="operand">1</e>
          <e type="operator" args="2">+</e>
          <e type="operator" args="2">:</e>
          <e type="function" args="2">el</e>
          <e type="operand">f</e>
          <e type="operand">X</e>
          <e type="operand">k</e>
          <e type="function" args="2">el</e>
          <e type="operand">X</e>
          <e type="operand">k</e>
          <e type="operand">1</e>
          <e type="operator" args="2">+</e>
          <e type="function" args="2">el</e>
          <e type="operand" style="unit">εy</e>
          <e type="operand" style="unit">εx</e>
          <e type="function" args="5">Bis</e>
          <e type="operator" args="2">:</e>
          <e type="function" args="2">if</e>
          <e type="function" args="3">for</e>
        </input>
      </math>
    </region>
    <region left="36" top="2403" width="616" height="208" color="#000000" bgColor="#ffffff" fontSize="10">
      <xyplot width="606" height="200" points="500" name="XYPlot">
        <chartstyle backcolor="White" bordercolor="Black" />
        <propertiessource index="1" sourcetype="PropertyGrid" />
        <grid gridcolor="LightGray" gridpattern="Dash" gridthickness="1" isxgrid="true" isygrid="true" isy2grid="false" />
        <xaxes xmin="-2.218934" xmax="60.88815" xtick="10" visible="True" decimalplaces="3" numberformat="General" />
        <yaxes ymin="-11.80097" ymax="9.597672" ytick="10" visible="True" decimalplaces="3" numberformat="General" />
        <y2axes isy2axis="false" y2min="-1" y2max="1" y2tick="0.5" visible="True" decimalplaces="3" numberformat="General" />
        <title2d title="" titlefont="Arial, 10pt" titlefontcolor="Black" />
        <xylabel labelfont="Arial, 10pt" labelfontcolor="Black" tickfont="Arial, 8pt" tickfontcolor="Black" xlabel="" ylabel="" y2label="y2" />
        <legend isbordervisible="true" islegendvisible="false" legendbackcolor="White" legendbordercolor="Black" legendfont="Arial, 8pt" legendposition="NorthEast" textcolor="Black" />
        <traces>
          <trace seriesname="" isy2data="false" isvisible="true" plotmethod="Lines" lineantialias="true" linecolor="Blue" linethickness="1" linepattern="Solid" symbolantialias="true" symbolsize="8" symboltype="None" symbolborderthickness="1" symbolbordercolor="Black" symbolfillcolor="White" />
          <trace seriesname="" isy2data="false" isvisible="true" plotmethod="Lines" lineantialias="true" linecolor="Red" linethickness="1" linepattern="Solid" symbolantialias="true" symbolsize="8" symboltype="None" symbolborderthickness="1" symbolbordercolor="Black" symbolfillcolor="White" />
          <trace seriesname="" isy2data="false" isvisible="true" plotmethod="Labels" lineantialias="true" linecolor="Green" linethickness="1" linepattern="Solid" symbolantialias="true" symbolsize="8" symboltype="None" symbolborderthickness="1" symbolbordercolor="Black" symbolfillcolor="White" />
          <trace seriesname="" isy2data="false" isvisible="true" plotmethod="Labels" lineantialias="true" linecolor="Fuchsia" linethickness="1" linepattern="Solid" symbolantialias="true" symbolsize="8" symboltype="None" symbolborderthickness="1" symbolbordercolor="Black" symbolfillcolor="White" />
        </traces>
        <input>
          <e type="operand">x</e>
          <e type="function" args="1">F</e>
          <e type="operand">x</e>
          <e type="function" args="1">f</e>
          <e type="operand">xo</e>
          <e type="operand">xo</e>
          <e type="function" args="1">F</e>
          <e type="function" args="1">vectorize</e>
          <e type="operand" style="string">o</e>
          <e type="function" args="3">augment</e>
          <e type="operand">uo</e>
          <e type="operand">uo</e>
          <e type="function" args="1">f</e>
          <e type="function" args="1">vectorize</e>
          <e type="operand" style="string">x</e>
          <e type="function" args="3">augment</e>
          <e type="operand">4</e>
          <e type="operand">1</e>
          <e type="function" args="6">sys</e>
        </input>
      </xyplot>
    </region>
    <region left="18" top="2736" width="60" height="24" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand">Alvaro</e>
        </input>
      </math>
    </region>
    <region left="315" top="2736" width="249" height="26" color="#000000" fontSize="10">
      <math>
        <input>
          <e type="operand">4</e>
          <e type="function" args="1">appVersion</e>
        </input>
        <result action="numeric">
          <e type="operand" style="string">1.0.8253.4763</e>
        </result>
      </math>
    </region>
  </regions>
</worksheet>