--- a/main/odepkg/doc/mfunref.texi	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/doc/mfunref.texi	Wed Feb 08 21:02:51 2012 +0000
@@ -29,11 +29,22 @@
 
 If this function is called with no return argument then plot the solution over time in a figure window while solving the set of DDEs that are defined in a function and specified by the function handle @var{@@fun}. The second input argument @var{slot} is a double vector that defines the time slot, @var{init} is a double vector that defines the initial values of the states, @var{lags} is a double vector that describes the lags of time, @var{hist} is a double matrix and describes the history of the DDEs, @var{opt} can optionally be a structure array that keeps the options created with the command @command{odeset} and @var{par1}, @var{par2}, @dots{} can optionally be other input arguments of any type that have to be passed to the function defined by @var{@@fun}.
 
+In other words, this function will solve a problem of the form
+@example
+dy/dt = fun (t, y(t), y(t-lags(1), y(t-lags(2), @dots{})))
+y(slot(1)) = init
+y(slot(1)-lags(1)) = hist(1), y(slot(1)-lags(2)) = hist(2), @dots{} 
+@end example
+
 If this function is called with one return argument then return the solution @var{sol} of type structure array after solving the set of DDEs. The solution @var{sol} has the fields @var{x} of type double column vector for the steps chosen by the solver, @var{y} of type double column vector for the solutions at each time step of @var{x}, @var{solver} of type string for the solver name and optionally the extended time stamp information @var{xe}, the extended solution information @var{ye} and the extended index information @var{ie} all of type double column vector that keep the informations of the event function if an event function handle is set in the option argument @var{opt}.
 
 If this function is called with more than one return argument then return the time stamps @var{t}, the solution values @var{y} and optionally the extended time stamp information @var{xe}, the extended solution information @var{ye} and the extended index information @var{ie} all of type double column vector.
 
-For example, solve an anonymous implementation of a chaotic behavior
+For example:
+@itemize @minus
+@item
+the following code solves an anonymous implementation of a chaotic behavior
+
 @example
 fcao = @@(vt, vy, vz) [2 * vz / (1 + vz^9.65) - vy];
 
@@ -43,6 +54,29 @@
 vlag = interp1 (vsol.x, vsol.y, vsol.x - 2);
 plot (vsol.y, vlag); legend ("fcao (t,y,z)");
 @end example
+
+@item
+to solve the following problem with two delayed state variables
+
+@example
+d y1(t)/dt = -y1(t)
+d y2(t)/dt = -y2(t) + y1(t-5)
+d y3(t)/dt = -y3(t) + y2(t-10)*y1(t-10)
+@end example
+
+one might do the following
+
+@example
+function f = fun (t, y, yd)
+f(1) = -y(1);                   %% y1' = -y1(t)
+f(2) = -y(2) + yd(1,1);         %% y2' = -y2(t) + y1(t-lags(1))
+f(3) = -y(3) + yd(2,2)*yd(1,2); %% y3' = -y3(t) + y2(t-lags(2))*y1(t-lags(2))
+endfunction
+T = [0,20]
+res = ode23d (@@fun, T, [1;1;1], [5, 10], ones (3,2));
+@end example
+
+@end itemize
 @end deftypefn
 
 @deftypefn  {Function File} {[@var{}] =} ode45 (@var{@@fun}, @var{slot}, @var{init}, [@var{opt}], [@var{par1}, @var{par2}, @dots{}])
@@ -76,20 +110,54 @@
 
 If this function is called with no return argument then plot the solution over time in a figure window while solving the set of DDEs that are defined in a function and specified by the function handle @var{@@fun}. The second input argument @var{slot} is a double vector that defines the time slot, @var{init} is a double vector that defines the initial values of the states, @var{lags} is a double vector that describes the lags of time, @var{hist} is a double matrix and describes the history of the DDEs, @var{opt} can optionally be a structure array that keeps the options created with the command @command{odeset} and @var{par1}, @var{par2}, @dots{} can optionally be other input arguments of any type that have to be passed to the function defined by @var{@@fun}.
 
