Example: Augmented Jacobian for Stiffb and Stiffr

Functions > Solving and Optimization > Differential Equation Solvers > Example: Augmented Jacobian for Stiffb and Stiffr

Use the Jacob function to compute the augmented Jacobian matrix for an ordinary differential equation (ODE), and then provide it as input to solvers Stiffb and Stiffr.

1. Define a system of four unknowns:

The yi variables are functions of x.

2. Define the constants in the system.

3. Define the initial values.

4. Define a vector function D(x,y) corresponding to the right-hand side of the system.

5. Use the augment function to create the augmented Jacobian:

<region id="ID0EH2U5" top="883.20000000000016" left="38.400000000000006" height="130.50400000000002" width="912.82666666666671" xmlns="http://schemas.mathsoft.com/worksheet50">
 <math resultRef="162">
 <define xmlns="http://schemas.mathsoft.com/math50">
 <function>
 <id xml:space="preserve" labels="VARIABLE">AJ</id>
 <boundVars>
 <id xml:space="preserve" labels="VARIABLE">x</id>
 <id xml:space="preserve" labels="VARIABLE">y</id>
 </boundVars>
 </function>
 <symEval>
 <apply>
 <id xml:space="preserve" labels="FUNCTION" label-is-contextual="true">augment</id>
 <sequence>
 <apply>
 <id xml:space="preserve">Jacob</id>
 <sequence>
 <apply>
 <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">D</id>
 <sequence>
 <id xml:space="preserve">x</id>
 <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">y</id>
 </sequence>
 </apply>
 <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">x</id>
 </sequence>
 </apply>
 <apply>
 <id xml:space="preserve">Jacob</id>
 <sequence>
 <apply>
 <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">D</id>
 <sequence>
 <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">x</id>
 <id xml:space="preserve">y</id>
 </sequence>
 </apply>
 <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">y</id>
 </sequence>
 </apply>
 </sequence>
 </apply>
 <command>
 <placeholder />
 </command>
 <symResult>
 <matrix rows="4" cols="5">
 <real>0</real>
 <real>0</real>
 <real>0</real>
 <real>0</real>
 <apply>
 <minus />
 <apply>
 <mult />
 <real>-1.34</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>2</real>
 </apply>
 </apply>
 <apply>
 <mult />
 <real>8000</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>1</real>
 </apply>
 </apply>
 </apply>
 <apply>
 <plus />
 <apply>
 <mult />
 <real>1.34</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>2</real>
 </apply>
 </apply>
 <apply>
 <mult />
 <real>8000</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>1</real>
 </apply>
 </apply>
 </apply>
 <apply>
 <mult />
 <real>-1.34</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>2</real>
 </apply>
 </apply>
 <apply>
 <mult />
 <real>8000</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>1</real>
 </apply>
 </apply>
 <apply>
 <mult />
 <real>-8000</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>0</real>
 </apply>
 </apply>
 <apply>
 <minus />
 <apply>
 <plus />
 <apply>
 <mult />
 <apply>
 <neg />
 <parens>
 <apply>
 <mult />
 <real>1.6</real>
 <apply>
 <pow />
 <real>10</real>
 <real>9</real>
 </apply>
 </apply>
 </parens>
 </apply>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>2</real>
 </apply>
 </apply>
 <apply>
 <mult />
 <real>8000</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>0</real>
 </apply>
 </apply>
 </apply>
 <apply>
 <mult />
 <real>160000000</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>1</real>
 </apply>
 </apply>
 </apply>
 <apply>
 <mult />
 <apply>
 <neg />
 <parens>
 <apply>
 <mult />
 <real>1.6</real>
 <apply>
 <pow />
 <real>10</real>
 <real>9</real>
 </apply>
 </apply>
 </parens>
 </apply>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>2</real>
 </apply>
 </apply>
 <apply>
 <mult />
 <real>8000</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>0</real>
 </apply>
 </apply>
 <apply>
 <mult />
 <real>-1.34</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>0</real>
 </apply>
 </apply>
 <apply>
 <plus />
 <apply>
 <mult />
 <real>1.34</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>0</real>
 </apply>
 </apply>
 <apply>
 <mult />
 <apply>
 <neg />
 <parens>
 <apply>
 <mult />
 <real>1.6</real>
 <apply>
 <pow />
 <real>10</real>
 <real>9</real>
 </apply>
 </apply>
 </parens>
 </apply>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>1</real>
 </apply>
 </apply>
 </apply>
 <apply>
 <plus />
 <apply>
 <mult />
 <real>-1.34</real>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>0</real>
 </apply>
 </apply>
 <apply>
 <mult />
 <apply>
 <neg />
 <parens>
 <apply>
 <mult />
 <real>1.6</real>
 <apply>
 <pow />
 <real>10</real>
 <real>9</real>
 </apply>
 </apply>
 </parens>
 </apply>
 <apply>
 <indexer />
 <id xml:space="preserve" labels="*">y</id>
 <real>1</real>
 </apply>
 </apply>
 </apply>
 <real>0</real>
 <real>0</real>
 <real>0</real>
 <real>1</real>
 <real>-1</real>
 </matrix>
 </symResult>
 </symEval>
 </define>
 </math>
