我正在尝试像这样解决这个指数方程:
my ($l,$r);
for (1 .. 100) -> $x {
$l = $x * e ** $x;
$r = 5 * (e ** $x - 1);
say $x if $l == $r;
}
但是它不起作用。如何以简明易懂的方式全面解决此问题?
我正在尝试像这样解决这个指数方程:
my ($l,$r);
for (1 .. 100) -> $x {
$l = $x * e ** $x;
$r = 5 * (e ** $x - 1);
say $x if $l == $r;
}
但是它不起作用。如何以简明易懂的方式全面解决此问题?
抱歉回答有些重复。
但这里是一种完全不同且更简单的方法,使用Raku解决了问题。
(它可能可以更加优雅地表达。)
#!/usr/bin/env raku
sub solver ($equ, $acc, $lower0, $upper0) {
my Real $lower = $lower0;
my Real $upper = $upper0;
my Real $middle = ($lower + $upper) / 2;
# zero must be in between
sign($equ($lower)) != sign($equ($upper)) || die 'Bad interval!';
for ^$acc { # accuracy steps
if sign($equ($lower)) != sign($equ($middle))
{ $upper = $middle }
else
{ $lower = $middle }
$middle = ($upper + $lower) / 2;
}
return $middle;
}
my $equ = -> $x { $x * e ** $x - 5 * (e ** $x - 1) }; # left side - right side
my $acc = 64; # 64 bit accuracy
my Real $lower = 1; # start search here
my Real $upper = 100; # end search here
my $solution = solver $equ, $acc, $lower, $upper;
say 'result is ', $solution;
say 'Inserted in equation calculates to ', $equ($solution), ' (hopefully nearly zero)'
对于Perl 5,有Math::GSL::Roots - 找到任意一维函数的根
https://metacpan.org/pod/Math::GSL::Roots
Raku支持使用Perl 5代码或直接访问GSL C库,对吗?
$fspec = sub {
my ( $x ) = shift;
# here the function has to be inserted in the format
# return leftside - rightside;
return ($x + $x**2) - 4;
};
gsl_root_fsolver_alloc($T); # where T is the solver algorithm, see link for the 6 type constants, e.g. $$gsl_root_fsolver_brent
gsl_root_fsolver_set( $s, $fspec, $x_lower, $x_upper ); # [$x_lower; $x_upper] is search interval
gsl_root_fsolver_iterate($s);
gsl_root_fsolver_iterate($s);
gsl_root_fsolver_iterate($s);
gsl_root_fsolver_iterate($s);
gsl_root_fsolver_iterate($s);
my $result = gsl_root_fsolver_root($s);
gsl_root_fsolver_free (s);
如果函数的导数可用,可以使用增强算法(gsl_root_fdfsolver_*)。
有关一般用法,请参见https://www.gnu.org/software/gsl/doc/html/roots.html#examples。
Math::Symbolic
(https://github.com/raydiak/Math-Symbolic)。看起来.new
和.isolate
可能会解决问题。 - raiph