Ошибка в коде?

@Ri · Регистрация: 15.09.2011

Author24 — интернет-сервис помощи студентам

Всем доброго времени суток! Буду очень благодарна,если кто-нибудь поможет разобраться с моей программой (потугами на нее)))) Суть дела: программа должна вычислить D. Задаются параметры. Скрины формул и данных прикладываю. Код вот он:

C++

#include <iostream>
#include <cstdlib>
#include <math.h>
#include <stdio.h>
using namespace std;
 
int main()
{
double b0=0;
double Ym=0;
double f=0;
double hm=0;
double D=0;
printf("\nPolutolshina sloya Ym = "); scanf("%lf", &Ym);
printf("\nRabochaa chastota f = "); scanf("%lf", &f);
printf("\nVisota sloya hm = "); scanf("%lf", &hm);
printf("\nUgol vihoda v gradusax b0 = "); scanf("%lf", &b0);
double temp=80.6*pow(10,12);
double f0=sqrt(temp);
double F=f/f0; //отношение рабочей частоты к критической
int re=6700;
double Rm=hm+re; //значение r, при котором концентрация максимальна
double r0=Rm-Ym; //начало слоя
double g=acos((re/r0)*cos(b0));
double A=1-(1/pow(F,2))+pow(r0/(F*Ym),2);
double B=-((2*Rm*pow(r0,2))/(pow(F,2)*pow(Ym,2)));
double C=pow((Rm*r0)/(F*Ym),2)-pow(re,2)*pow(cos(b0),2);
D=2*re*((g-b0)-((re*cos(b0)/2*sqrt(C)*log((B*B-4*A*C)/(4*C*((sin(g)+(sqrt(C)/r0)+(B/(2*sqrt(C))))*(sin(g)+(sqrt(C)/r0)+(B/(2*sqrt(C)))))))))); //расстояние скачка вдоль поверхности
printf("D = %lf \n",D);
return 0;
}

Заранее спасибо за то,что обратили внимание на мою тему))Не сердитесь,я просто ооочень давно не занималась написанием подобного( Должно получиться то,что в табличке в столбце D. При соответствующих данных. Но результат либо огромное число,либо NaN.

@~~IrineK~~ · 17.09.2015, 22:31

В строке 17 угол в градусах.
В строке 24 от него берется косинус. Там должны быть радианы.

Ну и дальше нужно смотреть.
Вы запустите в дебаге и отслеживайте, что выходит.

@Ri · 17.09.2015, 22:44 **[ТС]**

Спасибо,сейчас попробую. Про радианы сначала сделала,но потом сама же удалила. То есть в 24й строчке под косинусом в радианы переводить?Верно?

C++

1	double g=acos((re/r0)cos(b03.14/180));

Добавлено через 6 минут
Во всех местах где этот в градусах под косинусом видимо в радианы.

@~~IrineK~~ · 17.09.2015, 22:55

Все тригонометрические функции берутся от радиан.

@Ri · 17.09.2015, 23:02 **[ТС]**

Да,спасибо) Я то исправила это. Но,видимо,ошибка главная в другом) Ищу)

@~~IrineK~~ · 17.09.2015, 23:55

Ищите:
- деление на 0
- корень из отрицательного числа
- возведение отрицательного числа в степень и другие проблемы с
pow:

If the base is finite negative and the exponent is finite but not an integer value, it causes a domain error.
If both base and exponent are zero, it may also cause a domain error on certain implementations.
If base is zero and exponent is negative, it may cause a domain error or a pole error (or none, depending on the library implementation).
The function may also cause a range error if the result is too great or too small to be represented by a value of the return type

@zss · 18.09.2015, 08:29

Сообщение от Ri

double temp=80.6*pow(10,12);

Вы не знаете, как записываются константы типа double?

C++

1	double temp=80.6e12;

Зачем #include <iostream> если программа написана на Си?

Сообщение от Ri

pow(F,2)

Это принято писать как F*F

формулу в строке 28 разбейте на части:

C++

1
2
3

double D=B*B-4.0*A*C;
double t1=sin(g)+sqrt(C)/r0+B/(2.*sqrt(C));
....

@~~IrineK~~ Заблокирован
	17.09.2015, 22:31	2
	В строке 17 угол в градусах. В строке 24 от него берется косинус. Там должны быть радианы. Ну и дальше нужно смотреть. Вы запустите в дебаге и отслеживайте, что выходит. 0

@~~IrineK~~ Заблокирован
	17.09.2015, 22:55	4
	Все тригонометрические функции берутся от радиан. 0

@Ri 48 / 11 / 1 Регистрация: 15.09.2011 Сообщений: 87
	17.09.2015, 23:02 [ТС]	5
	Да,спасибо) Я то исправила это. Но,видимо,ошибка главная в другом) Ищу) 0

@~~IrineK~~ Заблокирован
	17.09.2015, 23:55	6
	Ищите: - деление на 0 - корень из отрицательного числа - возведение отрицательного числа в степень и другие проблемы с pow: If the base is finite negative and the exponent is finite but not an integer value, it causes a domain error. If both base and exponent are zero, it may also cause a domain error on certain implementations. If base is zero and exponent is negative, it may cause a domain error or a pole error (or none, depending on the library implementation). The function may also cause a range error if the result is too great or too small to be represented by a value of the return type 0

	18.09.2015, 08:29