Simpson method