solve6.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,Tt,Hh
 x(1)=0.5
 x(2)=0.1
 x(3)=0.1
 x(4)=0.1
 x(5)=0.1
 x(6)=0.1
 x(7)=29.4
 x(8)=0.001
 EPS=1.0E-10
 tt=0.1
 hh=0.1
 CALL DNETN(8,X,Y,EPS,FS,Tt,Hh,Aa,Bb,L,JS)
 rh2o=x(1)
 rco2=x(2)
 rh2=x(3)
 rco=x(4)
 rc2h2=x(5)
 rn2=x(6)
 Mgw=x(7)
 dmdt=x(8)
 dmco2=Mc/Mco2*ppa*Mgw/Mco2*rco2
 dmh2o=Mc/Mh2o*ppa*Mgw/Mh2o*rh2o
 dmh2=2*Mc/Mh2*ppb*Mgw/Mh2*rh2

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*cp)*x(1)-rh2oe+ppa*x(7)/Mh2o*x(1)/alf
 Y(3)=(1+x(8)/alf*cp)*x(2)-rco2e+ppa*x(7)/Mco2*x(2)/alf
 Y(4)=(1+x(8)/alf*cp)*x(3)-rh2e-ppa*Mh2*x(7)/Mh2o/Mh2o*x(1)/alf+ppa*x(7)/Mh2*x(3)/alf
 Y(5)=(1+x(8)/alf*cp)*x(4)-rcoe-Mco/Mh2o*ppa*x(7)/Mh2o*x(1)/alf-2*Mco/Mco2*ppa*x(7)/Mco2*x(2)/alf
 Y(6)=(1+x(8)/alf*cp)*x(5)-rc2h2e-Mc2h2/Mh2*x(7)/Mh2*ppa*x(3)/alf
 Y(7)=(1+x(8)/alf*cp)*x(6)-rn2e
 Y(8)=x(8)-ppa*Mc*x(7)*(x(2)/Mco2/Mco2+x(1)/Mh2o/Mh2o)-ppb*2*Mc*x(7)*x(3)/Mh2/Mh2
 return
end subroutine