33_trig.pm

   1 # trigonometry for DCT, by Shiar
   2
   3 # 1.11.1 200410282330 - cardial mode setting; rad/deg to switch to radians/degrees
   4 #                     - convert from/to radians for trig commands if rad mode set
   5 # 1.11.0 200410152320 - a?(sin|cos|tan)h? actions from math; links in main submenu trig
   6
   7 use strict;
   8 use warnings;
   9
  10 my $pi = atan2(1, 1) * 4;
  11
  12 $set{card} = 1;  # degrees radians grades
  13
  14 %action = (
  15         %action,
  16
  17         'pi'   => [0, sub { $pi }], # pi constant
  18
  19         'deg'  => [-1, sub { $set{card} = 1 }], # set degrees
  20         'rad'  => [-1, sub { $set{card} = 2 }], # set radians
  21
  22         # trigonometric
  23         'sin'  => [1, sub { sin $_[0] }], # sine
  24         'asin' => [1, sub { atan2($_[0], sqrt(1 - $_[0]*$_[0])) }], # inverse sine
  25         'cos'  => [1, sub { cos $_[0] }], # cosine
  26         'acos' => [1, sub { atan2(sqrt(1 - $_[0]*$_[0]), $_[0]) }], # inverse cosine
  27         'tan'  => [1, sub { sin($_[0]) / cos($_[0]) }], # tangent
  28         'atan' => [1, sub { atan2($_[0], 1) }], # arctangent
  29
  30         # hyperbolic
  31         'sinh' => [1, sub { (exp($_[0]) - exp(-$_[0])) / 2 }], # hyperbolic sine
  32         'cosh' => [1, sub { (exp($_[0]) + exp(-$_[0])) / 2 }], # hyperbolic cosine
  33         'tanh' => [1, sub { (exp($_[0]) - exp(-$_[0])) / (exp($_[0]) + exp(-$_[0])) }], # hyperbolic tangent (sinh/cosh)
  34         'asinh'=> [1, sub { log(sqrt($_[0]**2+1) + $_[0]) }], # inverse hyperbolic sine
  35         'acosh'=> [1, sub { log(sqrt($_[0]**2-1) + $_[0]) }], # inverse hyperbolic cosine
  36         'atanh'=> [1, sub { log((1+$_[0]) / (1-$_[0])) / 2 }], # inverse hyperbolic tangent
  37 ); # action
  38
  39 push @{$hook{preaction}}, sub {
  40         return unless $set{card}==2;
  41         # convert user input from radians if necessary
  42         $stack[0] /= 360/$pi if $_[1] =~ /^(?:sin|cos|tan)h?$/;
  43 }; # preaction
  44 push @{$hook{postaction}}, sub {
  45         return unless $set{card}==2;
  46         # convert command output to radians if necessary
  47         $stack[0] *= 360/$pi if $_[1] =~ /^a(?:sin|cos|tan)h?$/;
  48 }; # postaction
  49
  50 addmenu(["main", 0], "trig", #todo: in math, not in main
  51         qw(sin cos tan asin acos atan),
  52         qw(sinh cosh tanh asinh acosh atanh),
  53         qw(expm lnp1),
  54 );
  55
  56 return {
  57         author  => "Shiar",
  58         title   => "trigonometry",
  59         version => "1.11.1",
  60 };
  61