1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
110
111
112
113
114
115
116
117
118
119
|
#!perl
use strict;
use warnings;
use Test::More tests => 7;
BEGIN {
use_ok('Math::Symbolic');
}
if ($ENV{TEST_YAPP_PARSER}) {
require Math::Symbolic::Parser::Yapp;
$Math::Symbolic::Parser = Math::Symbolic::Parser::Yapp->new();
}
use Math::Symbolic::ExportConstants qw/:all/;
my $var = Math::Symbolic::Variable->new();
my $a = $var->new( 'x' => 2 );
my $c = Math::Symbolic::Constant->zero();
my $two = $c->new(2);
print "Vars: x=" . $a->value() . " (Value is optional)\n\n";
my $op = Math::Symbolic::Operator->new();
my $sin;
undef $@;
eval <<'HERE';
$sin = $op->new('cosh', $op->new('*', $two, $a));
HERE
ok( !$@, 'hyperbolic cosine creation' );
my $asin;
undef $@;
eval <<'HERE';
$asin = $op->new('acosh', $op->new('*', $two, $a));
HERE
ok( !$@, 'area hyperbolic cosine creation' );
print "Expression: cosh(2*x) and acosh(2*x)\n\n";
print "prefix notation and evaluation:\n";
undef $@;
eval <<'HERE';
print $sin->to_string('prefix') . "\n\n";
HERE
ok( !$@, 'h. cosine to_string' );
undef $@;
eval <<'HERE';
print $asin->to_string('prefix') . "\n\n";
HERE
ok( !$@, 'area h. cosine to_string' );
print "Now, we derive this partially to x: (prefix again)\n";
my $n_tree = $op->new(
{
type => U_P_DERIVATIVE,
operands => [ $sin, $a ],
}
);
my $n_tree2 = $op->new(
{
type => U_P_DERIVATIVE,
operands => [ $asin, $a ],
}
);
print $n_tree->to_string('prefix') . "\n\n";
print $n_tree2->to_string('prefix') . "\n\n";
print "Now, we apply the derivative to the term: (infix)\n";
my $derived;
undef $@;
eval <<'HERE';
$derived = $n_tree->apply_derivatives();
HERE
ok( !$@, 'h. cosine derivative' );
my $derived2;
undef $@;
eval <<'HERE';
$derived2 = $n_tree2->apply_derivatives();
HERE
ok( !$@, 'area h. cosine derivative' );
print "$derived\n\n";
print "$derived2\n\n";
print "Finally, we simplify the derived term as much as possible:\n";
$derived = $derived->simplify();
$derived2 = $derived2->simplify();
print "$derived\n\n";
print "$derived2\n\n";
print "Now, we do this two more times:\n";
for ( 1 .. 2 ) {
$derived = $op->new(
{
type => U_P_DERIVATIVE,
operands => [ $derived, $a ],
}
)->apply_derivatives()->simplify();
$derived2 = $op->new(
{
type => U_P_DERIVATIVE,
operands => [ $derived2, $a ],
}
)->apply_derivatives()->simplify();
}
print "$derived\n\n";
print "$derived2\n\n";
|