+In other words, this function will solve a problem of the form
+@example
+dy/dt = fun (t, y(t), y(t-lags(1), y(t-lags(2), @dots{})))
+y(slot(1)) = init
+y(slot(1)-lags(1)) = hist(1), y(slot(1)-lags(2)) = hist(2), @dots{} 
+@end example
+
 If this function is called with one return argument then return the solution @var{sol} of type structure array after solving the set of DDEs. The solution @var{sol} has the fields @var{x} of type double column vector for the steps chosen by the solver, @var{y} of type double column vector for the solutions at each time step of @var{x}, @var{solver} of type string for the solver name and optionally the extended time stamp information @var{xe}, the extended solution information @var{ye} and the extended index information @var{ie} all of type double column vector that keep the informations of the event function if an event function handle is set in the option argument @var{opt}.
 
 If this function is called with more than one return argument then return the time stamps @var{t}, the solution values @var{y} and optionally the extended time stamp information @var{xe}, the extended solution information @var{ye} and the extended index information @var{ie} all of type double column vector.
 
-For example, solve an anonymous implementation of a chaotic behavior
+For example:
+@itemize @minus
+@item
+the following code solves an anonymous implementation of a chaotic behavior
+
 @example
 fcao = @@(vt, vy, vz) [2 * vz / (1 + vz^9.65) - vy];
 
-vopt = odeset ("NormControl", "on", "RelTol", 1e-4);
+vopt = odeset ("NormControl", "on", "RelTol", 1e-3);
 vsol = ode45d (fcao, [0, 100], 0.5, 2, 0.5, vopt);
 
 vlag = interp1 (vsol.x, vsol.y, vsol.x - 2);
 plot (vsol.y, vlag); legend ("fcao (t,y,z)");
 @end example
+
+@item
+to solve the following problem with two delayed state variables
+
+@example
+d y1(t)/dt = -y1(t)
+d y2(t)/dt = -y2(t) + y1(t-5)
+d y3(t)/dt = -y3(t) + y2(t-10)*y1(t-10)
+@end example
+
+one might do the following
+
+@example
+function f = fun (t, y, yd)
+f(1) = -y(1);                   %% y1' = -y1(t)
+f(2) = -y(2) + yd(1,1);         %% y2' = -y2(t) + y1(t-lags(1))
+f(3) = -y(3) + yd(2,2)*yd(1,2); %% y3' = -y3(t) + y2(t-lags(2))*y1(t-lags(2))
+endfunction
+T = [0,20]
+res = ode45d (@@fun, T, [1;1;1], [5, 10], ones (3,2));
+@end example
+
+@end itemize
 @end deftypefn
 
 @deftypefn  {Function File} {[@var{}] =} ode54 (@var{@@fun}, @var{slot}, @var{init}, [@var{opt}], [@var{par1}, @var{par2}, @dots{}])
@@ -123,11 +191,22 @@
 
 If this function is called with no return argument then plot the solution over time in a figure window while solving the set of DDEs that are defined in a function and specified by the function handle @var{@@fun}. The second input argument @var{slot} is a double vector that defines the time slot, @var{init} is a double vector that defines the initial values of the states, @var{lags} is a double vector that describes the lags of time, @var{hist} is a double matrix and describes the history of the DDEs, @var{opt} can optionally be a structure array that keeps the options created with the command @command{odeset} and @var{par1}, @var{par2}, @dots{} can optionally be other input arguments of any type that have to be passed to the function defined by @var{@@fun}.
 
