Найти значение интеграла методом Симпсона с точностью 0,0001

@Sardorxon · Регистрация: 17.05.2016

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

считая а (алфа) = 1,2,3..20,определить при каких значениях параметра а (алфа) числовые значения интеграла будут min и max ?
у меня прога готова )))

есть начальное и конечное значение (a;b) если в место " a " вводишь 1,2,3..n-1,n то прога работает и вдоет ответь .. но когда в место "a" вводиш 0 (нол) ответа нету ((
с значением "b" нет проблем

Pascal
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
uses crt;
label 1,2;
var
a,b,eps,h,s,s1,alfa :real;
i,n,c1,c2 :longint;
rt :text;
function f(x,alfa:real):real;
begin
f:=(sqrt(exp(ln(x)*alfa))*ln(alfa*x+1)/sqrt(exp(ln(x)*2*alfa))+alfa*exp(ln(x)*alfa)+1);
end;
{---------}
begin clrscr;
assign(rt,'int.txt'); rewrite(rt);
writeln('vveite nachalnoe znachenie integrala a=');
read(a);
writeln('vveite konechnoe znachenie integrala b=');
read(b);
writeln('vveite tochnost eps=');
read(eps);
alfa:=1;
repeat
n:=20;
n:=2*n;
h:=abs(b-a)/n;
s:=f(a,alfa)+f(b,alfa);
c1:=1; c2:=-1;
for i:=2 to n-1 do
begin
c1:=c1*c2;
s:=s+(3-c1)*f(a+i*h,alfa);
end;
s:=h/3*s;
2: c1:=1; c2:=-1;
n:=2*n;
h:=abs(b-a)/n;
s1:=f(a,alfa)+f(b,alfa);
for i:=2 to n-1 do
begin
c1:=c1*c2;
s1:=s1+(3-c1)*f(a+i*h,alfa);
end;
s1:=h/3*s1;
if abs(s-s1)<=eps then
goto 1 else begin s:=s1; goto 2;end;
1:writeln('alfa= ',alfa:5:0,' znachenie integrala=',s1:10:2);
writeln(rt,'alfa= ',alfa:5:0,' znachenie integrala=',s1:10:2);
alfa:=alfa+1;
until alfa>10;
close(rt);
end.

оба они a,b: real но где то не то внутри программы.. где именно исправить, не как не смог найти неполадку .. помогите пожалуйста ((

@Sardorxon · 17.05.2016, 19:54 **[ТС]**

скриншот могу вложит

ФедосеевПавел · 17.05.2016, 20:06

Логарифм 0 - неопределён.
Измените формулы вычисления степени или попробуйте штатную для PABC.NET.

@Sardorxon · 17.05.2016, 20:12 **[ТС]**

еще по подробнее можете объяснить.. а то я не такой программист

Добавлено через 1 минуту
изменит эти формулы ?

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

exp(ln(x)*alfa

ФедосеевПавел · 17.05.2016, 20:13

exp(ln(x)*alfa) замените Power(x, alpha)
И в нескольких местах аналогично

@Sardorxon · 17.05.2016, 20:20 **[ТС]**

все равно не па могло прога работает но ответа нету (((

f:=(sqrt(Power(x,alfa))*ln(alfa*x+1)/sqrt(Power(x,2*alfa))+alfa*Power(x,alfa) +1);

Добавлено через 1 минуту
alfa= 1 znachenie integrala= NaN
alfa= 2 znachenie integrala= NaN
alfa= 3 znachenie integrala= NaN
alfa= 4 znachenie integrala= NaN
alfa= 5 znachenie integrala= NaN
alfa= 6 znachenie integrala= NaN
alfa= 7 znachenie integrala= NaN
alfa= 8 znachenie integrala= NaN
alfa= 9 znachenie integrala= NaN
alfa= 10 znachenie integrala= NaN

ФедосеевПавел · 17.05.2016, 20:34

Проверьте формулу.
при х=0 получается, что
f:=(sqrt(Power(x,alfa))*ln(alfa*x+1)/sqrt(Power(x,2*alfa))+alfa*Power(x,alfa)+1);
sqrt(Power(x,2*alfa))=0
Деление на 0.

@Sardorxon · 17.05.2016, 20:42 **[ТС]**

проверил ! естественно .. в место / (деление) поставил *
то есть изменил формулу , ответ есть

Добавлено через 22 секунды
f:=(sqrt(Power(x,alfa))*ln(alfa*x+1)*sqrt(Power(x,2*alfa))+alfa*Power(x,alfa) +1);

@Sardorxon · 17.05.2016, 20:45 **[ТС]**

теперь что будем делать ответ есть но функция изменился, но это уже не решения задачи

ФедосеевПавел · 17.05.2016, 20:47

Почему это не решение задачи?

@Sardorxon · 17.05.2016, 20:51 **[ТС]**

вот мой вариант ..

@Sardorxon · 17.05.2016, 20:53 **[ТС]**

а мы изменили задачу..нам надо было изменит команды внутри программы

ФедосеевПавел · 17.05.2016, 20:57

Скобки не на тех местах.

Добавлено через 1 минуту
Не нам, а вам.

Эта задача мне глубоко безразлична. Хотя бы из-за отсутствия какой либо заинтересованности.

@Sardorxon · 17.05.2016, 20:58 **[ТС]**

конечно .. благодарю ) ша поищу лишние скобки ..

ФедосеевПавел · 17.05.2016, 20:59

Не лишние, а не на своих местах.

@Sardorxon · 17.05.2016, 21:02 **[ТС]**

ураа ! спасибо огромное !
правильная формула
f:=(sqrt(Power(x,alfa))*ln(alfa*x+1)/sqrt(Power(x,2*alfa)+alfa*Power(x,alfa)+ 1));

Новые блоги и статьи Все статьи Все блоги /
Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Programma_Boinc 23.12.2025 Thinkpad X220 Tablet — это лучший бюджетный ноутбук для учёбы, точка. Рецензия / Мнение Это мой обзор планшета X220 с точки зрения школьника. Недавно я решила попытаться уменьшить свой. . .	PhpStorm 2025.3: WSL Terminal всегда стартует в ~ and_y87 14.12.2025 PhpStorm 2025. 3: WSL Terminal всегда стартует в ~ (home), игнорируя директорию проекта Симптом: После обновления до PhpStorm 2025. 3 встроенный терминал WSL открывается в домашней директории. . .	Как объединить две одинаковые БД Access с разными данными VikBal 11.12.2025 Помогите пожалуйста !! Как объединить 2 одинаковые БД Access с разными данными.	Новый ноутбук volvo 07.12.2025 Всем привет. По скидке в "черную пятницу" взял себе новый ноутбук Lenovo ThinkBook 16 G7 на Амазоне: Ryzen 5 7533HS 64 Gb DDR5 1Tb NVMe 16" Full HD Display Win11 Pro	Музыка, написанная Искусственным Интеллектом volvo 04.12.2025 Всем привет. Некоторое время назад меня заинтересовало, что уже умеет ИИ в плане написания музыки для песен, и, собственно, исполнения этих самых песен. Стихов у нас много, уже вышли 4 книги, еще 3. . .
От async/await к виртуальным потокам в Python IndentationError 23.11.2025 Армин Ронахер поставил под сомнение async/ await. Создатель Flask заявляет: цветные функции - провал, виртуальные потоки - решение. Не threading-динозавры, а новое поколение лёгких потоков. Откат?. . .	Поиск "дружественных имён" СОМ портов Argus19 22.11.2025 Поиск "дружественных имён" СОМ портов На странице: https:/ / norseev. ru/ 2018/ 01/ 04/ comportlist_windows/ нашёл схожую тему. Там приведён код на С++, который показывает только имена СОМ портов, типа,. . .	Сколько Государство потратило денег на меня, обеспечивая инсулином. Programma_Boinc 20.11.2025 Сколько Государство потратило денег на меня, обеспечивая инсулином. Вот решила сделать интересный приблизительный подсчет, сколько государство потратило на меня денег на покупку инсулинов. . . .	Ломающие изменения в C#.NStar Alpha Etyuhibosecyu 20.11.2025 Уже можно не только тестировать, но и пользоваться C#. NStar - писать оконные приложения, содержащие надписи, кнопки, текстовые поля и даже изображения, например, моя игра "Три в ряд" написана на этом. . .	Мысли в слух kumehtar 18.11.2025 Кстати, совсем недавно имел разговор на тему медитаций с людьми. И обнаружил, что они вообще не понимают что такое медитация и зачем она нужна. Самые базовые вещи. Для них это - когда просто люди. . .

@Sardorxon 0 / 0 / 0 Регистрация: 17.05.2016 Сообщений: 11
	17.05.2016, 19:54 [ТС]
	скриншот могу вложит Миниатюры 0

ФедосеевПавел Модератор 8645 / 4480 / 1669 Регистрация: 01.02.2015 Сообщений: 13,888 Записей в блоге: 11
	17.05.2016, 20:06
	Логарифм 0 - неопределён. Измените формулы вычисления степени или попробуйте штатную для PABC.NET. 1

ФедосеевПавел Модератор 8645 / 4480 / 1669 Регистрация: 01.02.2015 Сообщений: 13,888 Записей в блоге: 11
	17.05.2016, 20:13
	exp(ln(x)*alfa) замените Power(x, alpha) И в нескольких местах аналогично 0

@Sardorxon 0 / 0 / 0 Регистрация: 17.05.2016 Сообщений: 11
	17.05.2016, 20:20 [ТС]
	все равно не па могло прога работает но ответа нету ((( f:=(sqrt(Power(x,alfa))ln(alfax+1)/sqrt(Power(x,2alfa))+alfaPower(x,alfa) +1); Добавлено через 1 минуту alfa= 1 znachenie integrala= NaN alfa= 2 znachenie integrala= NaN alfa= 3 znachenie integrala= NaN alfa= 4 znachenie integrala= NaN alfa= 5 znachenie integrala= NaN alfa= 6 znachenie integrala= NaN alfa= 7 znachenie integrala= NaN alfa= 8 znachenie integrala= NaN alfa= 9 znachenie integrala= NaN alfa= 10 znachenie integrala= NaN 0

ФедосеевПавел Модератор 8645 / 4480 / 1669 Регистрация: 01.02.2015 Сообщений: 13,888 Записей в блоге: 11
	17.05.2016, 20:34
	Проверьте формулу. при х=0 получается, что f:=(sqrt(Power(x,alfa))ln(alfax+1)/sqrt(Power(x,2alfa))+alfaPower(x,alfa)+1); sqrt(Power(x,2*alfa))=0 Деление на 0. 0

@Sardorxon 0 / 0 / 0 Регистрация: 17.05.2016 Сообщений: 11
	17.05.2016, 20:42 [ТС]
	проверил ! естественно .. в место / (деление) поставил * то есть изменил формулу , ответ есть Добавлено через 22 секунды f:=(sqrt(Power(x,alfa))ln(alfax+1)sqrt(Power(x,2alfa))+alfa*Power(x,alfa) +1); 0

@Sardorxon 0 / 0 / 0 Регистрация: 17.05.2016 Сообщений: 11
	17.05.2016, 20:45 [ТС]
	теперь что будем делать ответ есть но функция изменился, но это уже не решения задачи Миниатюры 0

ФедосеевПавел Модератор 8645 / 4480 / 1669 Регистрация: 01.02.2015 Сообщений: 13,888 Записей в блоге: 11
	17.05.2016, 20:47
	Почему это не решение задачи? 0

@Sardorxon 0 / 0 / 0 Регистрация: 17.05.2016 Сообщений: 11
	17.05.2016, 20:51 [ТС]
	вот мой вариант .. Миниатюры 0

@Sardorxon 0 / 0 / 0 Регистрация: 17.05.2016 Сообщений: 11
	17.05.2016, 20:53 [ТС]
	а мы изменили задачу..нам надо было изменит команды внутри программы 0

ФедосеевПавел Модератор 8645 / 4480 / 1669 Регистрация: 01.02.2015 Сообщений: 13,888 Записей в блоге: 11
	17.05.2016, 20:57
	Скобки не на тех местах. Добавлено через 1 минуту Не нам, а вам. Эта задача мне глубоко безразлична. Хотя бы из-за отсутствия какой либо заинтересованности. 0

@Sardorxon 0 / 0 / 0 Регистрация: 17.05.2016 Сообщений: 11
	17.05.2016, 20:58 [ТС]
	конечно .. благодарю ) ша поищу лишние скобки .. 0

ФедосеевПавел Модератор 8645 / 4480 / 1669 Регистрация: 01.02.2015 Сообщений: 13,888 Записей в блоге: 11
	17.05.2016, 20:59
	Не лишние, а не на своих местах. 0

@Sardorxon 0 / 0 / 0 Регистрация: 17.05.2016 Сообщений: 11
	17.05.2016, 21:02 [ТС]
	ураа ! спасибо огромное ! правильная формула f:=(sqrt(Power(x,alfa))ln(alfax+1)/sqrt(Power(x,2alfa)+alfaPower(x,alfa)+ 1)); 0