public override void CreateDerivativeTree(BaseNode parent, bool isLeft = true)
{
MultiplicationNode multiplicationNode1 = new MultiplicationNode(Plotter.CloneTree(this.left), Plotter.CloneTree(this.right), null);
MultiplicationNode multiplicationNode2 = new MultiplicationNode(Plotter.CloneTree(this.left), Plotter.CloneTree(this.right), null);
SubstractionNode substraction = new SubstractionNode(multiplicationNode1, multiplicationNode2, null);
PowerNode power = new PowerNode(Plotter.CloneTree(this.right), new NumberNode(null, 2), null);
DivisionNode node = new DivisionNode(substraction, power, null);
if (parent != null)
{
if (isLeft)
{
parent.left = node;
}
else
{
parent.right = node;
}
}
node.left.left.left.CreateDerivativeTree(node.left.left);
node.left.right.right.CreateDerivativeTree(node.left.right, false);
//Plotter.SetDerivativeRoot (node);
SetDerivativeRoot(node);
}
/// <summary>
/// Produces a Lagrange Item, which looks like this: (x-Xj)/(Xterm - Xj), where term != j
/// Example: ((x-X1)(x-X2)...(x-Xn)) / ((Xterm - x1)(Xterm - x2)...(Xterm - Xn))
///
/// </summary>
/// <param name="points">Given set of points</param>
/// <param name="term">For which term you want to create a so-called Lagrange item</param>
/// <returns>Lagrange item of the type MultiplicationNode</returns>
public static MultiplicationNode ProduceLagrange(DataPoint[] points, int term)
{
MultiplicationNode item = new MultiplicationNode(null, null, null);
item.left = new NumberNode(null, 1);
for (int i = 0; i < points.Length; i++)
{
if (i == term)
{
continue;
}
SubstractionNode enumerator = new SubstractionNode(
new BasicFunctionXNode(""), // left
new NumberNode(null, points[i].X), // right
null // parent
);
NumberNode denominator = new NumberNode(null, points[term].X - points[i].X);
DivisionNode division = new DivisionNode(enumerator, denominator, null);
item.LagrangePutToRightNode(division);
}
return(item);
}
/// <summary>
/// Creates a MacLaurien series based off the function
/// </summary>
/// <param name="mcLaurienRoot">Where the McLaurien Series will be outputted</param>
/// <param name="order">Nth order of a series</param>
public void CreateMcLaurienSeries(out BaseNode mcLaurienRoot, int order = 5)
{
if (derivativeRoot == null)
{
CreateDerivativeTree(); derivativeRoot.Simplify();
}
// we made sure that there is a derivative
BaseNode myDerivative = Plotter.CloneTree(derivativeRoot);
myDerivative.derivativeRoot = myDerivative;
double[] values = new double[order + 1]; // values for functions (f(0), derivative of f(0), second derivative of f(0), etc..)
values[0] = root.Calculate(0); // set up a value for f(0)
if (values.Length >= 2)
{
values[1] = derivativeRoot.Calculate(0); // set up a value of the first derivative of f(0)
}
if (values.Length >= 3)
{
for (int i = 2; i < values.Length; i++)
{
myDerivative.CreateDerivativeTree(null);
values[i] = myDerivative.derivativeRoot.Calculate(0);
myDerivative = myDerivative.derivativeRoot;
}
}
List <BaseNode> mcLaurienItems = new List <BaseNode> ();
SumNode result = new SumNode(null, null, null);
result.left = new NumberNode(null, values[0]);
for (int i = 1; i < values.Length; i++)
{
DivisionNode item = new DivisionNode(null, null, null);
FactorialNode denominator = new FactorialNode(new NumberNode(null, i), null); // not sure about this line
MultiplicationNode numerator = new MultiplicationNode(
new NumberNode(null, values[i]),
new PowerNode(new BasicFunctionXNode("", null), new NumberNode(null, i), null), null
);
item.left = numerator;
item.right = denominator;
mcLaurienItems.Add(item);
}
foreach (var item in mcLaurienItems)
{
result.PutToRightNode(item);
}
mcLaurienRoot = result;
}
public void LagrangePutToRightNode(BaseNode node)
{
if (this.right == null)
{
this.right = node;
}
else
{
MultiplicationNode multiplication = new MultiplicationNode(this.right, node, null);
this.right = multiplication;
}
}
public override void CreateDerivativeTree(BaseNode parent, bool isLeft = true)
{
CosNode cosNode = new CosNode(Plotter.CloneTree(this.left), null);
MultiplicationNode node = new MultiplicationNode(Plotter.CloneTree(this.left), cosNode, null);
if (parent != null)
{
if (isLeft)
{
parent.left = node;
}
else
{
parent.right = node;
}
}
node.left.CreateDerivativeTree(node);
//Plotter.SetDerivativeRoot (node);
SetDerivativeRoot(node);
}
/// <summary>
/// Get the McLaurien series based off already given values for derivatives
/// </summary>
/// <param name="mcLaurienRoot">Where the McLaurien series will be outputted</param>
/// <param name="order">Order of the McLaurien</param>
public void CreateMcLaurienSeriesByLimits(out BaseNode mcLaurienRoot, int order)
{
List <double> values = new List <double> ();
for (int i = 0; i <= order; i++)
{
var value = ProcessNthDerivative_Quotient(0, i, root);
values.Add(value);
}
List <BaseNode> mcLaurienItems = new List <BaseNode> ();
SumNode result = new SumNode(null, null, null);
result.left = new NumberNode(null, values[0]);
for (int i = 1; i < values.Count; i++)
{
DivisionNode item = new DivisionNode(null, null, null);
FactorialNode denominator = new FactorialNode(new NumberNode(null, i), null); // not sure about this line
MultiplicationNode numerator = new MultiplicationNode(
new NumberNode(null, values[i]),
new PowerNode(new BasicFunctionXNode("", null), new NumberNode(null, i), null), null
);
item.left = numerator;
item.right = denominator;
mcLaurienItems.Add(item);
}
foreach (var item in mcLaurienItems)
{
result.PutToRightNode(item);
}
mcLaurienRoot = result;
}
public BaseNode CreatePolynomialThroughPoints(DataPoint[] points)
{
/* Let X be the number of points user has selected
* Then the polynomial is going to be of the degree (X - 1)
*
* Example: X = 5 => ax^4 + bx^3 + cx^2 + dx + e
*/
SumNode polynomial = new SumNode(null, null, null);
var firstDivision = ProduceLagrange(points, 0);
var firstMultiplication = new MultiplicationNode(firstDivision, new NumberNode(null, points[0].Y), null);
polynomial.left = firstMultiplication;
for (int j = 1; j < points.Length; j++)
{
var division = ProduceLagrange(points, j);
MultiplicationNode node = new MultiplicationNode(division, new NumberNode(null, points[j].Y), null);
polynomial.PutToRightNode(node);
}
return(polynomial);
}
public override void CreateDerivativeTree(BaseNode parent, bool isLeft = true)
{
// derivative of ln(x) = 1/x
DivisionNode division = new DivisionNode(new NumberNode(null, 1), Plotter.CloneTree(this.left), null);
MultiplicationNode multiplication = new MultiplicationNode(division, Plotter.CloneTree(this.left), null);
if (parent != null)
{
if (isLeft)
{
parent.left = multiplication;
}
else
{
parent.right = multiplication;
}
}
multiplication.right.CreateDerivativeTree(multiplication, false);
//Plotter.SetDerivativeRoot (multiplication);
SetDerivativeRoot(multiplication);
}
public override void CreateDerivativeTree(BaseNode parent, bool isLeft = true)
{
// d(e^f(x))/dx = d(f(x))/dx * e^f(x)
// MIGHT HAVE SOME BUGS, CHECK IT
MultiplicationNode multiplication = new MultiplicationNode(Plotter.CloneTree(this.left), Plotter.CloneTree(this), null);
if (parent != null)
{
if (isLeft)
{
parent.left = multiplication;
}
else
{
parent.right = multiplication;
}
}
multiplication.left.CreateDerivativeTree(multiplication);
SetDerivativeRoot(multiplication);
}
/// <summary>
/// Clones a specified tree based on a given node 'root'
/// </summary>
/// <param name="root"></param>
/// <returns></returns>
public static BaseNode CloneTree(BaseNode root)
{
if (root == null)
{
return(null);
}
BaseNode newNode = null;
if (root is SubstractionNode)
{
newNode = new SubstractionNode(root.value);
}
else if (root is MultiplicationNode)
{
newNode = new MultiplicationNode(root.value);
}
else if (root is SumNode)
{
newNode = new SumNode(root.value);
}
else if (root is DivisionNode)
{
newNode = new DivisionNode(root.value);
}
else if (root is NumberNode)
{
newNode = new NumberNode(null, (root as NumberNode).RealValue);
}
else if (root is BasicFunctionXNode)
{
newNode = new BasicFunctionXNode(root.value);
}
else if (root is SinNode)
{
newNode = new SinNode(root.value);
}
else if (root is CosNode)
{
newNode = new CosNode(root.value);
}
else if (root is PowerNode)
{
newNode = new PowerNode(root.value);
}
else if (root is ExponentNode)
{
newNode = new ExponentNode(root.value);
}
else if (root is LnNode)
{
newNode = new LnNode(root.value);
}
else if (root is FactorialNode)
{
newNode = new FactorialNode(root.value);
}
newNode.left = CloneTree(root.left);
newNode.right = CloneTree(root.right);
return(newNode);
}
/// <summary>
/// Creates tree based on an input string and root node
/// </summary>
/// <param name="s"></param>
/// <param name="baseNode"></param>
public void CreateTree(string s, BaseNode baseNode)
{
// if the string is empty, we don't do anything. This is the base case to leave the recursion
if (s == string.Empty)
{
return;
}
// if it's 's', or '+', or whatever, we create a dedicated class (watch first case to see the logic)
if (s[0] == 's')
{
SinNode node = new SinNode(s, baseNode); // dedicated class
baseNode.Insert(node); // we insert it to the current head node
CreateTree(node.value, node); // we change the head node to the newly created one
}
else if (s[0] == 'c')
{
CosNode node = new CosNode(s, baseNode);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == '*')
{
// same as in the first 'if'
MultiplicationNode node = new MultiplicationNode(s, baseNode);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == '+')
{
// same as in the first 'if'
SumNode node = new SumNode(s, baseNode);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == '/')
{
DivisionNode node = new DivisionNode(s, baseNode);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == '-' && !(s[1] >= '0' && s[1] <= '9'))
{
SubstractionNode node = new SubstractionNode(s, baseNode);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == 'l')
{
LnNode node = new LnNode(s, baseNode);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == '^')
{
PowerNode node = new PowerNode(s, baseNode);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == 'e')
{
ExponentNode node = new ExponentNode(s, baseNode);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == '!')
{
FactorialNode node = new FactorialNode(s, baseNode);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == 'p' || (s[0] >= '0' && s[0] <= '9'))
{
// stuff below just parses number
string toParseIntoNumber = string.Empty;
int counter = 0;
if (s[0] == 'p')
{
toParseIntoNumber = "p";
}
else
{
while ((s[counter] >= '0' && s[counter] <= '9') || s[counter] == '.')
{
toParseIntoNumber += s[counter];
counter++;
}
}
if (toParseIntoNumber.Contains('.'))
{
toParseIntoNumber = toParseIntoNumber.Replace('.', ',');
}
string @newS = string.Empty;
for (int i = (s[0] == 'p' ? 1 : counter); i < s.Length; i++)
{
newS += s[i];
}
// same stuff as in the first 'if'
NumberNode node = new NumberNode(newS, baseNode, toParseIntoNumber);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == '-' && (s[1] >= '0' && s[1] <= '9'))
{
// negative number
s = Plotter.GetStringFromIndex(s, 1);
string toParseIntoNumber = string.Empty;
int counter = 0;
if (s[0] == 'p')
{
toParseIntoNumber = "p";
}
else
{
do
{
toParseIntoNumber += s[counter];
counter++;
} while (counter < s.Length && ((s[counter] >= '0' && s[counter] <= '9') || s[counter] == '.'));
}
if (toParseIntoNumber.Contains('.'))
{
toParseIntoNumber = toParseIntoNumber.Replace('.', ',');
}
string @newS = string.Empty;
for (int i = (s[0] == 'p' ? 1 : counter); i < s.Length; i++)
{
newS += s[i];
}
NumberNode node = new NumberNode(newS, baseNode, "-" + toParseIntoNumber);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == 'x')
{
// same as in the first 'if'
BasicFunctionXNode node = new BasicFunctionXNode(s, baseNode);
baseNode.Insert(node);
CreateTree(node.value, node);
}
else if (s[0] == '(' || s[0] == ' ')
{
s = GetStringFromIndex(s, 1); // practically delete that ( or ' '
CreateTree(s, baseNode);
}
else if (s[0] == ')')
{
// count how many times ')' appears, let this number be 'i', then our head node is gonna go 'i' levels up
int i = 0;
while (s[i] == ')' && (s[i] != ',' || s[i] != ' '))
{
i++;
if (i == s.Length)
{
break;
}
}
for (int j = 0; j < i; j++)
{
if (baseNode.parent != null)
{
baseNode = baseNode.parent;
}
else
{
throw new Exception("Eror in your input");
}
}
s = GetStringFromIndex(s, i);
CreateTree(s, baseNode);
}
else if (s[0] == ',')
{
if (baseNode.parent == null)
{
throw new Exception("Error in your input");
}
// go one level up
baseNode = baseNode.parent;
s = GetStringFromIndex(s, 1);
CreateTree(s, baseNode);
}
}
public override void CreateDerivativeTree(BaseNode parent, bool isLeft = true)
{
if (this.right is NumberNode && this.left is BasicFunctionXNode)
{
var lesser = (right as NumberNode).RealValue - 1;
BasicFunctionXNode x = new BasicFunctionXNode("", null);
MultiplicationNode multiplication = new MultiplicationNode(new NumberNode(null, (right as NumberNode).RealValue),
new PowerNode(x, new NumberNode(null, lesser), null), null);
if (parent != null)
{
if (isLeft)
{
parent.left = multiplication;
}
else
{
parent.right = multiplication;
}
}
SetDerivativeRoot(multiplication);
return;
}
else
{
if (this.right is NumberNode && this.left is NumberNode)
{
// if both this.left and this.right are numbers, return 0 for its just a number and it's anyway gon be 0
NumberNode node = new NumberNode(parent, 0);
if (parent != null)
{
if (isLeft)
{
parent.left = node;
}
else
{
parent.right = node;
}
}
//Plotter.SetDerivativeRoot (node);
SetDerivativeRoot(node);
return;
}
else if (this.right is NumberNode && !(this.left is NumberNode))
{
// f(x) ^ (some number)
// if left one some function
double nMinus1 = ((NumberNode)this.right).RealValue - 1;
var value = ((NumberNode)this.right).RealValue;
if (value == 1)
{
var node = Plotter.CloneTree(this);
if (parent != null)
{
if (isLeft)
{
parent.left = node;
}
else
{
parent.right = node;
}
}
node.left.CreateDerivativeTree(node);
SetDerivativeRoot(node);
return;
}
PowerNode power = new PowerNode(Plotter.CloneTree(this.left), new NumberNode(null, nMinus1), null);
MultiplicationNode multiplication = new MultiplicationNode(new NumberNode(null, value), Plotter.CloneTree(this.left), null);
// if the f(x) is more complicated than just 'x', we do additional calculation
if (!(multiplication.right is BasicFunctionXNode))
{
MultiplicationNode node = new MultiplicationNode(multiplication, Plotter.CloneTree(this.left), parent);
node.right.CreateDerivativeTree(multiplication, false);
if (parent != null)
{
if (isLeft)
{
parent.left = node;
}
else
{
parent.right = node;
}
}
//Plotter.SetDerivativeRoot (node);
SetDerivativeRoot(node);
return;
}
multiplication.parent = parent;
if (parent != null)
{
if (isLeft)
{
parent.left = multiplication;
}
else
{
parent.right = multiplication;
}
}
//Plotter.SetDerivativeRoot (multiplication);
SetDerivativeRoot(multiplication);
return;
}
else if (!(this.right is NumberNode) && (this.left is NumberNode))
{
// (some number) ^ f(x)
var value = ((NumberNode)this.left).RealValue;
if (this.right is BasicFunctionXNode)
{
// simple function
PowerNode power = new PowerNode(new NumberNode(null, value), new BasicFunctionXNode(""), null);
LnNode ln = new LnNode(new NumberNode(null, value), null);
MultiplicationNode node = new MultiplicationNode(power, ln, parent);
if (parent != null)
{
if (isLeft)
{
parent.left = node;
}
else
{
parent.right = node;
}
}
//Plotter.SetDerivativeRoot (node);
SetDerivativeRoot(node);
return;
}
else
{
// function is more complicated
PowerNode power = new PowerNode(new NumberNode(null, value), this.right, null);
LnNode ln = new LnNode(new NumberNode(null, value), null);
MultiplicationNode multiplication = new MultiplicationNode(power, ln, parent);
MultiplicationNode node = new MultiplicationNode(multiplication, Plotter.CloneTree(this.right), parent);
node.right.CreateDerivativeTree(node, false);
if (parent != null)
{
if (isLeft)
{
parent.left = node;
}
else
{
parent.right = node;
}
}
//Plotter.SetDerivativeRoot (node);
SetDerivativeRoot(node);
return;
}
}
else if (!(this.right is NumberNode) && !(this.left is NumberNode))
{
// neither is a number
// CASE: f(x) ^ g(x)
// d(f(x) ^ g(x))/dx = e^(g(x)*ln(f(x)) * d((g(x)*f(x)))/dx )
// this.left = f(x), this.right = g(x)
LnNode lnFx = new LnNode(Plotter.CloneTree(this.left), null); // create ln(f(x))
MultiplicationNode multiplication = new MultiplicationNode(Plotter.CloneTree(this.right), lnFx, null); // create g(x)*ln(f(x))
PowerNode ePower = new PowerNode(new NumberNode(null, Math.E), multiplication, null); // create e^(g(x)*ln(f(x)))
MultiplicationNode derivativeOfMultiplication = new MultiplicationNode(Plotter.CloneTree(multiplication.left), Plotter.CloneTree(multiplication.right), null); // do the derivative of g(x)*ln(f(x))
MultiplicationNode node = new MultiplicationNode(ePower, derivativeOfMultiplication, parent); // put it all together
node.right.CreateDerivativeTree(node, false); // take a derivative
if (parent != null)
{
if (isLeft)
{
parent.left = node;
}
else
{
parent.right = node;
}
}
//Plotter.SetDerivativeRoot (node);
SetDerivativeRoot(node);
return;
}
}
}
