# >>>>> new test file # ????? no tests in file # !!!!! test had an unexpected result # ----- test was skipped # ***** code for the test # Search for the unexpected results in the file # then page back to find the filename which caused it. # The result may be an unexpected failure (in which # case an error will be reported) or an unexpected # success (in which case no error will be reported). >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\inexact_solvers\inexact_newton.m ***** test warning ('off', 'OdePkg:InvalidArgument'); f = @(x) [((x(1)+1)^2 + 3*x(2)^7 + 7*(x(1)+1)); sum(x)+1]; opts.NewtonTol = 1.e-10; opts.MaxNewtonIterations = 100; sol = inexact_newton (f, [.1;.1], opts); assert (sol, [-1; 0], 1e-9) assert (f (sol), [0; 0], 1e-10) !!!!! test failed 'inexact_newton' undefined near line 6 column 9 ***** end opts.NewtonTol = 1.e-10; opts.MaxNewtonIterations = 100; opts.Jacobian = 'on'; [sol, it] = inexact_newton (@func, [.1;.1], opts); assert (sol, [-1; 0], 1e-9) assert (f (sol), [0; 0], 1e-10) !!!!! unknown test type! ***** end opts.NewtonTol = 1.e-10; opts.MaxNewtonIterations = 100; opts.Jacobian = 'on'; opts.Eta = 0.1; opts.Choice = 2; opts.Algorithm = 'bicgstab'; [sol, it] = inexact_newton (@func, [.1;.1], opts); assert (sol, [-1; 0], 1e-9) assert (func (sol), [0; 0], 1e-10) warning ('on', 'OdePkg:InvalidArgument'); !!!!! unknown test type! >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\integrate_functions\integrate_const.m ***** test warning ('off', 'OdePkg:InvalidArgument'); vsol = integrate_const (@fwe_richardson, @fpol, [0, 2], [2;0], 0.01, vopt); assert (vsol.t, [0:.01:2]', 1e-12); assert ([vsol.t(end), vsol.x(end,:)], [2, fref], 5e-2); !!!!! test failed @fwe_richardson: no function and no method found shared variables scalar structure containing the fields: vopt = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vhaveoutputfunction = 0 vhaveeventfunction = 0 vhavenonnegative = 0 vfunarguments = {}(0x0) vdirection = 1 ***** test vsol = integrate_const (@runge_kutta_23, @fpol, [0, 2], [2;0], 0.05, vopt); assert (vsol.t, [0:.05:2]', 1e-12); assert ([vsol.t(end), vsol.x(end,:)], [2, fref], 1e-2); !!!!! test failed 'integrate_const' undefined near line 3 column 10 shared variables scalar structure containing the fields: vopt = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vhaveoutputfunction = 0 vhaveeventfunction = 0 vhavenonnegative = 0 vfunarguments = {}(0x0) vdirection = 1 ***** test vopt.solver = @fsolve; vsol = integrate_const (@stormer_verlet, @hamilt, [0, 1], [1;0], 0.1, vopt); !!!!! test failed @stormer_verlet: no function and no method found shared variables scalar structure containing the fields: vopt = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vhaveoutputfunction = 0 vhaveeventfunction = 0 vhavenonnegative = 0 vfunarguments = {}(0x0) vdirection = 1 ***** test vopt.solver = @fsolve; vsol = integrate_const (@symplectic_euler, @hamilt, [0, 1], [1;0], 0.1, vopt); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed @symplectic_euler: no function and no method found shared variables scalar structure containing the fields: vopt = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vhaveoutputfunction = 0 vhaveeventfunction = 0 vhavenonnegative = 0 vfunarguments = {}(0x0) vdirection = 1 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\integrate_functions\integrate_n_steps.m ***** test warning ('off', 'OdePkg:InvalidArgument'); vsol = integrate_n_steps (@fwe_richardson, @fpol, 0, [2;0], 0.01, 2/0.01, vopt); assert (vsol.t, [0:.01:2]', 1e-12); assert ([vsol.t(end), vsol.x(end,:)], [2, fref], 5e-2); !!!!! test failed @fwe_richardson: no function and no method found shared variables scalar structure containing the fields: vopt = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vhaveoutputfunction = 0 vhaveeventfunction = 0 vhavenonnegative = 0 vfunarguments = {}(0x0) vdirection = 1 ***** test vsol = integrate_n_steps (@runge_kutta_23, @fpol, 0, [2;0], 0.05, 2/0.05, vopt); assert (vsol.t, [0:.05:2]', 1e-12); assert ([vsol.t(end), vsol.x(end,:)], [2, fref], 1e-2); !!!!! test failed 'integrate_n_steps' undefined near line 3 column 10 shared variables scalar structure containing the fields: vopt = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vhaveoutputfunction = 0 vhaveeventfunction = 0 vhavenonnegative = 0 vfunarguments = {}(0x0) vdirection = 1 ***** test vopt.solver = @fsolve; vsol = integrate_n_steps (@stormer_verlet, @hamilt, 0, [1;0], 0.1, 1/0.1, vopt); !!!!! test failed @stormer_verlet: no function and no method found shared variables scalar structure containing the fields: vopt = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vhaveoutputfunction = 0 vhaveeventfunction = 0 vhavenonnegative = 0 vfunarguments = {}(0x0) vdirection = 1 ***** test vopt.solver = @fsolve; vsol = integrate_n_steps (@symplectic_euler, @hamilt, 0, [1;0], 0.1, 2/0.1, vopt); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed @symplectic_euler: no function and no method found shared variables scalar structure containing the fields: vopt = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vhaveoutputfunction = 0 vhaveeventfunction = 0 vhavenonnegative = 0 vfunarguments = {}(0x0) vdirection = 1 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\bwe_inexact_newton.m ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** test warning ('off', 'OdePkg:InvalidArgument'); [t,y] = bwe_inexact_newton (@fpol, 0, [2;0], 0.05, opts); !!!!! test failed 'bwe_inexact_newton' undefined near line 4 column 10 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = @fjac MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) funarguments = {}(0x0) UseJacobian = no havejachandle = 1 InexactSolver = inexact_newton ***** test opts.UseJacobian = 'yes'; opts.InexactSolver = 'fsolve'; [t,y,x] = bwe_inexact_newton (@fpol, 0, [2;0], 0.05, opts); !!!!! test failed 'bwe_inexact_newton' undefined near line 5 column 12 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = @fjac MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) funarguments = {}(0x0) UseJacobian = no havejachandle = 1 InexactSolver = inexact_newton ***** test opts.UseJacobian = 'no'; opts.InexactSolver = 'inexact_newton'; opts.NewtonTol = 1.e-10; opts.MaxNewtonIterations = 1000; [t,y,x] = bwe_inexact_newton (@fpol, 0, [2;0], 0.05, opts); !!!!! test failed 'bwe_inexact_newton' undefined near line 6 column 12 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = @fjac MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) funarguments = {}(0x0) UseJacobian = no havejachandle = 1 InexactSolver = inexact_newton ***** test opts.NewtonTol = 1.e-10; opts.MaxNewtonIterations = 1000; opts.Choice = 2; [t,y,x] = bwe_inexact_newton (@fpol, 0, [2;0], 0.05, opts); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'bwe_inexact_newton' undefined near line 5 column 12 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = @fjac MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) funarguments = {}(0x0) UseJacobian = no havejachandle = 1 InexactSolver = inexact_newton >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\bwe_newton_raphson.m ***** end !!!!! unknown test type! ***** end %# Turn off output of warning messages for all tests, turn them on %# again if the last test is called !!!!! unknown test type! ***** test warning ('off', 'OdePkg:InvalidArgument'); opts.issimplified = false; opts.MaxNewtonIterations = 100; opts.NewtonTol = 1.e-6; [t,y] = bwe_newton_raphson (@fpol, 0, [2;0], 0.05, opts); !!!!! test failed 'bwe_newton_raphson' undefined near line 7 column 10 ***** test opts.issimplified = true; opts.MaxNewtonIterations = 1000; opts.NewtonTol = 1.e-10; [t,y,x] = bwe_newton_raphson (@fpol, 0, [2;0], 0.05, opts); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'bwe_newton_raphson' undefined near line 6 column 12 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\fwe_heun.m ***** end !!!!! unknown test type! ***** test [t,y] = fwe_heun (@fpol, 0, [2;0], 0.05, opts); !!!!! test failed 'fwe_heun' undefined near line 3 column 10 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) funarguments = {}(0x0) ***** test [t,y,x] = fwe_heun (@fpol, 0, [2;0], 0.05, opts); !!!!! test failed 'fwe_heun' undefined near line 3 column 12 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) funarguments = {}(0x0) >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\fwe_richardson.m ***** end !!!!! unknown test type! ***** test [t,y] = fwe_richardson (@fpol, 0, [2;0], 0.05, opts); !!!!! test failed 'fwe_richardson' undefined near line 3 column 10 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) funarguments = {}(0x0) ***** test [t,y,x] = fwe_richardson (@fpol, 0, [2;0], 0.1, opts); !!!!! test failed 'fwe_richardson' undefined near line 3 column 12 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) funarguments = {}(0x0) >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\rattle.m ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** test opts.ConstraintGradFcn = @grad_constr; [t,y] = rattle (@pendulum, 0, [5.0*cos(-pi/4);5.0*sin(-pi/4);0;0], 0.05, opts); !!!!! test failed 'rattle' undefined near line 4 column 10 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vfunarguments = {}(0x0) solver = @fsolve ConstraintFcn = @constraint ConstraintsNb = 1 ***** test opts.ConstraintGradFcn = @grad_constr; [t,y,x] = rattle (@pendulum, 0, [5.0*cos(-pi/4);5.0*sin(-pi/4);0;0], 0.1, opts); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'rattle' undefined near line 4 column 12 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vfunarguments = {}(0x0) solver = @fsolve ConstraintFcn = @constraint ConstraintsNb = 1 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\runge_kutta_23.m ***** end !!!!! unknown test type! >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\runge_kutta_45_fehlberg.m ***** end !!!!! unknown test type! ***** test [t,y] = runge_kutta_45_fehlberg (@fpol, 0, [2;0], 0.05, opts); !!!!! test failed 'runge_kutta_45_fehlberg' undefined near line 3 column 10 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vfunarguments = {}(0x0) ***** test [t,y,x] = runge_kutta_45_fehlberg (@fpol, 0, [2;0], 0.1, opts); !!!!! test failed 'runge_kutta_45_fehlberg' undefined near line 3 column 12 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vfunarguments = {}(0x0) >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\runge_kutta_78.m ***** end !!!!! unknown test type! ***** test [t,y] = runge_kutta_78 (@fpol, 0, [2;0], 0.05); !!!!! test failed 'runge_kutta_78' undefined near line 3 column 10 ***** test [t,y,x] = runge_kutta_78 (@fpol, 0, [2;0], 0.1); !!!!! test failed 'runge_kutta_78' undefined near line 3 column 12 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\spectral_var_int.m ***** end !!!!! unknown test type! ***** test vopt = odeset; vopt.vfunarguments = {}; vopt.solver = @fsolve; [t,y] = spectral_var_int (@hamilt, 0, [1;0], 0.05, vopt); !!!!! test failed 'spectral_var_int' undefined near line 6 column 10 ***** test vopt = odeset; vopt.vfunarguments = {}; vopt.solver = @fsolve; [t,y,x] = spectral_var_int (@hamilt, 0, [1;1], 0.1, vopt); !!!!! test failed 'spectral_var_int' undefined near line 6 column 12 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\stormer_verlet.m ***** end !!!!! unknown test type! ***** test warning ('off', 'OdePkg:InvalidArgument'); [t,y] = stormer_verlet (@hamilt, 0, [1;0], 0.05, opts); !!!!! test failed 'stormer_verlet' undefined near line 4 column 10 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vfunarguments = {}(0x0) solver = @fsolve Explicit = no ***** test [t,y,x] = stormer_verlet (@hamilt, 0, [1;1], 0.1, opts); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'stormer_verlet' undefined near line 3 column 12 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vfunarguments = {}(0x0) solver = @fsolve Explicit = no >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\symplectic_euler.m ***** end !!!!! unknown test type! ***** test warning ('off', 'OdePkg:InvalidArgument'); [t,y] = symplectic_euler (@hamilt, 0, [1;0], 0.05, opts); !!!!! test failed 'symplectic_euler' undefined near line 4 column 10 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vfunarguments = {}(0x0) solver = @fsolve Explicit = no ***** test [t,y,x] = symplectic_euler (@hamilt, 0, [1;1], 0.1, opts); !!!!! test failed 'symplectic_euler' undefined near line 3 column 12 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vfunarguments = {}(0x0) solver = @fsolve Explicit = no >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\velocity_verlet.m ***** end !!!!! unknown test type! ***** test warning ('off', 'OdePkg:InvalidArgument'); [t,y] = stormer_verlet (@hamilt, 0, [1;0], 0.05, opts); !!!!! test failed 'stormer_verlet' undefined near line 4 column 10 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vfunarguments = {}(0x0) solver = @fsolve Explicit = no ***** test [t,y,x] = stormer_verlet (@hamilt, 0, [1;1], 0.1, opts); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'stormer_verlet' undefined near line 3 column 12 shared variables scalar structure containing the fields: opts = scalar structure containing the fields: AbsTol = [](0x0) BDF = [](0x0) Events = [](0x0) InitialSlope = [](0x0) InitialStep = [](0x0) JConstant = [](0x0) JPattern = [](0x0) Jacobian = [](0x0) MStateDependence = [](0x0) Mass = [](0x0) MassSingular = [](0x0) MaxOrder = [](0x0) MaxStep = [](0x0) MvPattern = [](0x0) NonNegative = [](0x0) NormControl = [](0x0) OutputFcn = [](0x0) OutputSel = [](0x0) Refine = [](0x0) RelTol = [](0x0) Stats = [](0x0) Vectorized = [](0x0) vfunarguments = {}(0x0) solver = @fsolve Explicit = no >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\string_compare\fuzzy_compare.m ***** test warning ('off', 'OdePkg:InvalidArgument'); a = struct ('Green', 5, 'Yellow', 6, 'Orange', 6, 'Blue', 4, 'Red', 3); a_f = fieldnames (a); pos = fuzzy_compare ('rrd', a_f); assert (pos, 5); !!!!! test failed 'fuzzy_compare' undefined near line 6 column 9 ***** test a = ['Green'; 'Yellow'; 'Orange'; 'Blue'; 'Red']; pos = fuzzy_compare ('rrd', a); assert (pos, 5); !!!!! test failed 'fuzzy_compare' undefined near line 4 column 9 ***** test a = struct ('Green', 5, 'Yellow', 6, 'Orange', 6, 'Blue', 4, 'Red', 3); a_f = fieldnames (a); pos = fuzzy_compare ('rrd', a_f, 1); assert (pos, 5); !!!!! test failed 'fuzzy_compare' undefined near line 5 column 9 ***** test a = struct('Green', 5, 'Yellow', 6, 'Orange', 6, 'Blue', 4, 'Red', 3); a_f = fieldnames (a); pos = fuzzy_compare ('rrd', a_f, 4); assert (pos, [1 4 5]'); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'fuzzy_compare' undefined near line 5 column 9 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\string_compare\levenshtein.m ***** test ## one output argument warning ('off', 'OdePkg:InvalidArgument'); res=levenshtein('sitting','kitten'); assert (res == 3); !!!!! test failed 'levenshtein' undefined near line 4 column 7 ***** test ## two output arguments [res,matrix] = levenshtein('sitting','kitten'); ex_matrix = [0 1 2 3 4 5 6; 1 1 2 3 4 5 6; 2 2 1 2 3 4 5; 3 3 2 1 2 3 4; 4 4 3 2 1 2 3; 5 5 4 3 2 2 3; 6 6 5 4 3 3 2; 7 7 6 5 4 4 3]; assert (res == 3); assert (matrix == ex_matrix); !!!!! test failed 'levenshtein' undefined near line 3 column 17 ***** test ## two output arguments [res,matrix]=levenshtein('Sunday','Saturday'); ex_matrix = [0 1 2 3 4 5 6 7 8; 1 0 1 2 3 4 5 6 7; 2 1 1 2 2 3 4 5 6; 3 2 2 2 3 3 4 5 6; 4 3 3 3 3 4 3 4 5; 5 4 3 4 4 4 4 3 4; 6 5 4 4 5 5 5 4 3]; assert (res == 3); assert (matrix == ex_matrix); !!!!! test failed 'levenshtein' undefined near line 3 column 15 ***** test ## two output arguments case sensitiveness [res,matrix]=levenshtein('Sunday','saturday'); assert (res == 4); !!!!! test failed 'levenshtein' undefined near line 3 column 15 ***** test ## two output arguments case no-sensitiveness [res,matrix]=levenshtein('Sunday','saturday',1); assert (res > 1e6); [res,matrix]=levenshtein('Sunday','saturday',10); assert (res == 4); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'levenshtein' undefined near line 3 column 15 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\utilities\approx_Constr_Hess.m ***** end %# Turn off output of warning messages for all tests, turn them on %# again if the last test is called !!!!! unknown test type! ***** test warning ('off', 'OdePkg:InvalidArgument'); res = approx_Constr_Hess (@constraint_grad, [5.0*cos(-pi/4);5.0*sin(-pi/4)]); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'approx_Constr_Hess' undefined near line 4 column 9 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\utilities\approx_Constr_grad.m ***** end %# Turn off output of warning messages for all tests, turn them on %# again if the last test is called !!!!! unknown test type! ***** test warning ('off', 'OdePkg:InvalidArgument'); res = approx_Constr_grad (@constraint, [5.0*cos(-pi/4);5.0*sin(-pi/4)]); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'approx_Constr_grad' undefined near line 4 column 9 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\utilities\approx_Hamilt_Hess.m ***** end %# Turn off output of warning messages for all tests, turn them on %# again if the last test is called !!!!! unknown test type! ***** test warning ('off', 'OdePkg:InvalidArgument'); res = approx_Hamilt_Hess (@armonic_oscillator, 0, [1;0]); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'approx_Hamilt_Hess' undefined near line 4 column 9 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\utilities\golub_welsch.m ***** test [nodes,w,L,D] = golub_welsch(2,2); !!!!! test failed 'golub_welsch' undefined near line 3 column 18 ***** test [nodes,w,L,D] = golub_welsch(2,4); !!!!! test failed 'golub_welsch' undefined near line 3 column 18 ***** test [nodes,w,L,D] = golub_welsch(3,2); !!!!! test failed 'golub_welsch' undefined near line 3 column 18 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\ode12.m ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** test # two output arguments [t, y] = ode12 (@fpol, [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # not too many steps [t, y] = ode12 (@fpol, [0 2], [2 0]); assert (size (t) < 120); !!!!! test failed @fwe_heun: no function and no method found ***** test # anonymous function instead of real function fvdp = @(t,y) [y(2); (1 - y(1)^2) * y(2) - y(1)]; [t, y] = ode12 (fvdp, [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # string instead of function [t, y] = ode12 ("fpol", [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # extra input arguments passed through [t, y] = ode12 (@fpol, [0 2], [2 0], 12, 13, "KL"); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # empty ODEOPT structure *but* extra input arguments opt = odeset; [t, y] = ode12 (@fpol, [0 2], [2 0], opt, 12, 13, "KL"); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # Solve another anonymous function below zero vref = [0, 14.77810590694212]; [t, y] = ode12 (@(t,y) y, [-2 0], 2); assert ([t(end), y(end,:)], vref, 5e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # InitialStep option opt = odeset ("InitialStep", 1e-8); [t, y] = ode12 (@fpol, [0 0.2], [2 0], opt); assert ([t(2)-t(1)], [1e-8], 1e-9); !!!!! test failed @fwe_heun: no function and no method found ***** test # MaxStep option opt = odeset ("MaxStep", 1e-3); sol = ode12 (@fpol, [0 0.2], [2 0], opt); assert ([sol.x(5)-sol.x(4)], [1e-3], 1e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # Solve with intermediate step [t, y] = ode12 (@fpol, [0 1 2], [2 0]); assert (any((t-1) == 0)); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # Solve in backward direction starting at t=0 vref = [-1.205364552835178, 0.951542399860817]; sol = ode12 (@fpol, [0 -2], [2 0]); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # Solve in backward direction starting at t=2 vref = [-1.205364552835178, 0.951542399860817]; sol = ode12 (@fpol, [2 -2], fref); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # Solve in backward direction starting at t=2, with intermediate step vref = [-1.205364552835178, 0.951542399860817]; [t, y] = ode12 (@fpol, [2 0 -2], fref); idx = find(y < 0, 1, "first") - 1; assert ([t(idx), y(idx,:)], [0,2,0], 5e-1); assert ([t(end), y(end,:)], [-2, vref], 5e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # Solve another anonymous function in backward direction vref = [-1, 0.367879437558975]; sol = ode12 (@(t,y) y, [0 -1], 1); assert ([sol.x(end); sol.y(:,end)], vref', 1e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # Solve another anonymous function below zero vref = [0, 14.77810590694212]; sol = ode12 (@(t,y) y, [-2 0], 2); assert ([sol.x(end); sol.y(:,end)], vref', 5e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # Solve in backward direction starting at t=0 with MaxStep option vref = [-1.205364552835178, 0.951542399860817]; opt = odeset ("MaxStep", 1e-3); sol = ode12 (@fpol, [0 -2], [2 0], opt); assert ([abs(sol.x(8)-sol.x(7))], [1e-3], 1e-2); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # AbsTol option opt = odeset ("AbsTol", 1e-5); sol = ode12 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # AbsTol and RelTol option opt = odeset ("AbsTol", 1e-8, "RelTol", 1e-8); sol = ode12 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-2); !!!!! test failed @fwe_heun: no function and no method found ***** test # RelTol and NormControl option -- higher accuracy opt = odeset ("RelTol", 1e-8, "NormControl", "on"); sol = ode12 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 2e-3); !!!!! test failed @fwe_heun: no function and no method found ***** test # Keeps initial values while integrating opt = odeset ("NonNegative", 2); sol = ode12 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; 2; 0], 0.5); !!!!! test failed @fwe_heun: no function and no method found ***** test # Details of OutputSel and Refine can't be tested opt = odeset ("OutputFcn", @fout, "OutputSel", 1, "Refine", 5); sol = ode12 (@fpol, [0 2], [2 0], opt); !!!!! test failed @fwe_heun: no function and no method found ***** test # Stats must add further elements in sol opt = odeset ("Stats", "on"); sol = ode12 (@fpol, [0 2], [2 0], opt); assert (isfield (sol, "stats")); assert (isfield (sol.stats, "nsteps")); !!!!! test failed @fwe_heun: no function and no method found ***** test # Events option add further elements in sol opt = odeset ("Events", @feve); sol = ode12 (@fpol, [0 10], [2 0], opt); assert (isfield (sol, "ie")); assert (sol.ie(1), 2); assert (isfield (sol, "xe")); assert (isfield (sol, "ye")); !!!!! test failed @fwe_heun: no function and no method found ***** test # Events option, now stop integration warning ("off", "integrate_adaptive:unexpected_termination", "local"); opt = odeset ("Events", @fevn, "NormControl", "on"); sol = ode12 (@fpol, [0 10], [2 0], opt); assert ([sol.ie, sol.xe, sol.ye], [2.0, 2.496110, -0.830550, -2.677589], 6e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # Events option, five output arguments warning ("off", "integrate_adaptive:unexpected_termination", "local"); opt = odeset ("Events", @fevn, "NormControl", "on"); [t, y, vxe, ye, vie] = ode12 (@fpol, [0 10], [2 0], opt); assert ([vie, vxe, ye], [2.0, 2.496110, -0.830550, -2.677589], 6e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # Mass option as function opt = odeset ("Mass", @fmas); sol = ode12 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # Mass option as matrix opt = odeset ("Mass", eye (2,2)); sol = ode12 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # Mass option as sparse matrix opt = odeset ("Mass", sparse (eye (2,2))); sol = ode12 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # Mass option as function and sparse matrix opt = odeset ("Mass", @fmsa); sol = ode12 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # Mass option as function and MStateDependence opt = odeset ("Mass", @fmas, "MStateDependence", "strong"); sol = ode12 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_heun: no function and no method found ***** test # Check that imaginary part of solution does not get inverted sol = ode12 (@(x,y) 1, [0 1], 1i); assert (imag (sol.y), ones (size (sol.y))) [x, y] = ode12 (@(x,y) 1, [0 1], 1i); assert (imag (y), ones (size (y))) !!!!! test failed @fwe_heun: no function and no method found >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\ode23d.m >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\ode23s.m ***** test test1=@(t,y) t - y + 1; [vt, vy] = ode23s (test1, [0 10], [1]); assert ([vt(end), vy(end)], [10, exp(-10) + 10], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** test # two output arguments [t, y] = ode23s (@fpol, [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # anonymous function instead of real function fvdp = @(t,y) [y(2); 10 * (1 - y(1)^2) * y(2) - y(1)]; [t, y] = ode23s (fvdp, [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # extra input arguments passed through [t, y] = ode23s (@fpol, [0 2], [2 0], 12, 13, "KL"); assert ([t(end), y(end,:)], [2, fref], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # empty OdePkg structure *but* extra input arguments opt = odeset; [t, y] = ode23s (@fpol, [0 2], [2 0], opt, 12, 13, "KL"); assert ([t(end), y(end,:)], [2, fref], 1e-2); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # InitialStep option opt = odeset ("InitialStep", 1e-8); [t, y] = ode23s (@fpol, [0 0.2], [2 0], opt); assert ([t(2)-t(1)], [1e-8], 1e-9); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # MaxStep option opt = odeset ("MaxStep", 1e-3); sol = ode23s (@fpol, [0 0.2], [2 0], opt); assert ([sol.x(5)-sol.x(4)], [1e-3], 1e-4); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # AbsTol option opt = odeset ("AbsTol", 1e-5); sol = ode23s (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # AbsTol and RelTol option opt = odeset ("AbsTol", 1e-8, "RelTol", 1e-8); sol = ode23s (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # RelTol and NormControl option -- higher accuracy opt = odeset ("RelTol", 1e-8, "NormControl", "on"); sol = ode23s (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-4); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # Details of OutputSel and Refine can't be tested opt = odeset ("OutputFcn", @fout, "OutputSel", 1, "Refine", 5); sol = ode23s (@fpol, [0 2], [2 0], opt); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # Stats must add further elements in sol opt = odeset ("Stats", "on"); sol = ode23s (@fpol, [0 2], [2 0], opt); assert (isfield (sol, "stats")); assert (isfield (sol.stats, "nsteps")); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # Events option add further elements in sol opt = odeset ("Events", @feve); sol = ode23s (@fpol, [0 10], [2 0], opt); assert (isfield (sol, "ie")); assert (sol.ie(1), 2); assert (isfield (sol, "xe")); assert (isfield (sol, "ye")); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # Mass option as function opt = odeset ("Mass", @fmas); sol = ode23s (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # Mass option as matrix opt = odeset ("Mass", eye (2,2)); sol = ode23s (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # Mass option as sparse matrix opt = odeset ("Mass", sparse (eye (2,2))); sol = ode23s (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # Mass option as function and sparse matrix opt = odeset ("Mass", @fmsa); sol = ode23s (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # Jacobian option as function opt = odeset ('Jacobian', @jac); sol = ode23s (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # Sparse Jacobian jac = @(t, y) sparse ([0 1; -20*y(1)*y(2)-1, 10*(1-y(1)^2)]); opt = odeset ('Jacobian', jac); sol = ode23s (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_23s: no function and no method found ***** test # Check that imaginary part of solution does not get inverted sol = ode23s (@(x,y) 1, [0 1], 1i); assert (imag (sol.y), ones (size (sol.y))) [x, y] = ode23s (@(x,y) 1, [0 1], 1i); assert (imag (y), ones (size (y))) !!!!! test failed @runge_kutta_23s: no function and no method found >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\ode45d.m >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\ode78.m ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** test # two output arguments [t, y] = ode78 (@fpol, [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 1e-2); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # not too many steps [t, y] = ode78 (@fpol, [0 2], [2 0]); assert (size (t) < 15); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # anonymous function instead of real function fvdp = @(t,y) [y(2); (1 - y(1)^2) * y(2) - y(1)]; [t, y] = ode78 (fvdp, [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 1e-2); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # string instead of function [t, y] = ode78 ("fpol", [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 1e-2); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # extra input arguments passed through [t, y] = ode78 (@fpol, [0 2], [2 0], 12, 13, "KL"); assert ([t(end), y(end,:)], [2, fref], 1e-2); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # empty ODEOPT structure *but* extra input arguments opt = odeset; [t, y] = ode78 (@fpol, [0 2], [2 0], opt, 12, 13, "KL"); assert ([t(end), y(end,:)], [2, fref], 1e-2); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Solve another anonymous function below zero vref = [0, 14.77810590694212]; [t, y] = ode78 (@(t,y) y, [-2 0], 2); assert ([t(end), y(end,:)], vref, 1e-1); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # InitialStep option opt = odeset ("InitialStep", 1e-8); [t, y] = ode78 (@fpol, [0 0.2], [2 0], opt); assert ([t(2)-t(1)], [1e-8], 1e-9); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # MaxStep option opt = odeset ("MaxStep", 1e-3); sol = ode78 (@fpol, [0 0.2], [2 0], opt); assert ([sol.x(5)-sol.x(4)], [1e-3], 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Solve with intermediate step [t, y] = ode78 (@fpol, [0 1 2], [2 0]); assert (any((t-1) == 0)); assert ([t(end), y(end,:)], [2, fref], 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Solve in backward direction starting at t=0 vref = [-1.205364552835178, 0.951542399860817]; sol = ode78 (@fpol, [0 -2], [2 0]); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-2); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Solve in backward direction starting at t=2 vref = [-1.205364552835178, 0.951542399860817]; sol = ode78 (@fpol, [2 -2], fref); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-2); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Solve in backward direction starting at t=2, with intermediate step vref = [-1.205364552835178, 0.951542399860817]; [t, y] = ode78 (@fpol, [2 0 -2], fref); idx = find(y < 0, 1, "first") - 1; assert ([t(idx), y(idx,:)], [0,2,0], 1e-2); assert ([t(end), y(end,:)], [-2, vref], 1e-2); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Solve another anonymous function in backward direction vref = [-1, 0.367879437558975]; sol = ode78 (@(t,y) y, [0 -1], 1); assert ([sol.x(end); sol.y(:,end)], vref', 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Solve another anonymous function below zero vref = [0, 14.77810590694212]; sol = ode78 (@(t,y) y, [-2 0], 2); assert ([sol.x(end); sol.y(:,end)], vref', 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Solve in backward direction starting at t=0 with MaxStep option vref = [-1.205364552835178, 0.951542399860817]; opt = odeset ("MaxStep", 1e-3); sol = ode78 (@fpol, [0 -2], [2 0], opt); assert ([abs(sol.x(8)-sol.x(7))], [1e-3], 1e-3); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # AbsTol option opt = odeset ("AbsTol", 1e-5); sol = ode78 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # AbsTol and RelTol option opt = odeset ("AbsTol", 1e-8, "RelTol", 1e-8); sol = ode78 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # RelTol and NormControl option -- higher accuracy opt = odeset ("RelTol", 1e-8, "NormControl", "on"); sol = ode78 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-5); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Keeps initial values while integrating opt = odeset ("NonNegative", 2); sol = ode78 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; 2; 0], 0.5); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Details of OutputSel and Refine can't be tested opt = odeset ("OutputFcn", @fout, "OutputSel", 1, "Refine", 5); sol = ode78 (@fpol, [0 2], [2 0], opt); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Stats must add further elements in sol opt = odeset ("Stats", "on"); sol = ode78 (@fpol, [0 2], [2 0], opt); assert (isfield (sol, "stats")); assert (isfield (sol.stats, "nsteps")); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Events option add further elements in sol opt = odeset ("Events", @feve); sol = ode78 (@fpol, [0 10], [2 0], opt); assert (isfield (sol, "ie")); assert (sol.ie(1), 2); assert (isfield (sol, "xe")); assert (isfield (sol, "ye")); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Events option, now stop integration warning ("off", "integrate_adaptive:unexpected_termination", "local"); opt = odeset ("Events", @fevn, "NormControl", "on"); sol = ode78 (@fpol, [0 10], [2 0], opt); assert ([sol.ie, sol.xe, sol.ye], [2.0, 2.496110, -0.830550, -2.677589], 6e-1); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Events option, five output arguments warning ("off", "integrate_adaptive:unexpected_termination", "local"); opt = odeset ("Events", @fevn, "NormControl", "on"); [t, y, vxe, ye, vie] = ode78 (@fpol, [0 10], [2 0], opt); assert ([vie, vxe, ye], [2.0, 2.496110, -0.830550, -2.677589], 6e-1); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Mass option as function opt = odeset ("Mass", @fmas); sol = ode78 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Mass option as matrix opt = odeset ("Mass", eye (2,2)); sol = ode78 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Mass option as sparse matrix opt = odeset ("Mass", sparse (eye (2,2))); sol = ode78 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Mass option as function and sparse matrix opt = odeset ("Mass", @fmsa); sol = ode78 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Mass option as function and MStateDependence opt = odeset ("Mass", @fmas, "MStateDependence", "strong"); sol = ode78 (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-3); !!!!! test failed @runge_kutta_78: no function and no method found ***** test # Check that imaginary part of solution does not get inverted sol = ode78 (@(x,y) 1, [0 1], 1i); assert (imag (sol.y), ones (size (sol.y))) [x, y] = ode78 (@(x,y) 1, [0 1], 1i); assert (imag (y), ones (size (y))) !!!!! test failed @runge_kutta_78: no function and no method found >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\ode78d.m >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odeRATTLE.m ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** test %# two output arguments [vt, vy] = odeRATTLE (@pendulum, @constraint, [0 0.4], [L*cos(-pi/4) L*sin(-pi/4) 0 0]); assert ([vt(end), sqrt(vy(end,1)^2+vy(end,2)^2)], [0.4, L], 1e-3); !!!!! test failed @rattle: no function and no method found shared variables scalar structure containing the fields: m = 2 L = 5 g = 9.8060 ***** test %# extra input arguments passed through [vt, vy] = odeRATTLE (@pendulum, @constraint, [0 0.4], [L*cos(-pi/4) L*sin(-pi/4) 0 0], 12, 13, 'KL'); assert ([vt(end), sqrt(vy(end,1)^2+vy(end,2)^2)], [0.4, L], 1e-2); !!!!! test failed @rattle: no function and no method found shared variables scalar structure containing the fields: m = 2 L = 5 g = 9.8060 ***** test %# empty OdePkg structure *but* extra input arguments opt = odeset(); [vt, vy] = odeRATTLE (@pendulum, @constraint, [0 0.4], [L*cos(-pi/4) L*sin(-pi/4) 0 0], opt, 12, 13, 'KL'); assert ([vt(end), sqrt(vy(end,1)^2+vy(end,2)^2)], [0.4, L], 1e-2); !!!!! test failed @rattle: no function and no method found shared variables scalar structure containing the fields: m = 2 L = 5 g = 9.8060 ***** test # Solve with intermidiate step [vt, vy] = odeRATTLE (@pendulum, @constraint, [0 0.2 0.4], [L*cos(-pi/4) L*sin(-pi/4) 0 0]); assert ([vt(end), sqrt(vy(end,1)^2+vy(end,2)^2)], [0.4, L], 1e-2); !!!!! test failed @rattle: no function and no method found shared variables scalar structure containing the fields: m = 2 L = 5 g = 9.8060 ***** test %# use_jacobian option opt = odeset(); opt.UseJacobian = 'yes'; [vt, vy] = odeRATTLE (@pendulum_with_H, @constraint, [0 0.4], [L*cos(-pi/4) L*sin(-pi/4) 0 0], opt); assert ([vt(end), sqrt(vy(end,1)^2+vy(end,2)^2)], [0.4, L], 1e-2); !!!!! test failed @rattle: no function and no method found shared variables scalar structure containing the fields: m = 2 L = 5 g = 9.8060 ***** test %# use approximated jacobian opt = odeset(); opt.UseJacobian = 'yes'; [vt, vy] = odeRATTLE (@pendulum, @constraint, [0 0.4], [L*cos(-pi/4) L*sin(-pi/4) 0 0], opt); assert ([vt(end), sqrt(vy(end,1)^2+vy(end,2)^2)], [0.4, L], 1e-2); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed @rattle: no function and no method found shared variables scalar structure containing the fields: m = 2 L = 5 g = 9.8060 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odeSE.m ***** end !!!!! unknown test type! ***** end %# Turn off output of warning messages for all tests, turn them on %# again if the last test is called !!!!! unknown test type! ***** test %# two output arguments [vt, vy] = odeSE (@armonic, [0 1], [1 0]); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 5e-2); !!!!! test failed @symplectic_euler: no function and no method found ***** test %# extra input arguments passed through [vt, vy] = odeSE (@armonic, [0 1], [1 0], 12, 13, 'KL'); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 5e-2); !!!!! test failed @symplectic_euler: no function and no method found ***** test %# empty OdePkg structure *but* extra input arguments vopt = odeset; [vt, vy] = odeSE (@armonic, [0 1], [1 0], vopt, 12, 13, 'KL'); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 5e-2); !!!!! test failed @symplectic_euler: no function and no method found ***** test %# explicit option vopt = odeset('Explicit','yes'); [vt, vy] = odeSE (@armonic, [0 1], [1 0], vopt); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 5e-2); !!!!! test failed @symplectic_euler: no function and no method found ***** test %# use_jacobian option vopt = odeset('Explicit','no','UseJacobian','yes'); [vt, vy] = odeSE (@armonic_with_H, [0 1], [1 0], vopt); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 5e-2); !!!!! test failed @symplectic_euler: no function and no method found ***** test %# two output arguments [vt, vy] = odeSE (@armonic, [0 1], [1 0]); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 5e-2); !!!!! test failed @symplectic_euler: no function and no method found ***** test # Solve with intermidiate step [vt, vy] = odeSE (@armonic, [0 .5 1], [1 0]); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 5e-2); !!!!! test failed @symplectic_euler: no function and no method found ***** test %# InitialStep option vopt = odeset ('InitialStep', 1e-8); [vt, vy] = odeSE (@armonic, [0 0.2], [1 0], vopt); assert ([vt(2)-vt(1)], [1e-8], 1e-9); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed @symplectic_euler: no function and no method found >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odeSPVI.m ***** end !!!!! unknown test type! ***** end %# Turn off output of warning messages for all tests, turn them on %# again if the last test is called !!!!! unknown test type! ***** test %# two output arguments [vt, vy] = odeSPVI (@armonic, [0 1], [1 0] ,odeset('PolynomialDegree',2,'QuadratureOrder',3)); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed 'golub_welsch' undefined near line 325 column 41 ***** test %# extra input arguments passed through [vt, vy] = odeSPVI (@armonic, [0 1], [1 0], 12, 13, 'KL'); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed 'golub_welsch' undefined near line 325 column 41 ***** test %# empty OdePkg structure *but* extra input arguments vopt = odeset(); [vt, vy] = odeSPVI (@armonic, [0 1], [1 0], vopt, 12, 13, 'KL'); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed 'golub_welsch' undefined near line 325 column 41 ***** test %# use_jacobian option vopt = odeset('PolynomialDegree',2,'QuadratureOrder',3,'UseJacobian','yes'); [vt, vy] = odeSPVI (@armonic_with_H, [0 1], [1 0], vopt); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed 'golub_welsch' undefined near line 325 column 41 ***** test %# use approximated jacobian vopt = odeset('PolynomialDegree',2,'QuadratureOrder',3,'UseJacobian','yes'); [vt, vy] = odeSPVI (@armonic, [0 1], [1 0], vopt); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed 'golub_welsch' undefined near line 325 column 41 ***** test %# InitialStep option vopt = odeset ('InitialStep', 1e-8); [vt, vy] = odeSPVI (@armonic, [0 0.2], [1 0], vopt); assert ([vt(2)-vt(1)], [1e-8], 1e-9); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed 'golub_welsch' undefined near line 325 column 41 >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odeSV.m ***** end !!!!! unknown test type! ***** end %# Turn off output of warning messages for all tests, turn them on %# again if the last test is called !!!!! unknown test type! ***** test %# two output arguments [vt, vy] = odeSV (@armonic, [0 1], [1 0]); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed @stormer_verlet: no function and no method found ***** test %# extra input arguments passed through [vt, vy] = odeSV (@armonic, [0 1], [1 0], 12, 13, 'KL'); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed @stormer_verlet: no function and no method found ***** test %# empty OdePkg structure *but* extra input arguments vopt = odeset; [vt, vy] = odeSV (@armonic, [0 1], [1 0], vopt, 12, 13, 'KL'); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed @stormer_verlet: no function and no method found ***** test %# explicit option vopt = odeset('Explicit','yes'); [vt, vy] = odeSV (@armonic, [0 1], [1 0], vopt); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed @stormer_verlet: no function and no method found ***** test %# use_jacobian option vopt = odeset('Explicit','no','UseJacobian','yes'); [vt, vy] = odeSV (@armonic_with_H, [0 1], [1 0], vopt); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed @stormer_verlet: no function and no method found ***** test %# two output arguments [vt, vy] = odeSV (@armonic, [0 1], [1 0]); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed @stormer_verlet: no function and no method found ***** test %# InitialStep option vopt = odeset ('InitialStep', 1e-8); [vt, vy] = odeSV (@armonic, [0 0.2], [1 0], vopt); assert ([vt(2)-vt(1)], [1e-8], 1e-9); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed @stormer_verlet: no function and no method found >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odeVV.m ***** end %# Turn off output of warning messages for all tests, turn them on %# again if the last test is called !!!!! unknown test type! ***** test %# two output arguments [vt, vy] = odeVV (@armonic, [0 1], [1 0]); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed @velocity_verlet: no function and no method found ***** test %# extra input arguments passed through [vt, vy] = odeVV (@armonic, [0 1], [1 0], 12, 13, 'KL'); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed @velocity_verlet: no function and no method found ***** test %# empty OdePkg structure *but* extra input arguments vopt = odeset; [vt, vy] = odeVV (@armonic, [0 1], [1 0], vopt, 12, 13, 'KL'); assert ([vt(end), vy(end,1)], [1, cos(vt(end))], 1e-2); !!!!! test failed @velocity_verlet: no function and no method found ***** test %# InitialStep option vopt = odeset ('InitialStep', 1e-8); [vt, vy] = odeVV (@armonic, [0 0.2], [1 0], vopt); assert ([vt(2)-vt(1)], [1e-8], 1e-9); warning ('on', 'OdePkg:InvalidArgument'); !!!!! test failed @velocity_verlet: no function and no method found >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odebwe.m ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** test # two output arguments [t, y] = odebwe (@fpol, [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # not too many steps [t, y] = odebwe (@fpol, [0 2], [2 0]); assert (size (t) < 120); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # anonymous function instead of real function fvdp = @(t,y) [y(2); (1 - y(1)^2) * y(2) - y(1)]; [t, y] = odebwe (fvdp, [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # string instead of function [t, y] = odebwe ("fpol", [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # extra input arguments passed through [t, y] = odebwe (@fpol, [0 2], [2 0], 12, 13, "KL"); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # empty ODEOPT structure *but* extra input arguments opt = odeset; [t, y] = odebwe (@fpol, [0 2], [2 0], opt, 12, 13, "KL"); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Solve another anonymous function below zero vref = [0, 14.77810590694212]; [t, y] = odebwe (@(t,y) y, [-2 0], 2); assert ([t(end), y(end,:)], vref, 5e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # InitialStep option opt = odeset ("InitialStep", 1e-8); [t, y] = odebwe (@fpol, [0 0.2], [2 0], opt); assert ([t(2)-t(1)], [1e-8], 1e-9); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # MaxStep option opt = odeset ("MaxStep", 1e-3); sol = odebwe (@fpol, [0 0.2], [2 0], opt); assert ([sol.x(5)-sol.x(4)], [1e-3], 1e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Solve with intermediate step [t, y] = odebwe (@fpol, [0 1 2], [2 0]); assert (any((t-1) == 0)); assert ([t(end), y(end,:)], [2, fref], 5e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Solve in backward direction starting at t=0 vref = [-1.205364552835178, 0.951542399860817]; sol = odebwe (@fpol, [0 -2], [2 0]); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Solve in backward direction starting at t=2 vref = [-1.205364552835178, 0.951542399860817]; sol = odebwe (@fpol, [2 -2], fref); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Solve in backward direction starting at t=2, with intermediate step vref = [-1.205364552835178, 0.951542399860817]; [t, y] = odebwe (@fpol, [2 0 -2], fref); idx = find(y < 0, 1, "first") - 1; assert ([t(idx), y(idx,:)], [0,2,0], 5e-1); assert ([t(end), y(end,:)], [-2, vref], 5e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Solve another anonymous function in backward direction vref = [-1, 0.367879437558975]; sol = odebwe (@(t,y) y, [0 -1], 1); assert ([sol.x(end); sol.y(:,end)], vref', 1e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Solve another anonymous function below zero vref = [0, 14.77810590694212]; sol = odebwe (@(t,y) y, [-2 0], 2); assert ([sol.x(end); sol.y(:,end)], vref', 5e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Solve in backward direction starting at t=0 with MaxStep option vref = [-1.205364552835178, 0.951542399860817]; opt = odeset ("MaxStep", 1e-3); sol = odebwe (@fpol, [0 -2], [2 0], opt); assert ([abs(sol.x(8)-sol.x(7))], [1e-3], 1e-2); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # AbsTol option opt = odeset ("AbsTol", 1e-5); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # AbsTol and RelTol option opt = odeset ("AbsTol", 1e-8, "RelTol", 1e-8); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # RelTol and NormControl option -- higher accuracy opt = odeset ("RelTol", 1e-8, "NormControl", "on"); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 2e-3); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Details of OutputSel and Refine can't be tested opt = odeset ("OutputFcn", @fout, "OutputSel", 1, "Refine", 5); sol = odebwe (@fpol, [0 2], [2 0], opt); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Stats must add further elements in sol opt = odeset ("Stats", "on"); sol = odebwe (@fpol, [0 2], [2 0], opt); assert (isfield (sol, "stats")); assert (isfield (sol.stats, "nsteps")); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Events option add further elements in sol opt = odeset ("Events", @feve); sol = odebwe (@fpol, [0 10], [2 0], opt); assert (isfield (sol, "ie")); assert (sol.ie(1), 2); assert (isfield (sol, "xe")); assert (isfield (sol, "ye")); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Events option, now stop integration warning ("off", "integrate_adaptive:unexpected_termination", "local"); opt = odeset ("Events", @fevn, "NormControl", "on"); sol = odebwe (@fpol, [0 10], [2 0], opt); assert ([sol.ie, sol.xe, sol.ye], [2.0, 2.496110, -0.830550, -2.677589], 6e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Events option, five output arguments warning ("off", "integrate_adaptive:unexpected_termination", "local"); opt = odeset ("Events", @fevn, "NormControl", "on"); [t, y, vxe, ye, vie] = odebwe (@fpol, [0 10], [2 0], opt); assert ([vie, vxe, ye], [2.0, 2.496110, -0.830550, -2.677589], 6e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Mass option as function opt = odeset ("Mass", @fmas); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Mass option as matrix opt = odeset ("Mass", eye (2,2)); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Mass option as sparse matrix opt = odeset ("Mass", sparse (eye (2,2))); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Mass option as function and sparse matrix opt = odeset ("Mass", @fmsa); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test # Mass option as function and MStateDependence opt = odeset ("Mass", @fmas, "MStateDependence", "strong"); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test %# Jacobian option opt = odeset ('Jacobian', @fjac); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 5e-2); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test %# Jacobian option and sparse return value opt = odeset ('Jacobian', @fjcc); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @bwe_newton_raphson: no function and no method found ***** test %# inexact_newton opt = odeset ('InexactSolver', 'inexact_newton', 'NewtonTol', 1e-8, 'Algorithm', 'gmres'); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 5e-2); !!!!! test failed @bwe_inexact_newton: no function and no method found ***** test %# inexact_newton and pcg opt = odeset ('InexactSolver', 'inexact_newton', 'NewtonTol', 1e-8, 'Algorithm', 'gmres'); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 5e-2); !!!!! test failed @bwe_inexact_newton: no function and no method found ***** test %# inexact_newton and UseJacobian opt = odeset ('InexactSolver', 'inexact_newton', 'Jacobian', @fjac, 'UseJacobian', 'yes'); sol = odebwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 5e-2); %# test for JPattern option is missing %# test for Vectorized option is missing !!!!! test failed @bwe_inexact_newton: no function and no method found ***** test # Check that imaginary part of solution does not get inverted sol = odebwe (@(x,y) 1, [0 1], 1i); assert (imag (sol.y), ones (size (sol.y))) [x, y] = odebwe (@(x,y) 1, [0 1], 1i); assert (imag (y), ones (size (y))) !!!!! test failed @bwe_newton_raphson: no function and no method found >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odefwe.m ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** end !!!!! unknown test type! ***** test # two output arguments [t, y] = odefwe (@fpol, [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # not too many steps [t, y] = odefwe (@fpol, [0 2], [2 0]); assert (size (t) < 120); !!!!! test failed @fwe_richardson: no function and no method found ***** test # anonymous function instead of real function fvdp = @(t,y) [y(2); (1 - y(1)^2) * y(2) - y(1)]; [t, y] = odefwe (fvdp, [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # string instead of function [t, y] = odefwe ("fpol", [0 2], [2 0]); assert ([t(end), y(end,:)], [2, fref], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # extra input arguments passed through [t, y] = odefwe (@fpol, [0 2], [2 0], 12, 13, "KL"); assert ([t(end), y(end,:)], [2, fref], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # empty ODEOPT structure *but* extra input arguments opt = odeset; [t, y] = odefwe (@fpol, [0 2], [2 0], opt, 12, 13, "KL"); assert ([t(end), y(end,:)], [2, fref], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Solve another anonymous function below zero vref = [0, 14.77810590694212]; [t, y] = odefwe (@(t,y) y, [-2 0], 2); assert ([t(end), y(end,:)], vref, 6e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # InitialStep option opt = odeset ("InitialStep", 1e-8); [t, y] = odefwe (@fpol, [0 0.2], [2 0], opt); assert ([t(2)-t(1)], [1e-8], 1e-9); !!!!! test failed @fwe_richardson: no function and no method found ***** test # MaxStep option opt = odeset ("MaxStep", 1e-3); sol = odefwe (@fpol, [0 0.2], [2 0], opt); assert ([sol.x(5)-sol.x(4)], [1e-3], 1e-2); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Solve with intermediate step [t, y] = odefwe (@fpol, [0 1 2], [2 0]); assert (any((t-1) == 0)); assert ([t(end), y(end,:)], [2, fref], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Solve in backward direction starting at t=0 vref = [-1.205364552835178, 0.951542399860817]; sol = odefwe (@fpol, [0 -2], [2 0]); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Solve in backward direction starting at t=2 vref = [-1.205364552835178, 0.951542399860817]; sol = odefwe (@fpol, [2 -2], fref); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Solve in backward direction starting at t=2, with intermediate step vref = [-1.205364552835178, 0.951542399860817]; [t, y] = odefwe (@fpol, [2 0 -2], fref); idx = find(y < 0, 1, "first") - 1; assert ([t(idx), y(idx,:)], [0,2,0], 5e-1); assert ([t(end), y(end,:)], [-2, vref], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Solve another anonymous function in backward direction vref = [-1, 0.367879437558975]; sol = odefwe (@(t,y) y, [0 -1], 1); assert ([sol.x(end); sol.y(:,end)], vref', 1e-2); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Solve in backward direction starting at t=0 with MaxStep option vref = [-1.205364552835178, 0.951542399860817]; opt = odeset ("MaxStep", 1e-3); sol = odefwe (@fpol, [0 -2], [2 0], opt); assert ([abs(sol.x(8)-sol.x(7))], [1e-3], 1e-2); assert ([sol.x(end); sol.y(:,end)], [-2; vref'], 1e-2); !!!!! test failed @fwe_richardson: no function and no method found ***** test # AbsTol option opt = odeset ("AbsTol", 1e-5); sol = odefwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # AbsTol and RelTol option opt = odeset ("AbsTol", 1e-8, "RelTol", 1e-8); sol = odefwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-2); !!!!! test failed @fwe_richardson: no function and no method found ***** test # RelTol and NormControl option -- higher accuracy opt = odeset ("RelTol", 1e-8, "NormControl", "on"); sol = odefwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 2e-3); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Keeps initial values while integrating opt = odeset ("NonNegative", 2); sol = odefwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; 2; 0], 0.5); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Details of OutputSel and Refine can't be tested opt = odeset ("OutputFcn", @fout, "OutputSel", 1, "Refine", 5); sol = odefwe (@fpol, [0 2], [2 0], opt); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Stats must add further elements in sol opt = odeset ("Stats", "on"); sol = odefwe (@fpol, [0 2], [2 0], opt); assert (isfield (sol, "stats")); assert (isfield (sol.stats, "nsteps")); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Events option add further elements in sol opt = odeset ("Events", @feve); sol = odefwe (@fpol, [0 10], [2 0], opt); assert (isfield (sol, "ie")); assert (sol.ie(1), 2); assert (isfield (sol, "xe")); assert (isfield (sol, "ye")); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Events option, now stop integration warning ("off", "integrate_adaptive:unexpected_termination", "local"); opt = odeset ("Events", @fevn, "NormControl", "on"); sol = odefwe (@fpol, [0 10], [2 0], opt); assert ([sol.ie, sol.xe, sol.ye], [2.0, 2.496110, -0.830550, -2.677589], 6e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Events option, five output arguments warning ("off", "integrate_adaptive:unexpected_termination", "local"); opt = odeset ("Events", @fevn, "NormControl", "on"); [t, y, vxe, ye, vie] = odefwe (@fpol, [0 10], [2 0], opt); assert ([vie, vxe, ye], [2.0, 2.496110, -0.830550, -2.677589], 6e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Mass option as function opt = odeset ("Mass", @fmas); sol = odefwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Mass option as matrix opt = odeset ("Mass", eye (2,2)); sol = odefwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Mass option as sparse matrix opt = odeset ("Mass", sparse (eye (2,2))); sol = odefwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Mass option as function and sparse matrix opt = odeset ("Mass", @fmsa); sol = odefwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Mass option as function and MStateDependence opt = odeset ("Mass", @fmas, "MStateDependence", "strong"); sol = odefwe (@fpol, [0 2], [2 0], opt); assert ([sol.x(end); sol.y(:,end)], [2; fref'], 1e-1); !!!!! test failed @fwe_richardson: no function and no method found ***** test # Check that imaginary part of solution does not get inverted sol = odefwe (@(x,y) 1, [0 1], 1i); assert (imag (sol.y), ones (size (sol.y))) [x, y] = odefwe (@(x,y) 1, [0 1], 1i); assert (imag (y), ones (size (y))) !!!!! test failed @fwe_richardson: no function and no method found >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_calcmescd.m ***** assert (odepkg_testsuite_calcmescd ([1, 2, 3], [1, 2, 3], NaN, NaN), Inf); !!!!! test failed ASSERT errors for: assert (odepkg_testsuite_calcmescd ([1, 2, 3], [1, 2, 3], NaN, NaN),Inf) Location | Observed | Expected | Reason () NaN Inf 'NaN' mismatch >>>>> processing C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_calcscd.m ***** assert (odepkg_testsuite_calcscd ([1, 2, 3], [1+1e-6, 2, 3], 1, 1), 6.5, 0.2); !!!!! test failed ASSERT errors for: assert (odepkg_testsuite_calcscd ([1, 2, 3], [1 + 1e-6, 2, 3], 1, 1),6.5,0.2) Location | Observed | Expected | Reason () 6 6.5 Abs err 0.5 exceeds tol 0.2 Files with no tests: C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\steppers\runge_kutta_23s.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\utilities\__dfxpdp__.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\bvp4c.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odeexamples.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odephas2.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odephas3.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_examples_dae.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_examples_dde.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_examples_ide.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_examples_ode.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_chemakzo.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_hires.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_implakzo.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_implrober.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_impltrans.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_oregonator.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_pollution.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_robertson.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odepkg_testsuite_transistor.m C:\Octave\OCTAVE~1.1\share\octave\packages\odepkg-0.9.1\odeprint.m