function R = bayes_risk_2normal(distribution1, distribution2, q)

Am = distribution1.Mean;
As = distribution1.Sigma;
Ap = distribution1.Prior;
Cm = distribution2.Mean;
Cs = distribution2.Sigma;
Cp = distribution2.Prior;

riskA = normcdf(q.t2,Am,As) - normcdf(q.t1,Am,As);
riskC = normcdf(q.t1,Cm,Cs) + 1 - normcdf(q.t2,Cm,Cs);
R = Ap*riskA + Cp*riskC;