При любых введённых значениях выводит все результаты в виде: -1.#IND

#pragma once
 
namespace Program
{
class Ulovitel
{
    public:
    //производительность по запыленному воздуху
    double Pr;
    //скорость воздуха в рабочей зоне аппарата 
    double Sk;
    //температура воздуха, поступающего на очистку
    double Te;
    //плотность частиц
    double Pl;
    //начальная запыленность воздуха
    double Za;
    //давление в аппарате
    double Da;
 
    Ulovitel(void);
    Ulovitel(double pr, double sk, double te, double pl, double za, double da);
 
    double getPr();
    void setPr(double pr);
 
    double getSk();
    void setSk(double sk); 
 
    double getTe();
    void setTe(double te);
 
    double getPl();
    void setPl(double pl);
 
    double getZa();
    void setZa(double za);
 
    double getDa();
    void setDa(double da);
};
class Raschet: public Ulovitel
{
    public: 
    //отношение сечений патрубков
    double OtnSech;
    //отношение скоростей газа 
    double OtnSkor;
 
    Raschet(void);
    Raschet(double dvp, double dsg, double dnp, double dvov, double vrz, double dsh, double pvvp, double pvnp, double dv, double otnsech, double otnskor);
    static const double W;
 
    double getOtnSech();
    void setOtnSech(double otnsech);
    
    double getOtnSkor();
    void setOtnSkor(double otnskor);
 
    double getDvp();
    double getDsg();
    double getDnp();
    double getDvov();
    double getVrz();
    double getDsh();
    double getPvvp();
    double getPvnp();
    double getDv();
 
    //коэффициенты
    double X1;
    double X2;
    double getX1();
    double getX2();
};
}

#include "VihrPyl.h"
#include <math.h>
 
namespace Program
{
Raschet::Raschet(void)
{
    //отношение сечений патрубков
    double OtnSech=0.0;
    //отношение скоростей газа 
    double OtnSkor=0.0;
};
const double Raschet::W=3600.0;
Raschet::Raschet(double dvp, double dsg, double dnp, double dvov, double vrz, double dsh, double pvvp, double pvnp, double dv, double otnsech, double otnskor)
{
    //отношение сечений патрубков
    OtnSech=otnsech;
    //отношение скоростей газа 
    OtnSkor=otnskor;
}
Ulovitel::Ulovitel(void)
{
    //производительность по запыленному воздуху
    double Pr=0.0;
    //скорость воздуха в рабочей зоне аппарата 
    double Sk=0.0;
    //температура воздуха, поступающего на очистку
    double Te=0.0;
    //плотность частиц
    double Pl=0.0;
    //начальная запыленность воздуха
    double Za=0.0;
    //давление в аппарате
    double Da=0.0;
};
Ulovitel::Ulovitel(double pr, double sk, double te, double pl, double za, double da)
{
    //производительность по запыленному воздуху
    Pr=pr;
    //скорость воздуха в рабочей зоне аппарата 
    Sk=sk;
    //температура воздуха, поступающего на очистку
    Te=te;
    //плотность частиц
    Pl=pl;
    //начальная запыленность воздуха
    Za=za;
    //давление в аппарате
    Da=da;
}
double Raschet::getOtnSech() {return OtnSech;}
void Raschet::setOtnSech(double otnsech)
{ 
    OtnSech=otnsech;
}
double Raschet::getOtnSkor() {return OtnSkor;}
void Raschet::setOtnSkor(double otnskor)
{ 
    OtnSkor=otnskor;
}
double Ulovitel::getPr() {return Pr;}
void Ulovitel::setPr(double pr)
{ 
    Pr=pr;
}
double Ulovitel::getSk() {return Sk;}
void Ulovitel::setSk(double sk)
{
    Sk=sk;
}
double Ulovitel::getTe() {return Te;}
void Ulovitel::setTe(double te)
{ 
    Te=te;
}
double Ulovitel::getPl() {return Pl;}
void Ulovitel::setPl(double pl)
{ 
    Pl=pl;
}
double Ulovitel::getZa() {return Za;}
void Ulovitel::setZa(double za)
{ 
    Za=za;
}
double Ulovitel::getDa() {return Da;}
void Ulovitel::setDa(double da)
{ 
    Da=da;
}
        //получаемые и рассчитываемые значения
        double Raschet::getX1() {return pow((OtnSkor/(1+OtnSech)),0.5);}
        double Raschet::getX2() {return pow(((1-OtnSech*OtnSkor)/(1+OtnSech)),0.5);}
        double Raschet::getDvp() {return pow(4*Pr/(3.14*W),0.5);}
        double Raschet::getDsg() {return pow(4*Pr/(W*3.14*getDvp()),0.5);}
        double Raschet::getDnp() {return getX1()*getDvp();}
        double Raschet::getDvov() {return getX2()*getDvp();}
        double Raschet::getVrz() {return 3.1*getDvp();}
        double Raschet::getDsh() {return 0.93*getDvp();}
        double Raschet::getPvvp() {return 3.14*pow(getDnp(),2)/4;}
        double Raschet::getPvnp() {return (3.14/4)*(pow(getDvp(),2)-pow(getDvov(),2));}
        double Raschet::getDv() {return 0.1*getDvp();}
    }