</region>

6. Call the Stiffb and Stiffr functions:

<region id="ID0EK5U5" top="1017.6000000000001" left="38.400000000000006" height="25.6" width="47.876666666666665" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="164"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE">r1</id> <real>0</real> </define> </math> </region>	<region id="ID0E65U5" top="1056.0000000000002" left="38.400000000000006" height="25.6" width="55.856666666666669" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="166"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE">r2</id> <real>20</real> </define> </math> </region>	<region id="ID0EU6U5" top="1094.4" left="38.400000000000006" height="25.6" width="100.86666666666667" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="168"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE">npoints</id> <real>150</real> </define> </math> </region>
<region id="ID0EPAV5" top="1132.8000000000002" left="38.400000000000006" height="27.791200000000003" width="275.35666666666668" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="170"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE">Yb</id> <apply> <id xml:space="preserve" labels="FUNCTION" label-is-contextual="true">Stiffb</id> <sequence> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">y<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">init</Subscript></Span> </id> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">r1</id> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">r2</id> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">npoints</id> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">D</id> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">AJ</id> </sequence> </apply> </define> </math> </region>
<region id="ID0EQBV5" top="1171.2000000000003" left="38.400000000000006" height="27.791200000000003" width="272.8633333333334" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="172"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE">Yr</id> <apply> <id xml:space="preserve" labels="FUNCTION" label-is-contextual="true">Stiffr</id> <sequence> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">y<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">init</Subscript></Span> </id> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">r1</id> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">r2</id> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">npoints</id> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">D</id> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">AJ</id> </sequence> </apply> </define> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" /> </resultFormat> </math> </region>

The returned matrices contain 5 columns corresponding to the number of points and the solutions for the four unknowns.

7. Extract the solutions for the four unknowns from the returned Stiffb and Stiffr matrices:

