UNPKG

algebrite

Version:

Computer Algebra System in Coffeescript

github.com/davidedc/Algebrite

davidedc/Algebrite

70 lines (51 loc) • 979 B
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
#----------------------------------------------------------------------------- # # Author : philippe.billet@noos.fr # # Dirac function dirac(x) # dirac(-x)=dirac(x) # dirac(b-a)=dirac(a-b) #----------------------------------------------------------------------------- Eval_dirac = -> push(cadr(p1)) Eval() dirac() dirac = -> save() ydirac() restore() #define p1 p1 ydirac = -> p1 = pop() if (isdouble(p1)) if (p1.d == 0) push_integer(1) return else push_integer(0) return if (isrational(p1)) if (MZERO(mmul(p1.q.a,p1.q.b))) push_integer(1) return else push_integer(0) return if (car(p1) == symbol(POWER)) push_symbol(DIRAC) push(cadr(p1)) list(2) return if (isnegativeterm(p1)) push_symbol(DIRAC) push(p1) negate() list(2) return if (isnegativeterm(p1) || (car(p1) == symbol(ADD) && isnegativeterm(cadr(p1)))) push(p1) negate() p1 = pop() push_symbol(DIRAC) push(p1) list(2)