function det=detesferico(d3,q2)
det=-d3^2*sin(q2);
end