#!/usr/athena/bin/perl

use strict;

my $eqn = 'r * x * (1 - x)';
my $frmin = 1;
my $frmax = 4;
my $frdensity = 2000;
my $fxdensity = 500;
my $fxdepth = 500;
my $levels = 5;
my $ratio = 2;
my $level;
my $t;
my $rmin;
my $rmax;

$t = ($frmax - $frmin);
$rmin = $frmin;
$rmax = $rmin + ($t / $ratio);

print "\t\**********\t $rmin < r < $rmax \t**********\n";

system ("./bifurcation.pl $rmin $rmax $frdensity $fxdensity $fxdepth '$eqn' > logistic_map_$rmin-$rmax");
system ("./make_plot.pl logistic_map_$rmin-$rmax dots");

for ($level = 1; $level < $levels; $level++) {
    $t = ($rmax - $rmin);
    $rmin = $rmax;
    $rmax = $rmin + ($t / $ratio);

    print "\t\**********\t $rmin < r < $rmax \t**********\n";

    system ("./bifurcation.pl $rmin $rmax $frdensity $fxdensity $fxdepth '$eqn' > logistic_map_$rmin-$rmax");
    system ("./make_plot.pl logistic_map_$rmin-$rmax dots");
    
}

$t = ($rmax - $rmin);
$rmin = $rmax;
$rmax = $frmax;

print "\t\**********\t $rmin < r < $rmax \t**********\n";

system ("./bifurcation.pl $rmin $rmax $frdensity $fxdensity $fxdepth '$eqn' > logistic_map_$rmin-$rmax");
system ("./make_plot.pl logistic_map_$rmin-$rmax dots");

exit;