石墨烯电导率，介电常数

#include <complex>
using namespace std;

typedef complex<double> Complex;

void compute()
{
	double T = 300;
	double uc = 0.8;
	double tor = 0.43e-3;
	double delta = 1e-9;

	double cs = 3e8;
	double ev = 1.602e-19;
	double Kb = 1.38e-23;
	double h = 6.626e-34;
	double eps0 = 8.85e-12;
	Complex j0(0,1);
	double pi = 3.1415926;
	double h1 = h/2/pi;

	double p = tor*ev/h*2*pi;

	double freq = 30e12;
	double w = 2*pi*freq;
	double k0 = w/cs;
	double lumda0 = cs/freq;

	double UC = uc*ev;


	//calculate interband:
	Complex sigmaInter = -j0*ev*ev/(4.0*pi*h1) * log( (2.0*abs(UC)-(w-j0*2.0*p)*h1)/(2.0*abs(UC)+(w-j0*2.0*p)*h1) );

	//calculate intraband:
	Complex sigmaIntra = -j0*ev*ev*Kb*T / (pi*(h1*h1)*(w-j0*2.0*p)) * (UC/Kb/T+2.0*log(exp(-UC/Kb/T)+1));
	
	Complex sigma = sigmaInter + sigmaIntra;

	Complex eps = -j0*sigma/delta/w/eps0;

	double tan_delta = eps.imag()/eps.real();

}