Технический форум - Найти ошибку в математических формулах. Pascal

Цитата:

Сообщение от molo4nik (Сообщение 1031168)

1)if x>0 then
p:=(a/(1-x)) + (b+ sqr(sin(x*x*x)/cos(x*x*x))/ (a*(exp(a*x) - exp(-x)) ));

if x=0 then
begin
z:=exp(ln(b+a*x) / 4);
P:= b*b + exp(ln(cos(18.6/z)*3) ) ;
end;

if x<0 then
P:=sqr((sin((x*x*x)/2*a))) - ln( ( x+b*b*b*b))/ sqrt(abs(x-1));

2)P:=(sqr(cos(x*x)/sin(x*x))/ 2*a + b);
P:=P+ exp(sqr(sin((a*x)*(a*x))));

S:=0;
h:=5*a;
n:=111-5*a;
for i:=h to n do
S:=(ln(abs(exp(ln(-1)*h)*(exp(ln(b)*5) +a)))) / 0.2*h;

P:=P*S;

Извольте.
Во-первых, давайте уберем эту чертову кучу ненужных пробелов. Они, в принципе, не мешают и к ошибкам не приводят, но зачем?

1а (случай х>0).
Неверная расстановка скобок. При такой записи во втором слагаемом на знаменатель делится не весь числитель, а только квадрат тангенса. Правильно:
if x>0 then
p:=(a/(1-x))+(b+sqr(sin(x*x*x)/cos(x*x*x)))/(a*(exp(a*x)-exp(-x)));

1b (случай х=0).
При вычислении z пропущен модуль. При вычислении Р - неверная расстановка скобок. Правильно так:
if x=0 then
begin
z:=exp(ln(Abs(b+a*x))/4);
P:= b*b+exp(ln(cos(18.6/z))*3);
end;

1c (случай х<0).
В обоих слагаемых - неверная расстановка скобок. Правильно так:
if x<0 then
P:=sqr(sin(x*x*x/(2*a)))-ln((x+b*b*b*b)/sqrt(abs(x-1)));

2. Прежде всего - запомните раз и навсегда: в области вещественных чисел функция логарифм (Ln) определена только для строго положительных аргументов. Запись типа Ln(-1) - бессмысленна и приведет к ошибке.
Как из этого выходить?
В некоторых версиях Паскаля есть функция возведения в степень, но в некоторых таковой функции нет. Поэтому найдем "универсальное" решение (попутно исправив, как и везде, неверную расстановку скобок и другие ошибки):
P:=sqr(cos(x*x)/sin(x*x))/(2*a+b);
M:=exp(sqr(sin((a*x)*(a*x))));

S:=0;
n1:=5*a;
n2:=111-5*a;
for h:=n1 to n2 do
begin
if (h mod 2)=0 then Q:=Ln(Abs(exp(ln(b)*5)+a)) else Q:=Ln(Abs(-exp(ln(b)*5)+a));
S:=S+Q/(0.2*h);
end;
P:=P+M*S;