Восстановить работоспособность программы! для курсача

#ifndef CALCENGINE_H
#define CALCENGINE__H
 
#include <math.h>
#include "gplots.h"
#define MAX_POINTS_Y 1000 //. Максимальное число точек для Y(Q)
//---------------------------------------------
class CEngine
{
private:
double Y,
 
Z1, // Производительность домны с учётом запаздывания 
dzl, // Скорость изменения производительности домны с учётом запаздывания
 Z, // Производительность домны с учётом запаздывания и инерционности 
dz, // Скорость изменения производительности домны
ZMaxPred, // Максимальная производительность домны, запомненная на предыдущем шаге
DeltaZ, // Приращение производительности домны на текущем шаге 
DeltaZMaxPred,  //  Максимальная разность между приращениями
производительности домны, запомненная на предыдущем шаге
                            U, // Управляющая функция 
                           Sigma, // Переключающая функция 
                           Kim; // Скорость ИМ
   public:
double Q, Q0, // Текущий и начальный расход природного газа
Qmin, Qmax, // Мах и Min расход природного газа по статике deltaQ, // Изменение расхода природного газа
Zn, // Зона нечувствительности
TAUimp, // Время импульса
 TAUz // Время запаздывания
 Tstep, // Время оного шага системы
 Tim, // Время полного хода ИМ
ТоЬ; // Постоянная времени объекта 
int               Tost, // Время остановки ИМ
Nrepeat; // Число повторных запусков системы после нахождения экстремума 
double а[5], // Коэффициенты уравнения 
CoefAlpha, CoefBeta, Tdr;
FPoint *QPoints, *ZQPoints,
 
*QTPoints, *ZTPoints, *ZlTPoints,
*YQdrPoints;
int ZQiterations; // Число итераций для цикла пересчёта
int MaxZQiterations; // Максимально возможное число итераций 
 
 
//—Конструктор и деструктор —
                      CEngine();
~CEngine();
//- 
void Calculate();
void ShowCalcResult(double &kim, double &tauimp) {kim = Kim; tauimp = TAUimp;}; 
private:    
                      void CalcY(int hour); // Вычисление Y(Q)
};
//==============================================
CEngine::CEngine()
{
QPoints = new FPoint[MAX_POINTS_Y]; // Создание массива структур FPoint размером MAX_POINTS__Y
ZQPoints = new FPoint[l];
QTPoints = new FPoint[l];
ZTPoints = new FPoint[l];
                    ZITPoints = new FPoint[l];
YQdrPoints = new FPoint[MAX_POINTS_Y];
MaxZQiterations = 40; // Максимальное число итераций 
                              
                                                       }
 
 
// - - - - - - - - - - - - - - 
CEngine::~CEngine() 
 
 
                                                       { 
delete []QPoints; delete []ZQPoints; delete []QTPoints; delete []ZTPoints; delete []ZlTPoints; delete []YQdrPoints; 
 
 
                                                        }
 
// - - - - - - - - - - - -
void CEngine::CalcY(int hour) 
                                            
double Qtmp = Q+CoefAlpha*hour;
Y = a[0];
int i(l);
while (i<5) {
Y += Qtmp*a[i];
 
Qtmp *= (Q+CoefAlpha*hour); // Увеличение степени в уравнении
i++;
}
 
