(fitteia's extended models library)
double R1CRsp1(double f, double Cnh, double tq, double fq1, double fq2, double Theta, double Phi)
double R1CRsp1(double f, double Cnh, double tq, double fq1, double fq2, double Theta, double Phi)
{
double af,c1,c2;
c1 = 1.0/3.0+pow(sin(Theta)*cos(Phi),2.0);
c2 = 1.0/3.0+pow(cos(Phi),2.0);
af= Cnh*(
c1*( tq/(1+pow(2*pi*(f-fq1)*tq,2.0)) + tq/(1+pow(2*pi*(f+fq1)*tq,2.0)) )
+ c1*( tq/(1+pow(2*pi*(f-fq2)*tq,2.0)) + tq/(1+pow(2*pi*(f+fq2)*tq,2.0)) )
+ c2*( tq/(1+pow(2*pi*(f-(fq1-fq2))*tq,2.0)) + tq/(1+pow(2*pi*(f+(fq1-fq2))*tq,2.0)) )
);
return af;
}