+In other words, this function will solve a problem of the form
+@example
+dy/dt = fun (t, y(t), y(t-lags(1), y(t-lags(2), @dots{})))
+y(slot(1)) = init
+y(slot(1)-lags(1)) = hist(1), y(slot(1)-lags(2)) = hist(2), @dots{} 
+@end example
+
 If this function is called with one return argument then return the solution @var{sol} of type structure array after solving the set of DDEs. The solution @var{sol} has the fields @var{x} of type double column vector for the steps chosen by the solver, @var{y} of type double column vector for the solutions at each time step of @var{x}, @var{solver} of type string for the solver name and optionally the extended time stamp information @var{xe}, the extended solution information @var{ye} and the extended index information @var{ie} all of type double column vector that keep the informations of the event function if an event function handle is set in the option argument @var{opt}.
 
 If this function is called with more than one return argument then return the time stamps @var{t}, the solution values @var{y} and optionally the extended time stamp information @var{xe}, the extended solution information @var{ye} and the extended index information @var{ie} all of type double column vector.
 
-For example, solve an anonymous implementation of a chaotic behavior
+For example:
+@itemize @minus
+@item
+the following code solves an anonymous implementation of a chaotic behavior
+
 @example
 fcao = @@(vt, vy, vz) [2 * vz / (1 + vz^9.65) - vy];
 
@@ -137,6 +216,29 @@
 vlag = interp1 (vsol.x, vsol.y, vsol.x - 2);
 plot (vsol.y, vlag); legend ("fcao (t,y,z)");
 @end example
+
+@item
+to solve the following problem with two delayed state variables
+
+@example
+d y1(t)/dt = -y1(t)
+d y2(t)/dt = -y2(t) + y1(t-5)
+d y3(t)/dt = -y3(t) + y2(t-10)*y1(t-10)
+@end example
+
+one might do the following
+
+@example
+function f = fun (t, y, yd)
+f(1) = -y(1);                   %% y1' = -y1(t)
+f(2) = -y(2) + yd(1,1);         %% y2' = -y2(t) + y1(t-lags(1))
+f(3) = -y(3) + yd(2,2)*yd(1,2); %% y3' = -y3(t) + y2(t-lags(2))*y1(t-lags(2))
+endfunction
+T = [0,20]
+res = ode54d (@@fun, T, [1;1;1], [5, 10], ones (3,2));
+@end example
+
+@end itemize
 @end deftypefn
 
 @deftypefn  {Function File} {[@var{}] =} ode78 (@var{@@fun}, @var{slot}, @var{init}, [@var{opt}], [@var{par1}, @var{par2}, @dots{}])
@@ -170,11 +272,22 @@
 
 If this function is called with no return argument then plot the solution over time in a figure window while solving the set of DDEs that are defined in a function and specified by the function handle @var{@@fun}. The second input argument @var{slot} is a double vector that defines the time slot, @var{init} is a double vector that defines the initial values of the states, @var{lags} is a double vector that describes the lags of time, @var{hist} is a double matrix and describes the history of the DDEs, @var{opt} can optionally be a structure array that keeps the options created with the command @command{odeset} and @var{par1}, @var{par2}, @dots{} can optionally be other input arguments of any type that have to be passed to the function defined by @var{@@fun}.
 
+In other words, this function will solve a problem of the form
+@example
+dy/dt = fun (t, y(t), y(t-lags(1), y(t-lags(2), @dots{})))
+y(slot(1)) = init
+y(slot(1)-lags(1)) = hist(1), y(slot(1)-lags(2)) = hist(2), @dots{} 
+@end example
+
 If this function is called with one return argument then return the solution @var{sol} of type structure array after solving the set of DDEs. The solution @var{sol} has the fields @var{x} of type double column vector for the steps chosen by the solver, @var{y} of type double column vector for the solutions at each time step of @var{x}, @var{solver} of type string for the solver name and optionally the extended time stamp information @var{xe}, the extended solution information @var{ye} and the extended index information @var{ie} all of type double column vector that keep the informations of the event function if an event function handle is set in the option argument @var{opt}.
 
 If this function is called with more than one return argument then return the time stamps @var{t}, the solution values @var{y} and optionally the extended time stamp information @var{xe}, the extended solution information @var{ye} and the extended index information @var{ie} all of type double column vector.
 
-For example, solve an anonymous implementation of a chaotic behavior
+For example:
+@itemize @minus
+@item
+the following code solves an anonymous implementation of a chaotic behavior
+
 @example
 fcao = @@(vt, vy, vz) [2 * vz / (1 + vz^9.65) - vy];
 
