Проблема с выводом данных в С++

[justso]

Пишу программу для сбора данных для дальнейшего построения графиков.

Вот вся программа (обратите внимание на выделение):

#include <iostream>
#include <cmath>
#include <fstream>
using namespace std;
#define PI 3.14159265

int main (void)
{
    int i=0, j=0;
    long double Tg, Ti, T, Tsafe, density, c, r, k, mu, mus, V, vis, time_step, V2, vel_inc, time=0, h, rate_of_change, change_in_temperature;
    double RED, Pr, muf;
    ofstream outfile1;  
    ofstream outfile2;  
    ofstream outfile3;
    
    outfile1.open("Project1VelH.out");
    outfile2.open("Project1TimeTemp.out"); 
    outfile3.open("Project1TotalTime.out");

    cout<< "This program calculates the parameters of cooling a copper sphere with the gas stream of given properties and velocity." <<endl;
    cout<< "\nPlease enter initial temperature of the sphere(C): --> ";
    cin>>Ti;
    cout<< "\nPlease enter gas temperature(C): --> ";
    cin>>Tg;
    cout<< "\nPlease enter safe handling temperature(C): --> ";
    cin>>Tsafe;
    cout<< "\nPlease enter sphere density (kg/m^3): --> ";
    cin>>density;
    cout<< "\nPlease enter sphere specific heat(J/kgK): --> ";
    cin>>c;
    cout<< "\nPlease enter sphere radius(m): --> ";
    cin>>r;
    cout<< "\nPlease enter desired time step for temperature history(s): --> ";
    cin>>time_step;
    cout<< "\nPlease enter gas properties:\nThermal Conductivity k (w/mK): -->";
    cin>>k;
    cout<< "\nPrandt Number for gas stream: --> ";
    cin>>Pr;
    cout<< "\nViscosity of Gas Stream at Tbulk (Ns/m^2): --> ";
    cin>>mu;
    cout<< "\nViscosity of Gas Stream at Tsurface (Ns/m^2): --> ";
    cin>>mus;
    cout<< "\nKinematic viscosity of gas stream (m^2/s): --> ";
    cin>>vis;
    cout<< "\nPlease enter gas stream velocities to test (initial, final, incriment (m/s)): --> ";
    cin>>V >> V2 >> vel_inc;
    
  
   
   double As = 4*PI*pow(r,2);
    double m = (density*4*PI*pow(r,3))/3;
    
     outfile1 << "\nVelocity        H\n\n";
     outfile3 << "Time needed to cool    Velocity\n\n";
   while (V<=V2)
   {T = Ti;

         RED = (V*2*r)/vis;
   muf = mu/mus;
   h = (k/(2*r))*(2+(0.4*pow(RED,0.5)+0.06*pow(RED,0.667))*pow(Pr,0.4)*pow(muf,0.25));
   cout << "H = " <<h;
  
   outfile1 << V << "        " << h <<endl;
   
   outfile2 << "Testing for Velocity " <<V <<"\n\n";
   outfile2 << "Time(s)       Temp(C)\n";
   while (T>=Tsafe+0.099)
   {
        cout << time << "     " << T << endl;
         outfile2 << time << "     " << T << endl;
   rate_of_change=(h*As*(Tg-T))/(m*c);
   change_in_temperature = time_step*rate_of_change;
   T = T + change_in_temperature;
    time = time + time_step;}
    outfile2 <<"\n\n\n\n\n\n\n";
   
     outfile3 << time <<"                " << V << endl;
    V=V+vel_inc;
    i=i+1;
    time=0;
    }
        system("pause");
    
    return 0;
}

Когда строка double m = (density*4*PI*pow(r,3))/3; находится в подобном состоянии, то я получаю один результат финального времени.

Если же я заменяю pow(r,3) на r*r*r, то результат иной. r = вводимая константа.
Объясните, что правильно и почему так. Уже 3 часа бьюсь.

Заранее благодарна.

[Stilet]

Во-первых: Нечего кросспостить. Сделала тему - все. Нужно будет модераторы перенесут.
Во-вторых:
Проверил так:

#include <cmath>

int _tmain(int argc, _TCHAR* argv[])
{
	double density=2,r=3,PI=3.1415926;
	double m = (density*4*PI*pow(r,3))/3;
	double q = (density*4*PI*(r*r*r))/3;
    
	return 0;
}

Результ одинаков полностью. так что что-то у тебя в другом месте не так.

Какие данные ты вводишь?

[justso]

r всегда одинакова и равняется 0.005. Остальное все меняется в зависимости от материала. Я проверяла, ошибка только в массе m, так как коеффициент теплопроводности h высчитывает правильно. Неверно лишь где-то тут:

rate_of_change=(h*As*(Tg-T))/(m*c);
   change_in_temperature = time_step*rate_of_change;
   T = T + change_in_temperature;
    time = time + time_step;

Так как именно тут используется m, я предположила что ошибка в ней.

Может где-то есть деление на int или что то вроде этого... это моя единственная мысль.

[NiCola999]

у вас за счет PI, при делении на int всегда будет double