static void Main(string[] args)
{
Console.CancelKeyPress += ResetTerm;
while (true)
{
Console.Write("y = ");
Expression expr = new Expression(Console.ReadLine());
DifferenceQuotient diff = new DifferenceQuotient(expr);
Derivative deriv = new Derivative(diff);
Console.WriteLine("y_0 = {0}", expr.Solve(0));
Console.WriteLine("y_1 = {0}", expr.Solve(1));
Console.WriteLine("y_2 = {0}", expr.Solve(2));
Console.WriteLine("dy/dx_x=0 = {0}", diff.SolveAtPoint(0));
Console.WriteLine("dy/dx_x=1 = {0}", diff.SolveAtPoint(1));
Console.WriteLine("dy/dx_x=2 = {0}", diff.SolveAtPoint(2));
Console.WriteLine("d^2y/dx^2_x=0 = {0}", diff.SolveAtPointRecursive(0, 2));
Console.WriteLine("d^2y/dx^2_x=1 = {0}", diff.SolveAtPointRecursive(1, 2));
Console.WriteLine("d^2y/dx^2_x=2 = {0}", diff.SolveAtPointRecursive(2, 2));
Console.WriteLine("dy/dx = {0}", deriv);
Console.WriteLine();
}
}
public Derivative(DifferenceQuotient differenceQuotient)
{
Terms = new List <Term>();
decimal factor = 1;
for (int i = 1; i < Accuracy; ++i)
{
Terms.Add(new Term {
Coefficient = differenceQuotient.SolveAtPointRecursive(Point, i) * factor,
Exponent = i - 1
});
factor /= i;
}
}