Подключение dll в C#

#include <complex>
#include <cmath>
using namespace std;
extern "C" __declspec(dllexport) double * count(double * data)
{
//Порядок входных данных - Ti, kp, Ak, Tk, d
const double eps = 1E-4;
double beta = 0.24, n = 2, Fob, Fob_op;

double * data_n = new double[6];
double const a = 1;
double * count;
count = new double[6];
double const Pi = 4 * atan(a);
double omega_k, Fre, Frs_k, Rrs_k, Rr_1;
double d = data[4], hys = data[4] / 2;
double x_1;
omega_k = (2 * Pi) / data[3];
Fre = -asin(hys / data[2]);
Frs_k = - Pi - Fre;
Rrs_k = (data[2] * Pi) / (d * 4);
complex<double> Wrs;
const complex<double> i(0,1);
Wrs = Rrs_k * exp(Frs_k * i);
complex<double> Wr1(1, -(1 / (omega_k * data[0])));
Rr_1 = abs(Wr1);
complex<double> Wob;
double Rob;
Wob = Wrs / Wr1;
Rob = abs(Wob);
Fob = arg(Wob);
double left = 0, right = 2;
double center = left;
while (right - left > eps)
{
if ((beta * center + atan(center) + atan(center * n) + Fob) * (beta * right + atan(right) + atan(right * n) + Fob) < 0)
{
left = center;
}
else
{
right = center;
}
center = (left + right) / 2;
}
x_1 = center;

double T1, T2, tau;
T1 = x_1 / omega_k;
T2 = T1 * n;
data_n[0] = T1;
data_n[1] = T2;
double Rob_1;
Rob_1 = 1 / (sqrt((1 + pow(x_1, 2)) * (1 + pow(x_1 * n, 2))));
double Kob;
Kob = Rrs_k / (Rr_1 * Rob_1);
tau = beta * T1;
data_n[2] = Kob;
data_n[3] = tau;
double Gs_op = -90, Rs_op = 1.55, B3_op = 2.7, omega_Ti = (2 * Pi) / B3_op, alpha = 0;
double Fs_op;
Fs_op = Gs_op * (Pi / 180);
complex<double> Ws_op;
Ws_op = Rs_op * exp(i * Fs_op);
double Rs;
Rs = abs(Ws_op);
const complex<double> one(1, 0);
complex<double> Wrs_op = Ws_op / (one - Ws_op);
double Rrs_op = abs(Wrs_op);
double Frs_op = arg(Wrs_op);
double omega_Td = omega_Ti * alpha;
double omega_Tf = 0.125 * omega_Td;
complex<double> Wf_op;
Wf_op = one / pow((i * omega_Tf + one), 2);
double Rf_op = abs(Wf_op);
complex<double> Wr_op;
Wr_op = one + (one / (i * omega_Ti)) + (i * omega_Td) * Wf_op;
double Rr_op = abs(Wr_op);
double Fr_op = arg(Wr_op);
Fob_op = Frs_op - Fr_op;
double x_0;
left = 0;
right = 2;
center = left;
while (right - left > eps)
{
if ((beta * center + atan(center) + atan(center * n) + Fob_op) * (beta * right + atan(right) + atan(right * n) + Fob_op) < 0)
{
left = center;
}
else
{
right = center;
}
center = (left + right) / 2;
}
x_0 = center;
double Rob_op;
Rob_op = 1 / sqrt((1 + pow(x_0, 2)) * (1 + pow(x_0 * n, 2)));
double Ks_op = Rrs_op / (Rr_op * Rob_op);
double Kr_n = Ks_op / Kob;
double Tres_n = data[3] * (x_1 / x_0);
double kp_n, Ti_n;
Ti_n = Tres_n / B3_op;
double Td_n = Ti_n * alpha;
kp_n = data[1] * Kr_n;
data_n[4] = kp_n;
data_n[5] = Ti_n;
//Порядок выходных данных - T1, T2, Kob, tau, kp_n, Ti_n
return data_n;
}