需要解以下方程，其中只有δ是未知数，其它都是已知的列向量。我只把方程写出来了，但不知道在matlab里具体怎么解，求教！多谢！


        yita(i)=(delta(i)/(1-delta(i)))*(1-N(i)/(D(i)*delta(i)*(1-delta(i))));

        N(i)=sum(b(i)*b(j)*(w(i)-w(j))^2/delta(i)*(1-delta(i)))+...
            sum((sita(j)*b(i)*b(j)/delta(j))*((w(i)/((delta(j)/y(j)-1)^0.5))+w(j)*(delta(j)/y(j)-1)^0.5)^2)=0;


        D(i)=(sum(b(j)/(delta(j)*(1-delta(j))))+sita(j)*b(j)/(delta(j)*(delta(j)/y(j)-1)))...
            *(sum(b(j)w(j)^2/(delta(j)*(1-delta(j))))+sita(j)*b(j)*w(j)^2*(delta(j)/y(j)-1)/(delta(j)))...
            -(b(j)*w(j)/(delta(j)*(1-delta(j)))-sita(j)*b(j)*w(j)/delta(j))^2=0;

         b(i)=x1(i)^2/(p(i)*q(i))=0;