use strict;
use warnings;

$, = "  ";
$\=$/;
my $h = 0.1;
my $aedge = 0;
my $bedge = 1;
my $n = int(($bedge - $aedge)/$h) + 1;
my (@lambda, @mu, @y, @A, @B);

#$mu[0]=$lambda[0]=0;

$A[0] = 1;
$B[0] = 0;
my $xi = $aedge;
for my $i(1..$n-1){
        $xi += $h;
        $A[$i] = 2 + $h*$h;
        $B[$i] = $h*$h*(92 + 45*$xi*(1-$xi));
}
$A[$n] = 1 + $h;
$B[$n] = -$h*(2*exp(1) + 43);

$lambda[1] = 0;
$mu[1]     = 0;

for my $i(1..$n){
        $lambda[$i+1] = 1/($A[$i] - $lambda[$i]);
        $mu[$i+1] = ($mu[$i]-$B[$i])/($A[$i]-$lambda[$i]);
}

$y[$n] = $mu[$n+1];
for my $i(reverse(1..$n)){
        print ($i);
        $y[$i-1] = $lambda[$i]*$y[$i] + $mu[$i];
}

open F, ">results.txt";
$xi = 0;
for my $i(0..$n){
        print F "($xi;$y[$i])";
        $xi += $h;
}

print @B;
print @lambda;
print @mu