solve.f90

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

subroutine solve()
use module1
real :: mtra(nr),mtrb(nr),mtrc(nr),mtrd(nr),mtrp(nr),mtrq(nr),cm(nr)

do i=1,nr
   mtra(i)=ap(i,kk)
   mtrb(i)=as(i,kk)
   mtrc(i)=an(i,kk)
   mtrd(i)=b(i,kk)
end do
!-----------------------------
mtrp(1)=mtrb(1)/mtra(1)
mtrq(1)=mtrd(1)/mtra(1)
do i=2,nr
   mtrp(i)=mtrb(i)/(mtra(i)-mtrc(i)*mtrp(i-1))
   mtrq(i)=(mtrd(i)+mtrc(i)*mtrq(i-1))/(mtra(i)-mtrc(i)*mtrp(i-1))
end do
!-------------------------------
mtrtd(nr,kk)=mtrq(nr)
do i=nr-1,1,-1
   mtrtd(i,kk)=mtrp(i)*mtrtd(i+1,kk)+mtrq(i)
end do
do i=1,nr
   if(mtrtd(i,kk)<=tini)then
     mtrtd(i,kk)=tini
   end if
end do

end subroutine