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
|
/* PIA - Lenguajes Modernos de Programación
* FACULTAD DE CIENCIAS FÍSICO MATEMÁTICAS
* Luis Sebastián Martínez Vega - LCC */
#include "include/differentiator.hpp"
#include "include/expressions.hpp"
Differentiator::Differentiator() {
differential = NULL;
ptr = NULL;
}
// Regla del constante. f(x)=c, f'(x)=0.
Expression *Differentiator::visit(Literal *expr) {
return new Literal(0);
}
// Regla de la variable. f(x)=x, f'(x)=1
Expression *Differentiator::visit(Variable *expr) {
return new Literal(1);
}
// Reglas de funciones trigonométricas.
Expression *Differentiator::visit(Function *expr) {
Expression *diff;
switch (expr->getFunctionName()) {
case i_sin:
diff = new MultiplicationExpression(new Function(expr->getArg()->copy(), i_cos), expr->getArg()->diff(this));
return diff;
case i_cos:
diff = new MultiplicationExpression(new NegationExpression(new Function(expr->getArg()->copy(), i_sin)), expr->getArg()->diff(this));
return diff;
case i_tan:
diff = new MultiplicationExpression(new PowerExpression(new Function(expr->getArg()->copy(), i_sec), new Literal(2)),
expr->getArg()->diff(this));
return diff;
case i_csc:
diff = new MultiplicationExpression(new MultiplicationExpression(new NegationExpression(new Function(expr->getArg()->copy(), i_csc)),
new Function(expr->getArg()->copy(), i_ctg)), expr->getArg()->diff(this));
return diff;
case i_sec:
diff = new MultiplicationExpression(new MultiplicationExpression(new Function(expr->getArg()->copy(), i_sec),
new Function(expr->getArg()->copy(), i_tan)), expr->getArg()->diff(this));
return diff;
case i_ctg:
diff = new MultiplicationExpression(new NegationExpression(new PowerExpression(new Function(expr->getArg()->copy(), i_csc),
new Literal(2))), expr->getArg()->copy());
return diff;
}
}
Expression *Differentiator::visit(NegationExpression *expr) {
return new NegationExpression(expr->diff(this));
}
Expression *Differentiator::visit(AddExpression *expr) {
return new AddExpression(expr->getLeft()->diff(this), expr->getRight()->diff(this));
}
Expression *Differentiator::visit(SubExpression *expr) {
return new SubExpression(expr->getLeft()->diff(this), expr->getRight()->diff(this));
}
// Regla de la potencia. f(x)=x^c, f'(x)=c*x^(c-1)}*f'(x).
Expression *Differentiator::visit(PowerExpression *expr) {
Literal *c = dynamic_cast<Literal *>(expr->getRight());
return new MultiplicationExpression(new MultiplicationExpression(c->copy(), new PowerExpression(expr->getLeft()->copy(), new Literal(c->getValue() - 1))),
expr->getLeft()->diff(this));
}
|