Numerical coupled integration
fn(q)= int K1(q,q’) * q’² * fn(q’) dq’ + K2(q,q’)*q’² fd(q’)
fd(q)= int k3(q,q’) * q’² * fn(q’) dq’
fn(q)= int K1(q,q’) * q’² * fn(q’) dq’ + K2(q,q’)*q’² fd(q’)
fd(q)= int k3(q,q’) * q’² * fn(q’) dq’