#!/usr/bin/perl -w
# Example 1: Solve  y'' + y = 0, y(0) = 0, y'(0) = 1
#            Solution: y = sin(x)
use Math::ODE;
my $o = new Math::ODE ( file => 'data',
                        step => 0.1,
                        initial => [0,1], 
                        ODE => [ \&DE1, \&DE2 ], 
                        t0 => 0,
                        tf => 10 );
$o->evolve;
system("gnuplot -persist gnuplot.1");
# plot 'data' using 1:2, sin(x)
# y'' + y = 0
sub DE1 { my ($t,$y) = @_; return $y->[1]; }
sub DE2 { my ($t,$y) = @_; return -$y->[0]; }
#