/////////////////////////////   计算发动机比冲  //////////////////////////////
#include <iostream>
#include <armadillo>
#include <time.h>
#include <stdio.h>
#include <stdlib.h>
#include "GPM.h"

using namespace std;
using namespace arma;

///////////////////////////////   计算发动机比冲  //////////////////////////////
int calp(vec Ma, vec phit,vec alpha, vec * CT, vec * IspDim)
{
	int n = size(Ma)(0);
	vec  Alpha = alpha / datum::pi * 180;
	*IspDim = 983.5 + 2187 * phit + 25.14*Ma - 35.11*phit%Ma - 562.5*phit%phit - 4.84*Ma%Ma - 6.782*phit%Ma%Ma + 31.55*phit%phit%Ma;;
	*CT = -0.3889 + 0.004263*Alpha + 0.3786*Ma + 0.005587*Alpha%Ma - 7.816e-9*Alpha%Alpha  - 0.008242*Ma%Ma - 2.86e-11*Alpha%Alpha%Ma;
	return 1;
}