void test_coth1()
{
	double Rohm=1, R1=1, alpha1=1, t1=1, f=1; 
	double yRe, yIm;
	complex cc1 = 1i;
	complex cc2 = coth(cc1);
	complex cc = Rohm+R1*coth((t1*1i*2*pi*f)^(alpha1/2))/(t1*1i*2*pi*f)^(alpha1/2);
	yRe = cc.m_re;
	yIm = cc.m_im;
}