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
|
#!/usr/bin/perl
use strict;
use warnings;
use lib '../lib/';
use Math::Symbolic qw/:all/;
use Math::Symbolic::VectorCalculus qw/:all/;
my @gradient = grad 'x*y + 2*z*x - y^2';
@gradient = map { $_->apply_derivatives()->simplify() } @gradient;
print "$_\n" foreach @gradient;
print "\n\n";
my @funcs = ( 'x*y+2*z*x-y^2', 'y+x+z', 'x*y*z' );
my $div = div @funcs;
print $div->apply_derivatives()->simplify();
print "\n\n";
my @rot = rot @funcs;
@rot = map { $_->apply_derivatives()->simplify() } @rot;
print "$_\n" foreach @rot;
print "\n\n";
my @matrix = Jacobi @funcs;
print "[\n";
foreach my $func (@matrix) {
print " [\n";
foreach my $var (@$func) {
$var = $var->apply_derivatives()->simplify();
print " $var\n";
}
print " ]\n";
}
print "]\n";
|
