I need the correct coding to solve the following.

theta=(2*i0*q0)/(g*n*D*q)

q=sqrt(((2*n*D)*(Av+(n*F)+R(L)))/(((n-2)*(1-(D/P))*hv)+hb+hv+(g*n^2*D*theta)+((p*I)/(1+I*t))))

TOT=((D*(Av+(n*F)+R(L)))/q)+(q/(2*n))*((n-2)*(1-(D/P))*hv)+hb+hv+((p*I)/(1+I*t))+(hb*K*sigma*sqrt(L))+((g*n*theta*q*D)/2)+(i0*q0*(ln(theta0/theta)))

n*(n*-1)≤Av(hb+hv-(2*hv*(1-(D/P)))+((p*I)/(1+It)))/((Av*g*D*theta)+(hv*(1-(D/P))*(F+R(L)))

≤ n*(n*+1)

For the data

P=7500; hb=50; hv=40; D=6000; Av=200; I=0.15; F=100; t=0.05; p=3;= 6, 000 unit/year; sigma=0.7; K=2.33; g=15; and q0(θ) = 400 ln(0.000009/θ).

Also, for i=0, 1, 2, 3. Li=8,6,4,3. and R(Li)=[0 5.6 22.4 57.4];

To find Q, theta, TOT, n*