From 7ebb7368fa53009acc87aaa8211fabee5608265b Mon Sep 17 00:00:00 2001 From: Karl Berry Date: Thu, 16 Apr 2009 16:08:57 +0000 Subject: mfpic4ode update (15apr09) git-svn-id: svn://tug.org/texlive/trunk@12734 c570f23f-e606-0410-a88d-b1316a301751 --- .../source/latex/mfpic4ode/mfpic4ode.dtx | 180 +++++++++++---------- 1 file changed, 99 insertions(+), 81 deletions(-) (limited to 'Master/texmf-dist/source/latex/mfpic4ode/mfpic4ode.dtx') diff --git a/Master/texmf-dist/source/latex/mfpic4ode/mfpic4ode.dtx b/Master/texmf-dist/source/latex/mfpic4ode/mfpic4ode.dtx index 0f712022f44..82b78640c93 100644 --- a/Master/texmf-dist/source/latex/mfpic4ode/mfpic4ode.dtx +++ b/Master/texmf-dist/source/latex/mfpic4ode/mfpic4ode.dtx @@ -22,7 +22,7 @@ %\NeedsTeXFormat{LaTeX2e}[1999/12/01] %\ProvidesPackage{mfpic4ode} %<*sty> - [2008/01/03 v0.2 mfpic4ode.dtx file] + [2009/04/15 v0.3 mfpic4ode.dtx file] % % %<*driver> @@ -38,7 +38,7 @@ % % \fi % -% \CheckSum{175} +% \CheckSum{185} % % \CharacterTable % {Upper-case \A\B\C\D\E\F\G\H\I\J\K\L\M\N\O\P\Q\R\S\T\U\V\W\X\Y\Z @@ -58,13 +58,14 @@ % % % \changes{v0.2}{2008/01/03}{First public version} +% \changes{v0.3}{2009/04/15}{Added connect environment to paths} % % \GetFileInfo{mfpic4ode.dtx} % % \DoNotIndex{\newcommand,\newenvironment} % -% \def\fileversion{0.2} -% \def\filedate{2008/01/03} +% \def\fileversion{0.3} +% \def\filedate{2009/04/15} % % \title{The \textsf{mfpic4ode} package\thanks{This document % corresponds to \textsf{mfpic4ode}~\fileversion, dated \filedate.}} @@ -164,13 +165,15 @@ % \end{verbatim} % The integral curve is drawn from short linear parts using |\ODEline| % command which expands to |\lines| command from |mfpic| package by -% default. A simple test is used to keep the arithmetics in reasonable -% bounds: if after the step the curve leaves the horizontal strip -% between |yneg| and |ypos| variables, then the evaluation is stopped -% (in fact, in this case we do not change the independent variable and -% we do the remaining steps with the same last point). Recall that -% |yneg| and |ypos| variables are set when you call |mfpicture| -% environment. If you call the environment as follows +% default. These linear parts are connected in connect environment and +% this allows to use prefixes like |\dotted| or |\dashed| to the +% trajectories. A simple test is used to keep the arithmetics in +% reasonable bounds: if after the step the curve leaves the horizontal +% strip between |yneg| and |ypos| variables, then the evaluation is +% stopped (in fact, in this case we do not change the independent +% variable and we do the remaining steps with the same last point). +% Recall that |yneg| and |ypos| variables are set when you call +% |mfpicture| environment. If you call the environment as follows % \begin{verbatim} % \begin{mfpic}[5][3]{-0.1}{1.5}{-0.1}{0.5} % ........... @@ -247,51 +250,60 @@ %%% ODEstep and the number of steps is ODEstepcount. The points are %%% stored in metapost variables x1,y1. \def\trajectory#1#2{ - \mfsrc{x1:=#1;y1:=#2; - for i=1 upto ODEstepcount: - x2:=x1+ODEstep; - y2:=y1+ODEstep*ODErhs(x1,y1);} - \ODEline{z1}{z2} - \mfsrc{ - if ((y2>yneg) and (y2yneg) and (y2yneg) and (y2yneg) and (y2yneg) and (y2yneg) and (y2yneg) and (y2yneg) and (y2yneg) and (y2xneg)): x1:=x2; y1:=y2 fi; - endfor - }} + \begin{connect} + \mfsrc{x1:=#1;y1:=#2; + TIMEsteps:=abs(TIMEend/TIMEstep); + TIME:=0; + for i=1 upto TIMEsteps: + k1:=ASf(x1,y1); + l1:=ASg(x1,y1); + k2:=ASf(x1+(TIMEstep*k1/2),y1+(TIMEstep*l1/2)); + l2:=ASg(x1+(TIMEstep*k1/2),y1+(TIMEstep*l1/2)); + k3:=ASf(x1+(TIMEstep*k2/2),y1+(TIMEstep*l2/2)); + l3:=ASg(x1+(TIMEstep*k2/2),y1+(TIMEstep*l2/2)); + k4:=ASf(x1+(TIMEstep*k3),y1+(TIMEstep*l3)); + l4:=ASg(x1+(TIMEstep*k3),y1+(TIMEstep*l3)); + k5:=((k1)/6)+((k2)/3)+((k3)/3)+((k4)/6); + l5:=(l1/6)+(l2/3)+(l3/3)+(l4/6); + x2:=x1+(TIMEstep*k5); + y2:=y1+(TIMEstep*l5);} + \ODEline{z1}{z2} + \mfsrc{ + if ((y2>yneg) and (y2xneg)): x1:=x2; y1:=y2 fi; + endfor + } + \end{connect} +} \catcode`\@12\relax % -- cgit v1.2.3