<region id="ID0ERFV5" top="1305.6000000000001" left="38.400000000000006" height="27.752000000000002" width="66.51666666666668" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="174"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE">n</id> <apply> <matcol /> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">Yb</id> <real>0</real> </apply> </define> </math> </region>
<region id="ID0ETGV5" top="1344.0000000000002" left="38.400000000000006" height="29.943200000000004" width="78.453333333333333" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="176"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">yb<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> <apply> <matcol /> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">Yb</id> <real>1</real> </apply> </define> </math> </region>	<region id="ID0EIHV5" top="1344.0000000000002" left="153.60000000000002" height="29.943200000000004" width="77.720000000000013" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="178"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">yr<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript></Span> </id> <apply> <matcol /> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">Yr</id> <real>1</real> </apply> </define> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" /> </resultFormat> </math> </region>
<region id="ID0EBIV5" top="1382.4" left="38.400000000000006" height="29.943200000000004" width="78.453333333333333" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="180"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">yb<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> <apply> <matcol /> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">Yb</id> <real>2</real> </apply> </define> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" /> </resultFormat> </math> </region>	<region id="ID0EWIV5" top="1382.4" left="153.60000000000002" height="29.943200000000004" width="77.720000000000013" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="182"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">yr<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript></Span> </id> <apply> <matcol /> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">Yr</id> <real>2</real> </apply> </define> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" /> </resultFormat> </math> </region>
<region id="ID0EPJV5" top="1420.8000000000002" left="38.400000000000006" height="29.943200000000004" width="78.453333333333333" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="184"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">yb<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">2</Subscript></Span> </id> <apply> <matcol /> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">Yb</id> <real>3</real> </apply> </define> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" /> </resultFormat> </math> </region>	<region id="ID0EEKV5" top="1420.8000000000002" left="153.60000000000002" height="29.943200000000004" width="77.720000000000013" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="186"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">yr<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">2</Subscript></Span> </id> <apply> <matcol /> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">Yr</id> <real>3</real> </apply> </define> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" /> </resultFormat> </math> </region>
<region id="ID0E3KV5" top="1459.2000000000003" left="38.400000000000006" height="29.943200000000004" width="78.453333333333333" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="188"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">yb<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">3</Subscript></Span> </id> <apply> <matcol /> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">Yb</id> <real>4</real> </apply> </define> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" /> </resultFormat> </math> </region>	<region id="ID0EQLV5" top="1459.2000000000003" left="153.60000000000002" height="29.943200000000004" width="77.720000000000013" xmlns="http://schemas.mathsoft.com/worksheet50"> <math resultRef="190"> <define xmlns="http://schemas.mathsoft.com/math50"> <id xml:space="preserve" labels="VARIABLE"> <Span xmlns="http://schemas.microsoft.com/winfx/2006/xaml/presentation">yr<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">3</Subscript></Span> </id> <apply> <matcol /> <id xml:space="preserve" labels="VARIABLE" label-is-contextual="true">Yr</id> <real>4</real> </apply> </define> <resultFormat> <matrix size="12,12" offset="0,0" show-indices="false" expand-nested-arrays="true" /> </resultFormat> </math> </region>

8. Plot and compare the returned solutions for each unknown from the two functions:

<region id="ID0EENV5" top="1516.8000000000002" left="38.400000000000006" height="211.200004863739" width="672.600008368492" xmlns="http://schemas.mathsoft.com/worksheet50">
 <plot background-type="white" origin-positioning="true">
 <xyPlot>
 <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
 <legend />
 <traces>
 <trace resultRef="192">
 <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
 </trace>
 <trace resultRef="193">
 <traceStyle color="#FF0000FF" symbol="none" line-weight="1" line-style="DashDotDot">lines</traceStyle>
 </trace>
 </traces>
 <graph-size width="572.600008368492" height="153.600004863739" />
 <axes>
 <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="20">
 <axisLine position="origin" positionticmark="5" legendWidth="61.7533333333333" />
 <axisGrid>
 <gridFrequency>5</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="196">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">n</id>
 </math>
 <math resultRef="197">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </startValue>
 <secondTickValue resultRef="195">
 <real xmlns="http://schemas.mathsoft.com/math50">5</real>
 </secondTickValue>
 <endValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </endValue>
 </xyDomain>
 </xAxis>
 <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-1.25" end="0.25">
 <axisLine position="origin" positionticmark="0" legendWidth="73.69" />
 <axisGrid>
 <gridFrequency>7</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="204">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">
 yb<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript>
 </id>
 </math>
 <math resultRef="205">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 <plotEquation>
 <math resultRef="206">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">
 yr<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">0</Subscript>
 </id>
 </math>
 <math resultRef="207">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue resultRef="199">
 <real xmlns="http://schemas.mathsoft.com/math50">-1.25</real>
 </startValue>
 <secondTickValue resultRef="201">
 <real xmlns="http://schemas.mathsoft.com/math50">-1</real>
 </secondTickValue>
 <endValue resultRef="203">
 <real xmlns="http://schemas.mathsoft.com/math50">0.25</real>
 </endValue>
 </xyDomain>
 </yAxis>
 </axes>
 </xyPlot>
 </plot>