@@ -184,6 +297,29 @@
 vlag = interp1 (vsol.x, vsol.y, vsol.x - 2);
 plot (vsol.y, vlag); legend ("fcao (t,y,z)");
 @end example
+
+@item
+to solve the following problem with two delayed state variables
+
+@example
+d y1(t)/dt = -y1(t)
+d y2(t)/dt = -y2(t) + y1(t-5)
+d y3(t)/dt = -y3(t) + y2(t-10)*y1(t-10)
+@end example
+
+one might do the following
+
+@example
+function f = fun (t, y, yd)
+f(1) = -y(1);                   %% y1' = -y1(t)
+f(2) = -y(2) + yd(1,1);         %% y2' = -y2(t) + y1(t-lags(1))
+f(3) = -y(3) + yd(2,2)*yd(1,2); %% y3' = -y3(t) + y2(t-lags(2))*y1(t-lags(2))
+endfunction
+T = [0,20]
+res = ode78d (@@fun, T, [1;1;1], [5, 10], ones (3,2));
+@end example
+
+@end itemize
 @end deftypefn
 
 @deftypefn  {Function File} {[@var{}] =} odebwe (@var{@@fun}, @var{slot}, @var{init}, [@var{opt}], [@var{par1}, @var{par2}, @dots{}])
@@ -340,6 +476,8 @@
 @example
 demo odepkg_testsuite_calcmescd
 @end example
+
+This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
 @end deftypefn
 
 @deftypefn {Function File} {[@var{scd}] =} odepkg_testsuite_calcscd (@var{solution}, @var{reference}, @var{abstol}, @var{reltol})
@@ -350,6 +488,8 @@
 @example
 demo odepkg_testsuite_calcscd
 @end example
+
+This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
 @end deftypefn
 
 @deftypefn {Function File} {[@var{solution}] =} odepkg_testsuite_chemakzo (@var{@@solver}, @var{reltol})
@@ -360,6 +500,8 @@
 @example
 demo odepkg_testsuite_chemakzo
 @end example
+
+This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
 @end deftypefn
 
 @deftypefn {Function File} {[@var{solution}] =} odepkg_testsuite_hires (@var{@@solver}, @var{reltol})
@@ -370,6 +512,8 @@
 @example
 demo odepkg_testsuite_hires
 @end example
+
+This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
 @end deftypefn
 
 @deftypefn {Function File} {[@var{solution}] =} odepkg_testsuite_implakzo (@var{@@solver}, @var{reltol})
@@ -380,6 +524,8 @@
 @example
 demo odepkg_testsuite_implakzo
 @end example
+
+This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
 @end deftypefn
 
 @deftypefn {Function File} {[@var{solution}] =} odepkg_testsuite_implrober (@var{@@solver}, @var{reltol})
@@ -390,6 +536,8 @@
 @example
 demo odepkg_testsuite_implrober
 @end example
+
+This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
 @end deftypefn
 
 @deftypefn {Function File} {[@var{solution}] =} odepkg_testsuite_oregonator (@var{@@solver}, @var{reltol})
@@ -400,6 +548,8 @@
 @example
 demo odepkg_testsuite_oregonator
 @end example
+
+This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
 @end deftypefn
 
 @deftypefn {Function File} {[@var{solution}] =} odepkg_testsuite_pollution (@var{@@solver}, @var{reltol})
@@ -410,6 +560,8 @@
 @example
 demo odepkg_testsuite_pollution
 @end example
+
+This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
 @end deftypefn
 
 @deftypefn {Function File} {[@var{solution}] =} odepkg_testsuite_robertson (@var{@@solver}, @var{reltol})
@@ -420,6 +572,8 @@
 @example
 demo odepkg_testsuite_robertson
 @end example
+
+This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
 @end deftypefn
 
 @deftypefn {Function File} {[@var{solution}] =} odepkg_testsuite_transistor (@var{@@solver}, @var{reltol})
