Esempio: soluzione di un sistema di equazioni ODE del primo ordine

Fare clic per copiare questa espressione

<region id="ID0E2Y3BB" top="249.6" left="28.8" height="448.762999962962" width="485.073333333333" xml:lang="it" xmlns="http://schemas.mathsoft.com/worksheet50">
  <solveblock xml:lang="it">
    <regions xml:lang="it">
      <region id="" top="38.4" left="19.2" xml:lang="it" solve-block-category="constraint">
        <math xml:lang="it" disable-calc="false">
          <apply xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
            <equal xml:lang="it" />
            <apply xml:lang="it">
              <derivative xml:lang="it" />
              <lambda xml:lang="it">
                <boundVars xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">u</id>
                </boundVars>
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y0</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
              </lambda>
              <degree xml:lang="it">
                <placeholder xml:lang="it" />
              </degree>
            </apply>
            <apply xml:lang="it">
              <plus xml:lang="it" />
              <apply xml:lang="it">
                <mult xml:lang="it" />
                <real xml:lang="it">-8</real>
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y0</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
              </apply>
              <apply xml:lang="it">
                <mult xml:lang="it" />
                <real xml:lang="it">8</real>
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y1</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
              </apply>
            </apply>
          </apply>
        </math>
      </region>
      <region id="" top="96" left="19.2" xml:lang="it" solve-block-category="constraint">
        <math xml:lang="it" disable-calc="false">
          <apply xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
            <equal xml:lang="it" />
            <apply xml:lang="it">
              <id xml:space="preserve" xml:lang="it">y0</id>
              <real xml:lang="it">0</real>
            </apply>
            <real xml:lang="it">-1</real>
          </apply>
        </math>
      </region>
      <region id="" top="153.6" left="19.2" xml:lang="it" solve-block-category="constraint">
        <math xml:lang="it" disable-calc="false">
          <apply xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
            <equal xml:lang="it" />
            <apply xml:lang="it">
              <derivative xml:lang="it" />
              <lambda xml:lang="it">
                <boundVars xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">u</id>
                </boundVars>
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y1</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
              </lambda>
              <degree xml:lang="it">
                <placeholder xml:lang="it" />
              </degree>
            </apply>
            <apply xml:lang="it">
              <minus xml:lang="it" />
              <apply xml:lang="it">
                <plus xml:lang="it" />
                <apply xml:lang="it">
                  <mult xml:lang="it" />
                  <real xml:lang="it">30</real>
                  <apply xml:lang="it">
                    <id xml:space="preserve" xml:lang="it">y0</id>
                    <id xml:space="preserve" xml:lang="it">u</id>
                  </apply>
                </apply>
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y1</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
              </apply>
              <apply xml:lang="it">
                <mult xml:lang="it" />
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y0</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y2</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
              </apply>
            </apply>
          </apply>
        </math>
      </region>
      <region id="" top="211.2" left="19.2" xml:lang="it" solve-block-category="constraint">
        <math xml:lang="it" disable-calc="false">
          <apply xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
            <equal xml:lang="it" />
            <apply xml:lang="it">
              <id xml:space="preserve" xml:lang="it">y1</id>
              <real xml:lang="it">0</real>
            </apply>
            <real xml:lang="it">0</real>
          </apply>
        </math>
      </region>
      <region id="" top="268.8" left="19.2" xml:lang="it" solve-block-category="constraint">
        <math xml:lang="it" disable-calc="false">
          <apply xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
            <equal xml:lang="it" />
            <apply xml:lang="it">
              <derivative xml:lang="it" />
              <lambda xml:lang="it">
                <boundVars xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">u</id>
                </boundVars>
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y2</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
              </lambda>
              <degree xml:lang="it">
                <placeholder xml:lang="it" />
              </degree>
            </apply>
            <apply xml:lang="it">
              <minus xml:lang="it" />
              <apply xml:lang="it">
                <mult xml:lang="it" />
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y0</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y1</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
              </apply>
              <apply xml:lang="it">
                <mult xml:lang="it" />
                <apply xml:lang="it">
                  <div xml:lang="it" />
                  <real xml:lang="it">8</real>
                  <real xml:lang="it">3</real>
                </apply>
                <apply xml:lang="it">
                  <id xml:space="preserve" xml:lang="it">y2</id>
                  <id xml:space="preserve" xml:lang="it">u</id>
                </apply>
              </apply>
            </apply>
          </apply>
        </math>
      </region>
      <region id="" top="316.800012" left="19.2000000000001" xml:lang="it" solve-block-category="constraint">
        <math xml:lang="it" disable-calc="false">
          <apply xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
            <equal xml:lang="it" />
            <apply xml:lang="it">
              <id xml:space="preserve" xml:lang="it">y2</id>
              <real xml:lang="it">0</real>
            </apply>
            <real xml:lang="it">1</real>
          </apply>
        </math>
      </region>
      <region id="" top="403.200037" left="19.2" xml:lang="it" solve-block-category="solver">
        <math xml:lang="it" disable-calc="false">
          <define xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
            <matrix rows="3" cols="1" xml:lang="it">
              <id xml:space="preserve" xml:lang="it" label="Variable">f</id>
              <id xml:space="preserve" xml:lang="it" label="Variable">g</id>
              <id xml:space="preserve" xml:lang="it" label="Variable">h</id>
            </matrix>
            <apply xml:lang="it">
              <id xml:space="preserve" xml:lang="it" label="Keyword">Odesolve </id>
              <sequence xml:lang="it">
                <matrix rows="3" cols="1" xml:lang="it">
                  <apply xml:lang="it">
                    <id xml:space="preserve" xml:lang="it">y0</id>
                    <id xml:space="preserve" xml:lang="it">u</id>
                  </apply>
                  <apply xml:lang="it">
                    <id xml:space="preserve" xml:lang="it">y1</id>
                    <id xml:space="preserve" xml:lang="it">u</id>
                  </apply>
                  <apply xml:lang="it">
                    <id xml:space="preserve" xml:lang="it">y2</id>
                    <id xml:space="preserve" xml:lang="it">u</id>
                  </apply>
                </matrix>
                <id xml:space="preserve" xml:lang="it" label-is-contextual="true" label="Variable">T1</id>
                <apply xml:lang="it">
                  <mult xml:lang="it" />
                  <real xml:lang="it">1</real>
                  <apply xml:lang="it">
                    <pow xml:lang="it" />
                    <real xml:lang="it">10</real>
                    <real xml:lang="it">3</real>
                  </apply>
                </apply>
              </sequence>
            </apply>
          </define>
        </math>
      </region>
    </regions>
  </solveblock>
