use Test::More tests => 1; use File::Spec; use lib File::Spec->catfile("..","lib"); use Math::ODE; use Data::Dumper; # analytic solution is y(x) = 2 e^{-x} my $o = new Math::ODE ( file => 'data', step => 0.1, initial => [2], ODE => [ \&DE1 ], t0 => 0, tf => 1 ); $o->evolve; my $res = abs($o->error - $o->{step}**4 ); ok( $res < 1e-12, "returns correct amount of error , res=$res" ); sub DE1 { my ($t,$y) = @_; return -$y->[0]; }