static void Main(string[] args)
{
decimal n = 1000;
decimal a = 0;
decimal b = 1;
decimal error = Convert.ToDecimal(0.0001);
string function = "ln((1+x)/(1-x))";
PostfixNotationExpression parser = new PostfixNotationExpression();
parser.ToPostfixNotation(function);
parser.CalculatePoint(a, b, n);
List<decimal> Xs = parser.GetXsList();
List<decimal> Ys = parser.GetYsList();
Differentiation differentiationComponent = new Differentiation(Xs, Ys);
List<decimal> Dys = differentiationComponent.CalculateDerivativeByFiniteDifferencies(Xs, Ys);
List<decimal> D2ys = differentiationComponent.CalculateDerivativeByFiniteDifferencies(Xs, Dys);
decimal result = 0;
decimal partitionCount = 0;
/// Trapezoidal Rule
TrapezoidalRule trapezoidalRuleComponent = new TrapezoidalRule();
partitionCount = trapezoidalRuleComponent.CalculatePartitionCount(a, b, error, D2ys.Max());
if (partitionCount == 0)
{
partitionCount = n;
}
parser.CalculatePoint(a, b, partitionCount);
List<decimal> FunctionValues = parser.GetYsList();
result = trapezoidalRuleComponent.Calculate(a, b, partitionCount, FunctionValues);
Console.WriteLine(String.Format("Result: {0}", result));
Console.Read();
}
private void calculateButton_Click(object sender, EventArgs e)
{
try
{
string method = MethodsComboBox.SelectedItem.ToString();
string function = UserFunctionTextBox.Text;
PostfixNotationExpression parser = new PostfixNotationExpression();
parser.ToPostfixNotation(function);
decimal a = Convert.ToDecimal(ParamATextBox.Text);
decimal b = Convert.ToDecimal(ParamBTextBox.Text);
decimal n = 1000;
decimal error = Convert.ToDecimal(ErrorTextBox.Text);
parser.CalculatePoint(a, b, n);
List<decimal> Xs = parser.GetXsList();
List<decimal> Ys = parser.GetYsList();
Differentiation differentiationComponent = new Differentiation(Xs, Ys);
List<decimal> Dys = differentiationComponent.GetFirstDervative();
//= differentiationComponent.CalculateDerivativeByFiniteDifferencies(Xs, Ys);
List<decimal> D2ys = differentiationComponent.GetSecondDervative();
//= differentiationComponent.CalculateDerivativeByFiniteDifferencies(Xs, Dys);
List<decimal> D3ys = differentiationComponent.GetThirdDervative();
//= differentiationComponent.CalculateDerivativeByFiniteDifferencies(Xs, D2ys);
List<decimal> D4ys = differentiationComponent.GetFourthDervative();
//= differentiationComponent.CalculateDerivativeByFiniteDifferencies(Xs, D3ys);
decimal result = 0;
decimal partitionCount = 0;
switch (method)
{
case "Rectangle Method":
{
/// Rectangle Method
RectangleMethod rectangleMethodComponent = new RectangleMethod();
partitionCount = rectangleMethodComponent.CalculatePartitionCount(a, b, error, D2ys.Max());
if (partitionCount == 0)
{
partitionCount = n;
}
parser.CalculatePoint(a, b, partitionCount);
List<decimal> FunctionHalfValues = parser.GetYsHalfList();
result = rectangleMethodComponent.Calculate(a, b, partitionCount, FunctionHalfValues);
MessageBox.Show(result.ToString(), "Rectangle Method");
rectangleMethodComponent.Dispose();
break;
}
case "Trapezoidal Rule":
{
/// Trapezoidal Rule
TrapezoidalRule trapezoidalRuleComponent = new TrapezoidalRule();
partitionCount = trapezoidalRuleComponent.CalculatePartitionCount(a, b, error, D2ys.Max());
if (partitionCount == 0)
{
partitionCount = n;
}
parser.CalculatePoint(a, b, partitionCount);
List<decimal> FunctionValues = parser.GetYsList();
result = trapezoidalRuleComponent.Calculate(a, b, partitionCount, FunctionValues);
MessageBox.Show(result.ToString(), "Trapezoidal Rule");
trapezoidalRuleComponent.Dispose();
break;
}
case "Simpson's Rule":
{
/// Simpson's Rule
SimpsonsRule simpsonsRuleComponent = new SimpsonsRule();
partitionCount = simpsonsRuleComponent.CalculatePartitionCount(a, b, error, D4ys.Max());
if (partitionCount == 0)
{
partitionCount = n;
}
parser.CalculatePoint(a, b, partitionCount);
List<decimal> FunctionHalfValues = parser.GetYsHalfList();
List<decimal> FunctionValues = parser.GetYsList();
result = simpsonsRuleComponent.Calculate(a, b, partitionCount, FunctionValues, FunctionHalfValues);
MessageBox.Show(result.ToString(), "Simpson's Rule");
simpsonsRuleComponent.Dispose();
break;
}
case "All of these methods":
{
Container container = new Container();
RectangleMethod rectangleMethodComponent = new RectangleMethod();
TrapezoidalRule trapezoidalRuleComponent = new TrapezoidalRule();
SimpsonsRule simpsonsRuleComponent = new SimpsonsRule();
container.Add(rectangleMethodComponent, "Rectangle Method");
container.Add(trapezoidalRuleComponent, "Trapezoidal Rule");
container.Add(simpsonsRuleComponent, "Simpson's Rule");
string message = "";
/// Rectangle Method
partitionCount = rectangleMethodComponent.CalculatePartitionCount(a, b, error, D2ys.Max());
if (partitionCount == 0)
{
partitionCount = n;
}
parser.CalculatePoint(a, b, partitionCount);
List<decimal> FunctionHalfValues = parser.GetYsHalfList();
result = rectangleMethodComponent.Calculate(a, b, partitionCount, FunctionHalfValues);
message += rectangleMethodComponent.Site.Name + ":\n" + result.ToString() + "\n";
/// Trapezoidal Rule
partitionCount = trapezoidalRuleComponent.CalculatePartitionCount(a, b, error, D2ys.Max());
if (partitionCount == 0)
{
partitionCount = n;
}
parser.CalculatePoint(a, b, partitionCount);
List<decimal> FunctionValues = parser.GetYsList();
result = trapezoidalRuleComponent.Calculate(a, b, partitionCount, FunctionValues);
message += trapezoidalRuleComponent.Site.Name + ":\n" + result.ToString() + "\n";
/// Simpson's Rule
partitionCount = simpsonsRuleComponent.CalculatePartitionCount(a, b, error, D4ys.Max());
if (partitionCount == 0)
{
partitionCount = n;
}
parser.CalculatePoint(a, b, partitionCount);
FunctionHalfValues = parser.GetYsHalfList();
FunctionValues = parser.GetYsList();
result = simpsonsRuleComponent.Calculate(a, b, partitionCount, FunctionValues, FunctionHalfValues);
message += simpsonsRuleComponent.Site.Name + ":\n" + result.ToString();
MessageBox.Show(message, "All methods");
container.Dispose();
break;
}
default:
{
MessageBox.Show("Choose method", "Warning");
break;
}
}
differentiationComponent.Dispose();
parser.Dispose();
}
catch (Exception exception)
{
MessageBox.Show(exception.ToString());
}
}