Exemple : Fonction jacobienne augmentée pour Stiffb et Stiffr

Fonctions > Résolution et optimisation > Solveurs d'équations différentielles > Exemple : Fonction jacobienne augmentée pour Stiffb et Stiffr

Utilisez la fonction Jacob pour calculer la matrice jacobienne augmentée pour une équation différentielle ordinaire (EOD), puis l'indiquer comme entrée aux solveurs Stiffb et Stiffr.

1. Définissez un système pour vos 4 inconnues :

Les variables yi sont des fonctions de x.

2. Définissez les constantes dans le système.

3. Définissez les valeurs initiales.

4. Définissez un D(x,y) de fonction vectorielle correspondant au côté droit du système.

5. Créez la matrice jacobienne augmentée à l'aide de la fonction augment :

<region id="ID0EWVU5" 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. Appelez les fonctions Stiffb et Stiffr :

<region id="ID0EYYU5" 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="ID0ENZU5" 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="ID0EC1U5" 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="ID0E41U5" 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="ID0E52U5" 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>

Les matrices renvoyées contiennent 5 colonnes correspondant au nombre de points, ainsi que les solutions pour les 4 inconnues.

7. Extrayez les solutions pour les 4 inconnues à partir des matrices Stiffb et Stiffr renvoyées :

<region id="ID0E46U5" 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="ID0E4AV5" 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="ID0ESBV5" 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="ID0ELCV5" 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="ID0EADV5" 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="ID0EZDV5" 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="ID0EOEV5" 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="ID0EGFV5" 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="ID0E1FV5" 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. Tracez et comparez les solutions renvoyées pour chaque inconnue à partir des deux fonctions :

<region id="ID0EOHV5" 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="ID0E3HV5" 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="ID0EKIV5" 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="ID0EYIV5" 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>

Les tracés montrent que les deux solveurs EDO renvoient des solutions identiques.

Rubriques associées

Solveurs d'EDO

Est-ce que cela a été utile ?