内弹道程序报告姓名：学号：专业：武器发射工程学院：能源与动力工程学院指导教师：1
程序功能对于一般的内弹道方程组，微分方程是非线性的，编制通过数值解法进行求解的程序
设计思路和整体方案（1）输入已知的数据，包括火炮构造，弹丸诸元，装药条件，起始条件，步长等参数
（2）进行一些常量的计算，方便在程序中进行使用
（3）计算内弹道方程中变量的初值
（4）将内弹道方程组进行量纲化，然后使用龙格库塔法计算微分方程，包括最大压力搜索和特征点判断，得出各个过程的数据
（5）将程序输出成表格和曲线
程序遇到的主要问题和解决方案（1）在程序中首先要编制一个四阶的龙格库塔子函数，同时还要编制一个微分方程右端子函数
在龙格库塔子函数编制过程中借鉴了之前导弹运动弹道程序的编制
voidrk（n,h）intn;doubleh;{externvoiddery（）;doublea[4],old_y[6],Y1[6],*dy;inti,j;dy=calloc(n,sizeof(double));a[0]=a[1]=h/2;a[2]=a[3]=h;dery(n,dy,Y);for(i=0;i
