#include <vector>
#include <nlopt.hpp>
#include <cmath>
#include <algorithm>
#include <iostream>

double f (const std::vector<double> &X, std::vector<double> &grad, void *params);

std::vector<double> argmax(double r, double T, double S0, double V, double K, int d);