phi=normal(0,1)
phi1=normal(0,1)
phi2=normal(0,1)
x1=d1*g1(x)*g1p(x)*phi1^2+d2*g2