一维激波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;
}
總結
- 上一篇: 数据结构小作业---家族族谱
- 下一篇: s3c2440移植MQTT