solve6try.f90

运载火箭用拉瓦尔喷管计算（此程序参数保密

subroutine solve6()
use module1
external FS
dimension X(8),Y(8),Aa(8,8),Bb(8),JS(8)
double precision X,Y,Aa,Bb,Tta,Hha
do jj=1,nz
   x(1)=rh2oe+0.001
   x(2)=rco2e-0.001
   x(3)=rh2e+0.001
   x(4)=rcoe-0.001
   x(5)=rc2h2e+0.001
   x(6)=rn2e-0.001
   x(7)=29.4
   x(8)=0.001
   eps=1.0E-10
   tta=0.1
   hha=0.1
   call DNETN(8,X,Y,EPS,FS,Tta,Hha,Aa,Bb,L,JS)
   rh2o(jj)=x(1)
   rco2(jj)=x(2)
   rh2(jj)=x(3)
   rco(jj)=x(4)
   rc2h2(jj)=x(5)
   rn2(jj)=x(6)
   Mgw(jj)=x(7)
   dmdt(jj)=x(8)
   dmco2(jj)=Mc/Mco2*ppa(jj)*Mgw(jj)/Mco2*rco2(jj)
   dmh2o(jj)=Mc/Mh2o*ppa(jj)*Mgw(jj)/Mh2o*rh2o(jj)
   dmh2(jj)=2*Mc/Mh2*ppb(jj)*Mgw(jj)/Mh2*rh2(jj)
end do
end subroutine

subroutine FS(X,Y,N)
use module1
dimension X(N),Y(N)
double precision X,Y
Y(1)=x(7)*(x(1)/mh2o+x(2)/Mco2+x(3)/Mh2+x(4)/Mco+x(5)/Mc2h2+x(6)/Mn2)-1
Y(2)=(1+x(8)/alf(jj)*cp)*x(1)-rh2oe+ppa(jj)*x(7)/Mh2o*x(1)/alf(jj)*cp
Y(3)=(1+x(8)/alf(jj)*cp)*x(2)-rco2e+ppa(jj)*x(7)/Mco2*x(2)/alf(jj)*cp
Y(4)=(1+x(8)/alf(jj)*cp)*x(3)-rh2e-ppa(jj)*Mh2*x(7)/Mh2o/Mh2o*x(1)/alf(jj)*cp
Y(5)=(1+x(8)/alf(jj)*cp)*x(4)-rcoe-Mco/Mh2o*ppa(jj)*x(7)/Mh2o*x(1)/alf(jj)*cp-2*Mco/Mco2*ppa(jj)*x(7)/Mco2*x(2)/alf(jj)*cp
Y(6)=(1+x(8)/alf(jj)*cp)*x(5)-rc2h2e
Y(7)=(1+x(8)/alf(jj)*cp)*x(6)-rn2e
Y(8)=x(8)-ppa(jj)*Mc*x(7)*(x(2)/Mco2/Mco2+x(1)/Mh2o/Mh2o)
return
end subroutine