</region>

<region id="ID0ESNV5" top="1747.2000000000003" left="38.400000000000006" height="230.400004863739" width="672.60000836849213" xmlns="http://schemas.mathsoft.com/worksheet50">
 <plot background-type="white" origin-positioning="true">
 <xyPlot>
 <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
 <legend />
 <traces>
 <trace resultRef="218">
 <traceStyle color="#FFED1D2F" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
 </trace>
 <trace resultRef="219">
 <traceStyle color="#FF2E3192" symbol="none" line-weight="1" line-style="DashDotDot">lines</traceStyle>
 </trace>
 </traces>
 <graph-size width="572.600008368492" height="172.800004863739" />
 <axes>
 <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="20">
 <axisLine position="origin" positionticmark="4" legendWidth="61.7533333333333" />
 <axisGrid>
 <gridFrequency>5</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="222">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">n</id>
 </math>
 <math resultRef="223">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </startValue>
 <secondTickValue resultRef="221">
 <real xmlns="http://schemas.mathsoft.com/math50">5</real>
 </secondTickValue>
 <endValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </endValue>
 </xyDomain>
 </xAxis>
 <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-1E-09" end="2.06795153138257E-25">
 <axisLine position="ticknumberlock" positionticmark="0" legendWidth="73.69" />
 <axisGrid>
 <gridFrequency>5</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="230">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">
 yb<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript>
 </id>
 </math>
 <math resultRef="231">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 <plotEquation>
 <math resultRef="232">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">
 yr<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">1</Subscript>
 </id>
 </math>
 <math resultRef="233">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue resultRef="225">
 <apply xmlns="http://schemas.mathsoft.com/math50">
 <neg />
 <apply>
 <pow />
 <real>10</real>
 <real>-9</real>
 </apply>
 </apply>
 </startValue>
 <secondTickValue resultRef="227">
 <apply xmlns="http://schemas.mathsoft.com/math50">
 <mult />
 <real>-7.5</real>
 <apply>
 <pow />
 <real>10</real>
 <real>-10</real>
 </apply>
 </apply>
 </secondTickValue>
 <endValue resultRef="229">
 <apply xmlns="http://schemas.mathsoft.com/math50">
 <mult />
 <real>-1</real>
 <apply>
 <pow />
 <real>10</real>
 <real>-11</real>
 </apply>
 </apply>
 </endValue>
 </xyDomain>
 </yAxis>
 </axes>
 </xyPlot>
 </plot>
</region>

<region id="ID0EAOV5" top="1996.8000000000002" left="38.400000000000006" height="211.200004863739" width="672.600008368492" xmlns="http://schemas.mathsoft.com/worksheet50">
 <plot background-type="white" origin-positioning="true">
 <xyPlot>
 <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
 <legend />
 <traces>
 <trace resultRef="244">
 <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
 </trace>
 <trace resultRef="245">
 <traceStyle color="#FF0000FF" symbol="none" line-weight="1" line-style="DashDotDot">lines</traceStyle>
 </trace>
 </traces>
 <graph-size width="572.600008368492" height="153.600004863739" />
 <axes>
 <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="20">
 <axisLine position="origin" positionticmark="1" legendWidth="61.7533333333333" />
 <axisGrid>
 <gridFrequency>5</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="248">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">n</id>
 </math>
 <math resultRef="249">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </startValue>
 <secondTickValue resultRef="247">
 <real xmlns="http://schemas.mathsoft.com/math50">5</real>
 </secondTickValue>
 <endValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </endValue>
 </xyDomain>
 </xAxis>
 <yAxis rank="1" legend-position="PlotBoundaryLeft" start="-0.5" end="2.5">
 <axisLine position="origin" positionticmark="0" legendWidth="73.69" />
 <axisGrid>
 <gridFrequency>7</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="256">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">
 yb<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">2</Subscript>
 </id>
 </math>
 <math resultRef="257">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 <plotEquation>
 <math resultRef="258">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">
 yr<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">2</Subscript>
 </id>
 </math>
 <math resultRef="259">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue resultRef="251">
 <real xmlns="http://schemas.mathsoft.com/math50">-0.5</real>
 </startValue>
 <secondTickValue resultRef="253">
 <real xmlns="http://schemas.mathsoft.com/math50">0</real>
 </secondTickValue>
 <endValue resultRef="255">
 <real xmlns="http://schemas.mathsoft.com/math50">2.5</real>
 </endValue>
 </xyDomain>
 </yAxis>
 </axes>
 </xyPlot>
 </plot>
