几年前，当自己还是个“热血”青年的时候

那时候对什么事情都看不惯

都是带着批判的眼光看待世界的一切

现在，当我看到那些小屁孩在否定世界批判世界的时候

又站在了“保守”的一边。

正如英语书上说的，60年代的年轻人潮压倒了成年人。现在许多成年人试图恢复稳定的局势。

主题：斜激波中的马赫数、偏转角、激波角程序！

作者：生来孤独（rockandtcl）

语言：Fortran

程序：Compaq VF6.0，Tecplot360

日期：2008-07-10

备注：版权所有，如果引用请注明出处！

上次发的时候没有很好的去修改，其实在最小角度方面实际上是有理论依据的。

假设偏转角很小。那么可以把激波角认为是马赫角！

于是，代码就变得很简单，同时也是很精确的！

 PARAMETER(NI=31,NJ=201)
 PARAMETER(PI=3.1415926535897)
 PARAMETER(GM=1.4)
 DIMENSION AM(NI)
 OPEN(1,FILE='OS.PLT')
 WRITE(1,*)'VARIABLES=<greek>d</greek>,<greek>s</greek>'
 
 DO I=1 , 10
  AM(I)=1.0+0.05*I
 ENDDO

 DO I=11 , 15
  AM(I)=AM(10)+0.1*(I-10)
 ENDDO

 DO I=16 , 25
  AM(I)=AM(15)+0.2*(I-15)
 ENDDO

 DO I=26 , 28
  AM(I)=AM(25)+0.5*(I-25)

主题：斜激波中的马赫数、偏转角、激波角程序！

作者：生来孤独（rockandtcl）

语言：Fortran

程序：Compaq VF6.0，Tecplot360

日期：2008-06-24

备注：版权所有，如果引用请注明出处！

今天无聊，上空动2的时候看到了这张表。于是就决定写一个程序来生成这个表！

代码如下：

 PARAMETER(NI=31,NJ=201)
 PARAMETER(PI=3.1415926535897)
 PARAMETER(GM=1.4)
 DIMENSION AM(NI)
 OPEN(1,FILE='OS.PLT')
 WRITE(1,*)'VARIABLES=d,s'
 
 DO I=1 , 10
  AM(I)=1.0+0.05*I
 ENDDO

 DO I=11 , 15
  AM(I)=AM(10)+0.1*(I-10)
 ENDDO

 DO I=16 , 25
  AM(I)=AM(15)+0.2*(I-15)
 ENDDO

 DO I=26 , 28
  AM(I)=AM(25)+0.5*(I-25)
 ENDDO

 AM(29)=8.0
 AM(30)=10.0
 AM(31)=20.0

 DO I=1,NI
  DB=PI/2