</region>

<region id="ID0EC23BB" top="873.6" left="9.6" height="249.600005435944" width="608.800014591217" xml:lang="it" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white" origin-positioning="true" xml:lang="it" show-axis-expressions="true">
    <xyPlot xml:lang="it">
      <title class="- topic/title " xml:lang="it" toc="yes" wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend xml:lang="it" />
      <traces xml:lang="it">
        <trace resultRef="14" xml:lang="it">
          <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="Solid" xml:lang="it">lines</traceStyle>
        </trace>
        <trace resultRef="15" xml:lang="it">
          <traceStyle color="#FF008000" symbol="none" line-weight="1" line-style="Solid" xml:lang="it">lines</traceStyle>
        </trace>
        <trace resultRef="16" xml:lang="it">
          <traceStyle color="#FF0000FF" symbol="none" line-weight="1" line-style="Solid" xml:lang="it">lines</traceStyle>
        </trace>
      </traces>
      <axes xml:lang="it">
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="10" xml:lang="it">
          <axisLine position="origin" positionticmark="4" legendWidth="67.3866666666667" xml:lang="it" />
          <axisGrid xml:lang="it">
            <gridFrequency xml:lang="it">11</gridFrequency>
            <gridLabels display="true" xml:lang="it" />
            <gridLines xml:lang="it" />
            <tickMarks display="true" xml:lang="it" />
          </axisGrid>
          <axisLabel xml:lang="it" />
          <markers xml:lang="it" />
          <numberFormat xml:lang="it">
            <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" xml:lang="it" />
          </numberFormat>
          <plotEquations xml:lang="it">
            <plotEquation xml:lang="it">
              <math xml:lang="it">
                <id xml:space="preserve" xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">z1</id>
              </math>
              <math xml:lang="it">
                <placeholder xml:lang="it" xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true" xml:lang="it">
            <startValue resultRef="17" xml:lang="it">
              <real xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">0</real>
            </startValue>
            <endValue xml:lang="it" resultRef="18">
              <id xml:space="preserve" xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">T1</id>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-40" end="60" xml:lang="it">
          <axisLine position="origin" positionticmark="0" legendWidth="98.1933333333333" xml:lang="it" />
          <axisGrid xml:lang="it">
            <gridFrequency xml:lang="it">11</gridFrequency>
            <gridLabels display="true" xml:lang="it" />
            <gridLines xml:lang="it" />
            <tickMarks display="true" xml:lang="it" />
          </axisGrid>
          <axisLabel xml:lang="it" />
          <markers xml:lang="it" />
          <numberFormat xml:lang="it">
            <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" xml:lang="it" />
          </numberFormat>
          <plotEquations xml:lang="it">
            <plotEquation xml:lang="it">
              <math xml:lang="it">
                <apply xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
                  <id xml:space="preserve" xml:lang="it">f</id>
                  <id xml:space="preserve" xml:lang="it">z1</id>
                </apply>
              </math>
              <math xml:lang="it">
                <placeholder xml:lang="it" xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation xml:lang="it">
              <math xml:lang="it">
                <apply xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
                  <id xml:space="preserve" xml:lang="it">g</id>
                  <id xml:space="preserve" xml:lang="it">z1</id>
                </apply>
              </math>
              <math xml:lang="it">
                <placeholder xml:lang="it" xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation xml:lang="it">
              <math xml:lang="it">
                <apply xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
                  <id xml:space="preserve" xml:lang="it">h</id>
                  <id xml:space="preserve" xml:lang="it">z1</id>
                </apply>
              </math>
              <math xml:lang="it">
                <placeholder xml:lang="it" xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true" xml:lang="it">
            <startValue resultRef="19" xml:lang="it">
              <real xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">-40</real>
            </startValue>
            <endValue xml:lang="it" resultRef="20">
              <real xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">60</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

