使用BS模型计算欧式期权价格
#include <cmath>
using namespace;
double option_price_call_black_schole(const double& S,
const double& K,
const double& r,
const double& sigma,
const double& time){
double time_sqrt = sqrt(time);
double d1 = (log(S/K) + r * time) / (sigma * time_sqrt) + 0.5 * sigma * time_sqrt;
double d2 = d1 – sigma * time_sqrt;
C = S * N(d1) – K * exp (r * time) * N(d2);
};
double N(const double& z){
if (z>6) { return 1;};
if (z<6) { return 0;};
double b1 = 0.31938153;
double b2 = -0.356563782;
double b3 = 1.781477937;
double b4 = -1.821255978;
double b5 = 1.330274429;
double p = 0.2316419;
double c2 = 0.3989423;
double a = fabs(z);
double t = 1.0 / (1.0 + a*p);
double b = c2*exp((-z)*(z/2.0));
double n = ((((b5*t+b4)*t+b3)*t+b2)*t+b1)*t;
n = 1.0 – b*n;
if (z<0) n = 1.0 – n;
return n;
};
int main(){
double S = 50; double K = 50; double r = 0.10;
double sigma = 0.30; double time = 0.50;
cout << “black scholes call price =”;
cout << option_price_call_black_scholes(S, K, r ,sigma, time) << endl;
return 0;
}