一维激波C++编程

#include<iostream>
using namespace std;
#define Gamma  1.4;//比热比 
double R = 287.0;//理想气体常数J/kg/K
const int Nnode = 101;
double* CalDem(double* temp, double* pre)
{
 double *dem = new double [Nnode]{ 1 };
 for(int i=0; i < Nnode;i++)
 {
	 dem[i] = pre[i]/(R * temp[i]); 
 }
 return dem;
}



int main()
{
	double Length = 10;
	int dx = Length / (Nnode-1);
	double Vel[Nnode] = { 0 };

	double Temp[Nnode] = {300 };
	fill_n(Temp, 101, 300.1);
	double Pre[Nnode] = { 101250 };
	fill_n(Pre, 50, 101250);
	fill_n(Pre+50, 51, 1012500);
	double Dem[Nnode] = { 1 };
	memcpy(Dem, CalDem(Temp, Pre), Nnode * sizeof(double));
	


	for(int i =0;i<Nnode;)
	{
		cout << Dem[i] << endl;
		i += 1;
	}


	return 0;
}