<region id="ID0ERF4BB" top="1804.8" left="48" height="25.6" width="244.413333333333" xml:lang="it" xmlns="http://schemas.mathsoft.com/worksheet50">
  <math xml:lang="it">
    <define xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">
      <id xml:space="preserve" xml:lang="it" label="Variable">S1</id>
      <apply xml:lang="it">
        <id xml:space="preserve" xml:lang="it" label-is-contextual="true" label="Function">Rkadapt</id>
        <sequence xml:lang="it">
          <id xml:space="preserve" xml:lang="it" label-is-contextual="true" label="Variable">Y0</id>
          <id xml:space="preserve" xml:lang="it" label-is-contextual="true" label="Variable">z0</id>
          <id xml:space="preserve" xml:lang="it" label-is-contextual="true" label="Variable">T1</id>
          <id xml:space="preserve" xml:lang="it" label-is-contextual="true" label="Variable">num</id>
          <id xml:space="preserve" xml:lang="it" label-is-contextual="true" label="Variable">D</id>
        </sequence>
      </apply>
    </define>
  </math>
</region>

È inoltre possibile utilizzare le funzioni rkfixed, Bulstoer o Radau.

<region id="ID0EGP4BB" top="2246.4" left="28.8" height="249.600005435944" width="589.800012588501" xml:lang="it" xmlns="http://schemas.mathsoft.com/worksheet50">
  <plot background-type="white" origin-positioning="true" xml:lang="it" show-axis-expressions="true">
    <xyPlot xml:lang="it">
      <title class="- topic/title " xml:lang="it" toc="yes" wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
      <legend xml:lang="it" />
      <traces xml:lang="it">
        <trace resultRef="21" xml:lang="it">
          <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="Solid" xml:lang="it">lines</traceStyle>
        </trace>
        <trace resultRef="22" xml:lang="it">
          <traceStyle color="#FF008000" symbol="none" line-weight="1" line-style="Solid" xml:lang="it">lines</traceStyle>
        </trace>
        <trace resultRef="23" xml:lang="it">
          <traceStyle color="#FF0000FF" symbol="none" line-weight="1" line-style="Solid" xml:lang="it">lines</traceStyle>
        </trace>
      </traces>
      <axes xml:lang="it">
        <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="10" xml:lang="it">
          <axisLine position="origin" positionticmark="4" legendWidth="61.21" xml:lang="it" />
          <axisGrid xml:lang="it">
            <gridFrequency xml:lang="it">11</gridFrequency>
            <gridLabels display="true" xml:lang="it" />
            <gridLines xml:lang="it" />
            <tickMarks display="true" xml:lang="it" />
          </axisGrid>
          <axisLabel xml:lang="it" />
          <markers xml:lang="it" />
          <numberFormat xml:lang="it">
            <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" xml:lang="it" />
          </numberFormat>
          <plotEquations xml:lang="it">
            <plotEquation xml:lang="it">
              <math xml:lang="it">
                <id xml:space="preserve" xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">u</id>
              </math>
              <math xml:lang="it">
                <placeholder xml:lang="it" xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true" xml:lang="it">
            <startValue resultRef="24" xml:lang="it">
              <id xml:space="preserve" xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">z0</id>
            </startValue>
            <endValue xml:lang="it" resultRef="25">
              <id xml:space="preserve" xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">T1</id>
            </endValue>
          </xyDomain>
        </xAxis>
        <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-40" end="60" xml:lang="it">
          <axisLine position="origin" positionticmark="0" legendWidth="68.72" xml:lang="it" />
          <axisGrid xml:lang="it">
            <gridFrequency xml:lang="it">11</gridFrequency>
            <gridLabels display="true" xml:lang="it" />
            <gridLines xml:lang="it" />
            <tickMarks display="true" xml:lang="it" />
          </axisGrid>
          <axisLabel xml:lang="it" />
          <markers xml:lang="it" />
          <numberFormat xml:lang="it">
            <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" xml:lang="it" />
          </numberFormat>
          <plotEquations xml:lang="it">
            <plotEquation xml:lang="it">
              <math xml:lang="it">
                <id xml:space="preserve" xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">y0</id>
              </math>
              <math xml:lang="it">
                <placeholder xml:lang="it" xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation xml:lang="it">
              <math xml:lang="it">
                <id xml:space="preserve" xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">y1</id>
              </math>
              <math xml:lang="it">
                <placeholder xml:lang="it" xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
            <plotEquation xml:lang="it">
              <math xml:lang="it">
                <id xml:space="preserve" xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">y2</id>
              </math>
              <math xml:lang="it">
                <placeholder xml:lang="it" xmlns="http://schemas.mathsoft.com/math50" />
              </math>
            </plotEquation>
          </plotEquations>
          <xyDomain scale-type="linear" auto-scale="true" xml:lang="it">
            <startValue resultRef="26" xml:lang="it">
              <real xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">-33.608</real>
            </startValue>
            <endValue xml:lang="it" resultRef="27">
              <real xml:lang="it" xmlns="http://schemas.mathsoft.com/math50">60</real>
            </endValue>
          </xyDomain>
        </yAxis>
      </axes>
    </xyPlot>
  </plot>
</region>

Il grafico delle soluzioni di odesolve è quasi identico a quello delle soluzioni di Rkadapt.