Y += CoefBeta*hour;
}
// - - - - - - - - - - - - - - - - - - -
void CEngine::Calculate()
{
ZQiterations = MaxZQiterations;
// Вывод первого графика статики Y(Q) =======================
double step = (Qmax - Qmin)/MAX_POINTS_Y;
int i(0);
for (double st = Qmin; st < Qmax; st += step, ++i) {
// Расчёт
Q = st;
CalcY(O); // Вызов функции void CEngine::CalcY()
QPoints[i].x=st; // Заполнение данных 1 графика (без дрейфа) QPoints[i].y = Y;
CalcY(Tdr);
YQdrPoints[i].x = st;
YQdrPoints[i].y = Y;
 
}
// Перевыделение памяти под графики
delete []ZQPoints;
ZQPoints = new FPoint[(Tstep+1 )*ZQiterations+1 ];
int qtSize = Tstep*ZQiterations+l;
delete []ZTPoints;
delete []ZlTPoints;
free(QTPoints);
QTPoints = 0;
QTPoints = new FPoint[qtSize+ZQiterations];
ZTPoints    =   new FPoint[qtSize];
ZITPoints = new FPoint[qtSize];
// Расчёт графиков Z1(t), Z(t), Q(t) = = = = = = = = = = = = = = = = = = = = = =
Q = Q0;
CalcY(0);
Z1 = Z = Y;
ZMaxPred = Z;
dzl = dz = DeltaZ = DeltaZMaxPred = 0;
Kim = (Qmax-Qmin)/Tim;
TAUimp = deltaQ/Kim;
Sigma = -1; // Задали произвольно направление движения ИМ на уменьшение Q
int Lines = 0;
int zpiterations = 0;
int qiterations = 0;
ZQPoints[zpiterations].x = Q;
ZQPoints[zpiterations].y = Z;
++zpiterations; // На следующую итерацию ZQPoints                                    QTPoints[qiterations].x = Lines-qtSize+1;
ZTPoints[Lines].x = ZlTPoints[Lines].x = Lines;
QTPoints[qiterations].y = Q;
ZTPoints[Lines].y = Z;
Z1TPoints[Lines].y = Z1;
++Lines; // На следующую итерацию
++qiterations;
int TostTemp = 0;// Переменная для расчёта времени остановки Tost
int ExtremumFindCount = 0; // Число повторов остановок ИМ
bool Pusk = false;
bool Stop = false;
double XI = 0, X2 = 0;
for (int j = 0; j < ZQiterations; ++j) { //
if (DeltaZ>=0) {
if ((DeltaZ-DeltaZMaxPred)+Zn>=0) U=1; else U=0;
//if (DeltaZ-fabs(DeltaZMaxPred)>=0) U = 1; else U=0;
} else {
if (DeltaZ+Zn<0) U = -1; //else U=1;
}
if (U!=l && Sigma)
 
if(!U) {
Sigma = 0;
TostTemp = Lines;
++ExtremumFindCount; // Число остановок ИМ
} else
Sigma = -Sigma;
}
if (= = -1) {   //  Когда  происходит    реверс ИМ при    смене  знака Sigma
ZMaxPred = Z;
 
DeltaZMaxPred = DeltaZ;
}else {
if (Z>ZMaxPred) ZMaxPred = Z; if(DeltaZ>DeltaZMaxPred) DeltaZMaxPred = DeltaZ;
}
if (Pusk) {
if (!Sigma) Sigma = -1; else Sigma = -Sigma;
 
Pusk = false;
ZMaxPred = Z;
DeltaZMaxPred = DeltaZ;
}   
// Устранение колебательного режима в районе максимума 
if (U== -l && Sigma==l) XI =Q;
if (U==-l && Sigma==-1) X2=Q;
if (Q == (Xl+X2)/2 && !Stop)
 {Sigma = 0;
Stop = true;
TostTemp = Lines;
++ExtremumFindCount; // Число остановок ИМ
}
if (!Pusk && !Stop) { 
if (Sigma) { //если Sigma<>0
Sigma = (Q == Qmax) ? -1 : ((Q == Qmin) ? 1 : Sigma);
Q += Sigma*Kim*TAUimp;
Q = (Q > Qmax) ? Qmax : ((Q<Qmin)?Qmin:Q); // Условие концевых
выключателей
}
}
ZQPoints[zpiterations].x = Q;
 
ZQPoints[zpiterations].y = Z;
QTPoints[qiterations].x = Lines-qtSize;
QTPoints[qiterations].y = Q;
++qiterations;
++zpiterations; // На следующую итерацию ZQPoints
int i = 0;
for (; i < Tstep; ++i) { // Цикл итераций //
if (Lines<=Tdr) CalcY(Lines); else CalcY(Tdr);
Z1 += dzl;
dzl = (Y-Zl)/TAUz;
Z+=dz;
dz = (Zl-Z)/Tob;
QTPoints[qiterations].x = Lines-qtSize+1;
ZTPoints[Lines].x = ZlTPoints[Lines].x = Lines;
QTPoints[qiterations].y = Q;
ZTPoints[Lines].y = Z;
ZlTPoints[Lines].y = Z1;
++Lines; // На следующую итерацию ZQPoints[zpiterations].x = Q;
ZQPoints[zpiterations].y = Z;
++zpiterations; // На следующую'итерацию ZQPoints ++qiterations;
if (TostTemp>0)
if (Lines-TostTemp==Tost) { // Проверка времени стояния: простояла система Tost или нет
Pusk = true;
Stop = false;
XI = 0;
X2-0;
if (ExtremumFindCount==Nrepeat+l) {
ZQiterations = j+1; // Остановка после стояния системы
}
break;
}
}
DeltaZ = Z-ZMaxPred;
}
ZPiterations= zpiterations;
}
//-------------------
#endif//_CALCENGINE_H