@@ -430,6 +584,8 @@
 @example
 demo odepkg_testsuite_transistor
 @end example
+
+This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
 @end deftypefn
 
 @deftypefn {Function File} {[@var{ret}] =} odeplot (@var{t}, @var{y}, @var{flag})

--- a/main/odepkg/doc/odepkg.texi	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/doc/odepkg.texi	Wed Feb 08 21:02:51 2012 +0000
@@ -9,7 +9,7 @@
 @setfilename odepkg.info
 @settitle OdePkg - A package for solving ordinary differential equations and more
 @afourpaper
-@set VERSION 0.8.0
+@set VERSION 0.8.2
 @c @afourwide
 @c %*** End of the HEADER
 
@@ -107,8 +107,10 @@
 @tab Added a new solver function @command{odekdi} for the direct method (not the Krylov method) of the @file{daskr.f} solver from the authors Peter N. Brown, Alan C. Hindmarsh, Linda R. Petzold and Clement W. Ulrich that is available under a modified BSD license (without advertising clause). Ongoing work with this manual.
 @item OdePkg Version 0.5.x
 @tab Added new solver functions @command{ode23d}, @command{ode45d}, @command{ode54d} and @command{ode78d} for solving non--stiff delay differential equations (non-stiff DDEs). These solvers are based on the Runge--Kutta solvers @command{ode23}..@command{ode78}. Tests and demos have been included for this type of solvers. Added new functions @command{odeexamples}, @command{odepkg_examples_ode}, @command{odepkg_examples_dae}, @command{odepkg_examples_ide} and @command{odepkg_examples_ide}. Ongoing work with this manual.
-@item @b{(current)} Version 0.6.x
-@tab A lot of compatibility tests. Final release before version 1.0.0.
+@item OdePkg Version 0.6.x
+@tab A lot of compatibility tests, improvements, bugfixes, etc.
+@item @b{(current)} Version 0.8.x
+@tab Final releases before version 1.0.0.
 @item @b{(future)} Version 1.0.0
 @tab Completed OdePkg release 1.0.0 with M--solvers and DLD--solvers.
 @end multitable

--- a/main/odepkg/inst/odepkg_testsuite_calcmescd.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_calcmescd.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_calcmescd
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}
 

--- a/main/odepkg/inst/odepkg_testsuite_calcscd.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_calcscd.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_calcscd
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}
 

--- a/main/odepkg/inst/odepkg_testsuite_chemakzo.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_chemakzo.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_chemakzo
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}
 

--- a/main/odepkg/inst/odepkg_testsuite_hires.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_hires.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_hires
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}
 

--- a/main/odepkg/inst/odepkg_testsuite_implakzo.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_implakzo.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_implakzo
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}
 

--- a/main/odepkg/inst/odepkg_testsuite_implrober.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_implrober.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_implrober
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}
 

--- a/main/odepkg/inst/odepkg_testsuite_impltrans.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_impltrans.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_impltrans
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}
 

--- a/main/odepkg/inst/odepkg_testsuite_oregonator.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_oregonator.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_oregonator
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}
 

--- a/main/odepkg/inst/odepkg_testsuite_pollution.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_pollution.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_pollution
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}
 

--- a/main/odepkg/inst/odepkg_testsuite_robertson.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_robertson.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_robertson
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}
 

--- a/main/odepkg/inst/odepkg_testsuite_transistor.m	Wed Feb 08 08:42:25 2012 +0000
+++ b/main/odepkg/inst/odepkg_testsuite_transistor.m	Wed Feb 08 21:02:51 2012 +0000
@@ -23,9 +23,9 @@
 %# @example
 %# demo odepkg_testsuite_transistor
 %# @end example
-%# @end deftypefn
 %#
 %# This function has been ported from the "Test Set for IVP solvers" which is developed by the INdAM Bari unit project group "Codes and Test Problems for Differential Equations", coordinator F. Mazzia.
+%# @end deftypefn
 %#
 %# @seealso{odepkg}