6 phi11(0.0), phi12(0.0), phi21(0.0), phi22(0.0), dphi11(0.0), dphi12(0.0),
7 dphi21(0.0), dphi22(0.0), frequency(0.0)
12 const double p12,
const double p21,
const double p22,
13 const double dp11,
const double dp12,
const double dp21,
15 phi11(p11), phi12(p12), phi21(p21), phi22(p22), dphi11(dp11),
16 dphi12(dp12), dphi21(dp21), dphi22(dp22), frequency(f)
void f(vector< double > &v1, vector< double > &v2, vector< double > &v3, vector< double > &v4)