# Propagate error for n by quadrature# Declare function variablesd=var('d')m=var('m')R=var('R')sigd=var('sigd')sigt=var('sigt')# Declare function for ntop=4*R*d*sin(t)bottom=R*d*sin(t)-4*mn=sqrt(top/bottom)# Take partial derivatives of ndt=diff(n,t)dd=diff(n,d)