文档介绍：附录程序
#define _CRT_SECURE_NO_DEPRECATE
#include <>
#include <iostream>
#include <>
#include <>
#include <>
#include <iomanip>//C head file
#include <>
#include <>
#include <string>
#include <fstream>
using namespace std;
/*初始参数*/
//////////燃烧室内参数
//***燃烧的参数
#define U1 //燃速系
#define an //燃速指数
#define omga //装药质量
#define am //弹丸质量
#define f 980000 //***力
#define F1 //阻力系数
#define STA //（k-1）
#define DP 1600 //***装填密度
#define ALF -3 //***燃气余容
double PG=30e6 ; //挤进压力
//***参数（单基9/7）
#define HE1 -3 //第二种***半弧厚
#define HD0 -3 //第二种***内孔直径
#define HD2C 12e-2 //第二种******长
//***形状特征量Π1、Q1、β
double XHT,XHQ,XHB;
//double XHB;
//碎粒燃烧结束时的Zb，***形状参数χ、λ和μ
double Zb0,CI,ALM,AMU;
//碎粒燃烧阶段的形状参数χs、λs
double CIS,ALMS;
//////////火炮参数
#define V0 -3
#define omga2 //第一种药质量
#define f2 980000 //第一种药***力
#define DP2 1600 //第一种药密度
#define ALF2 1e-3 //第一种药余容
//***参数（选用杆状药）
#define HE12 -3 //***半弧厚
#define HD02 -3 //内孔直径
#define HD2C2 -3 //***长
//#define f2 980000
double Z2,PSI2,dZ2,Zb2,Zm2;
double DLT2;
//***形状参数χ、λ和μ
double CI2,ALM2,AMU2;
double XHB2;
#define pi //圆周率
#define atm //大气压
#define h 1e-7