UNPKG

algebrite

Version:

Computer Algebra System in Coffeescript

github.com/davidedc/Algebrite

davidedc/Algebrite

108 lines (90 loc) • 1.68 kB
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
#----------------------------------------------------------------------------- # # Author : philippe.billet@noos.fr # # Gamma function gamma(x) # #----------------------------------------------------------------------------- Eval_gamma = -> push(cadr(p1)) Eval() gamma() gamma = -> save() gammaf() restore() gammaf = -> # double d p1 = pop() if (isrational(p1) && MEQUAL(p1.q.a, 1) && MEQUAL(p1.q.b, 2)) if evaluatingAsFloats push_double(Math.PI) else push_symbol(PI) push_rational(1,2) power() return if (isrational(p1) && MEQUAL(p1.q.a, 3) && MEQUAL(p1.q.b, 2)) if evaluatingAsFloats push_double(Math.PI) else push_symbol(PI) push_rational(1,2) power() push_rational(1,2) multiply() return # if (p1->k == DOUBLE) { # d = exp(lgamma(p1.d)) # push_double(d) # return # } if (isnegativeterm(p1)) if evaluatingAsFloats push_double(Math.PI) else push_symbol(PI) push_integer(-1) multiply() if evaluatingAsFloats push_double(Math.PI) else push_symbol(PI) push(p1) multiply() sine() push(p1) multiply() push(p1) negate() gamma() multiply() divide() return if (car(p1) == symbol(ADD)) gamma_of_sum() return push_symbol(GAMMA) push(p1) list(2) return gamma_of_sum = -> p3 = cdr(p1) if (isrational(car(p3)) && MEQUAL(car(p3).q.a, 1) && MEQUAL(car(p3).q.b, 1)) push(cadr(p3)) push(cadr(p3)) gamma() multiply() else if (isrational(car(p3)) && MEQUAL(car(p3).q.a, -1) && MEQUAL(car(p3).q.b, 1)) push(cadr(p3)) gamma() push(cadr(p3)) push_integer(-1) add() divide() else push_symbol(GAMMA) push(p1) list(2) return