Could you guys do me a favor for solving a equation set?

Z=d/4*(1-SIN(X)/X)

X=8q/(D^2*Y)+SIN(X)

Y=1/n*Z^(2/3)*i^(1/2)

In this equation set, X,Y&Z are the unkown parameters, the others say,

d, q, n&i are known. SO in python, how to program it to represent X, Y

and Z in the form of d, q, n and i?

Thanks very much.