</region>

<region id="ID0EOOV5" top="2227.2000000000003" left="38.400000000000006" height="230.40000486373904" width="671.600008368492" xmlns="http://schemas.mathsoft.com/worksheet50">
 <plot background-type="white" origin-positioning="true">
 <xyPlot>
 <title class="- topic/title " wwtype:type="Paragraph" xmlns:wwtype="urn:WebWorks-Type-Schema" />
 <legend />
 <traces>
 <trace resultRef="270">
 <traceStyle color="#FFFF0000" symbol="none" line-weight="1" line-style="Solid">lines</traceStyle>
 </trace>
 <trace resultRef="271">
 <traceStyle color="#FF000000" symbol="none" line-weight="1" line-style="DashDotDot">lines</traceStyle>
 </trace>
 </traces>
 <graph-size width="571.600008368492" height="172.800004863739" />
 <axes>
 <xAxis rank="1" legend-position="PlotBoundaryBottom" start="0" end="20">
 <axisLine position="ticknumberlock" positionticmark="0" legendWidth="61.7533333333333" />
 <axisGrid>
 <gridFrequency>5</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="274">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">n</id>
 </math>
 <math resultRef="275">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </startValue>
 <secondTickValue resultRef="273">
 <real xmlns="http://schemas.mathsoft.com/math50">5</real>
 </secondTickValue>
 <endValue>
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </endValue>
 </xyDomain>
 </xAxis>
 <yAxis rank="1" legend-position="PlotBoundaryLeft" start="0" end="8E-06">
 <axisLine position="ticknumberlock" positionticmark="0" legendWidth="73.69" />
 <axisGrid>
 <gridFrequency>3</gridFrequency>
 <gridLabels display="true" />
 <gridLines />
 <tickMarks display="true" />
 </axisGrid>
 <axisLabel />
 <markers />
 <numberFormat>
 <general precision="3" show-trailing-zeros="false" radix="dec" zero-threshold="15" imaginary-value="i" exponential-threshold="3" />
 </numberFormat>
 <plotEquations>
 <plotEquation>
 <math resultRef="282">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">
 yb<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">3</Subscript>
 </id>
 </math>
 <math resultRef="283">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 <plotEquation>
 <math resultRef="284">
 <id xml:space="preserve" xmlns="http://schemas.mathsoft.com/math50">
 yr<Subscript xmlns="clr-namespace:Ptc.Wpf;assembly=Ptc.Core">3</Subscript>
 </id>
 </math>
 <math resultRef="285">
 <placeholder xmlns="http://schemas.mathsoft.com/math50" />
 </math>
 </plotEquation>
 </plotEquations>
 <xyDomain scale-type="linear" auto-scale="true">
 <startValue resultRef="277">
 <real xmlns="http://schemas.mathsoft.com/math50">0</real>
 </startValue>
 <secondTickValue resultRef="279">
 <apply xmlns="http://schemas.mathsoft.com/math50">
 <mult />
 <real>4</real>
 <apply>
 <pow />
 <real>10</real>
 <real>-6</real>
 </apply>
 </apply>
 </secondTickValue>
 <endValue resultRef="281">
 <apply xmlns="http://schemas.mathsoft.com/math50">
 <mult />
 <real>8</real>
 <apply>
 <pow />
 <real>10</real>
 <real>-6</real>
 </apply>
 </apply>
 </endValue>
 </xyDomain>
 </yAxis>
 </axes>
 </xyPlot>
 </plot>
</region>

The plots show that the two ODE solvers return identical solutions.