#include“class.h“
#include
#include

constexpr?auto?gamma?=?1.4;//氣體比熱比
constexpr?auto?R?=?287.0;//空氣常數
constexpr?auto?cfl?=?0.5;//庫朗數
constexpr?auto?timeStep?=?1000;//時間步
constexpr?auto?deltaX?=?0.1;//x方向間距

using?namespace?std;


int?main（）?{
	cout?<	cellField?rho?V?T?Aa;//定義物理量（rho-密度，V-速度，T-溫度，A-截面積）
	cellField?drhodt?dVdt?dTdtrhopVpTp;//定義預估步偏導預估值
	cellField?drhopdt?dVpdt?dTpdtdrhodtavdVdtavdTdtav;//定義校正步偏導，校正值
	cellField?X?deltaTime;//定義空間步、時間步
	double?i?jk;//定義下標
	ofstream?Rfile;//定義文件類
	Rfile.open（“data.dat“）;//打開文件

	//*********************初始化**************************//
	for?（i?=?0;?i?<=?cellNumber;?i++）
	{
		X[i]?=?i?*?deltaX;
		A[i]?=?1.0?+?2.2?*?（X[i]?-?1.5）*?（X[i]?-?1.5）;//公式7-73，p219
		rho[i]?=?1.0?-?0.314?*?X[i];//公式7-74a，p220
		T[i]?=?1.0?-?0.2314?*?X[i];//公式7-74b，p220
		V[i]?=?（0.1?+?1.09?*?X[i]）?*?sqrt（T[i]）;//公式7-74c，p220
		a[i]?=?sqrt（T[i]）;//公式7-75，p222
	}
	//輸出初始場
	Rfile?<	Rfile.precision（3）;
	Rfile.setf（ios::fixed）;
	Rfile?<	for?（i?=?0;?i?<=?cellNumber;?i++）
	{
		Rfile?<			<	}
	//******************************************************//

	for?（k?=?0;?k?<=?timeStep;?k++）//開始計算timeStep個時間步
	{	
		//*******************求最小時間步*****************//
		cellField?Va?=?V?+?a;
		deltaTime.inverse（Va?cfl?*?deltaX）;//公式7-67，p217
		j?=?0;
		for?（i?=?0;?i?<=?cellNumber;?i++）
		{
			if?（deltaTime[j]?>?deltaTime[i]）
				j?=?i;//求最小時間間隔，將下標儲存為j，即deltaTime[j]為最小時間步
			//cout?<		}
		//*************************************************//

		//********************預估步***********************//
		for?（i?=?0;?i?		{
			drhodt[i]?=?-V[i]?*?（rho[i?+?1]?-?rho[i]）?/?deltaX
				-?rho[i]?*?（V[i?+?1]?-?V[i]）?/?deltaX
				-?rho[i]?*?V[i]?*?（log（A[i?+?1]）?-?log（A[i]））?/?deltaX;//7-51p215
			dVdt[i]?=?-V[i]?*?（V[i?+?1]?-?V[i]）?/?deltaX
				-?1?/?gamma?*?（（T[i?+?1]?-?T[i]）?/?deltaX?+?T[i]?/?rho[i]?*?（rho[i?+?1]?-?rho[i]）?/?deltaX）