Интегральный синус

Psilon · Регистрация: 10.07.2011

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

День добрый.
Собственно, элементарная программа рассчета интегрального синуса, но что-то у меня руки-крюки

Хотел использовать рекуррентную формулу
$https://www.cyberforum.ru/cgi-bin/latex.cgi?{a}_{0} = x$
$https://www.cyberforum.ru/cgi-bin/latex.cgi?{a}_{n + 1} = -{a}_{n} * {x}^{2}* \frac{2n+1}{{2(2n+3)}^{2}(n+1)}$
Но программа не работает. Подскажите, где косяк?

Заранее спасибо.

Prolog
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
predicates  
    si(real, real, integer, real,real).
    si(real, integer, real).
clauses
    si(X, N, R) :- si(X, X, N, 0, R).
    si(_, _, 0, R, R) :- !.
    si(X, A, N, Temp, R) :-     X2 = X*X,
                        A1 = 2*N+1,                     
                        A2 = 2*N+3,
                        N1 = N + 1,
                        NM1 = N - 1,
                        A3 = 2*A2*A2*N1,
                        
                        AL = -A*X2*A1/A3,
                        R = Temp + AL,
                        si(X,AL,NM1,R,R).

Добавлено через 13 минут
Хотя я тут посмотрю, вообще косяков миллион, например считаем старшие N для младших степеней...

Есть какой-то более адекватный способ, чем лепить предикат с миллионом аргументов?..

Добавлено через 18 часов 53 минуты
Исправленная программа в первом приближении работает ожидаемо, но вот при N > 2 почему-то отказывается выдавать результат. Как дебажить пролог-программы понятия не имею. Есть варианты?..

Prolog
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
predicates  
    si(real, real, integer, integer, real, real).
    si(real, integer, real).
clauses
    si(X, N, R) :- si(X, X, 1, N, X, R).
    si(_, _, N, N, R, R) :- !.
    si(X, A, CurN, N, PrevR, R) :-  X2 = X*X,
                        A1 = 2*CurN+1,                      
                        A2 = 2*CurN+3,
                        N1 = CurN + 1,
                        A3 = 2*A2*A2*N1,
                        
                        AL = -A*X2*A1/A3,
                        R = PrevR + AL,
                        si(X,AL,N1,N,R,R).

@Catstail · 16.05.2014, 20:05

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

Как дебажить

- в ТП можно попробовать вставить отладочную печать. Или включить трассировку.

Psilon · 16.05.2014, 20:47 **[ТС]**

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

спасибо, действительно, надо попробовать.

Psilon · 16.05.2014, 22:09 **[ТС]**

Кто может подсказать, пишу

Prolog
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
predicates  
    si(real, real, integer, integer, real, real).
    si(real, integer, real).
clauses
    si(X, N, R) :- si(X, X, 0, N, X, R).
    si(_, _, N, N, R, R) :- !.
    si(X, A, CurN, N, PrevR, R) :-      
                        X2 = X*X, write("X2= ") ,write(X2),write(" "),              
                        A1 = 2*CurN+1, write("A1= ")    ,write(A1),write(" ")   ,               
                        A2 = 2*CurN+3,write("A2= ") ,write(A2),write(" ")   ,   
                        N1 = CurN + 1,write("N1= ") ,write(N1),write(" ")   ,   
                        A3 = 2*A2*A2*N1,write("A3= ")   ,write(A3),write(" "),      
                        
                        AL = -A*X2*A1/A3,write("AL= ")  ,write(AL),write(" ")   ,   
                        R = PrevR + AL,write("R= ") ,   write(R), nl,
                        si(X,AL,N1,N,R,R).

вывод такой:

почему он дальше не считает и говорит NO SOLUTION?..

Psilon · 16.05.2014, 22:18 **[ТС]**

То есть теоретически он должен сложить верхний R и нижний AL, получаем нижний R.

Добавлено через 3 минуты
Даже такой код:

Prolog
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
predicates  
    si(real, real, integer, integer, real, real).
    si(real, integer, real).
clauses
    si(X, N, R) :- si(X, X, 0, N, X, R).
    si(_, _, N, N, R, R) :- !.
    si(X, A, CurN, N, PrevR, R) :-      
                        X2 = X*X,           
                        A1 = 2*CurN+1,          
                        A2 = 2*CurN+3,  
                        N1 = CurN + 1,
                        A3 = 2*A2*A2*N1,
                        
                        AL = -A*X2*A1/A3,
                        write("PrevR= ")    ,write(PrevR),write(" || ") ,   
                        Temp = PrevR + AL, write("Temp= ")  ,write(Temp),
                        write(" |---| ")    ,
                        R = Temp, write("R= ")  ,   write(R), nl,
                        si(X,AL,N1,N,R,R).

он корректно считает Temp, но не может сделать R = temp...

@Black Fregat · 17.05.2014, 00:40

Вы перепутали "итоговый" R, который должен получиться в результате, который сопоставляется один раз в первом правиле по достижении конца цикла, и "новый" R, который вычисляется заново на каждом шаге. Они обозначены одной и той же буквой R, а это неправильно. Нужно разнести примерно так:

Prolog
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
predicates  
    si(real, real, integer, integer, real, real).
    si(real, integer, real).
clauses    si(X, N, R) :- si(X, X, 0, N, X, R).
    si(_, _, N, N, R, R) :- !.
    si(X, A, CurN, N, PrevR, R) :-      
                        X2 = X*X, write("X2= ") ,write(X2),write(" "),              
                        A1 = 2*CurN+1, write("A1= ")    ,write(A1),write(" ")   ,               
                        A2 = 2*CurN+3,write("A2= ") ,write(A2),write(" ")   ,   
                        N1 = CurN + 1,write("N1= ") ,write(N1),write(" ")   ,   
                        A3 = 2*A2*A2*N1,write("A3= ")   ,write(A3),write(" "),      
                        
                        AL = -A*X2*A1/A3,write("AL= ")  ,write(AL),write(" ")   ,   
                        NewR = PrevR + AL,write("NewR= ") ,   write(NewR), nl,
                        si(X,AL,N1,N,NewR,R).

Psilon · 17.05.2014, 11:22 **[ТС]**

Black Fregat, огромное спасибо, буду знать! )

Добавлено через 1 минуту
посмотрел код, понял ошибку. Еще раз спасибо!

Новые блоги и статьи Все статьи Все блоги /
Как я обхитрил таблицу Word Alexander-7 21.03.2026 Когда мигает курсор у внешнего края таблицы, и нам надо перейти на новую строку, а при нажатии Enter создается новый ряд таблицы с ячейками, то мы вместо нервных нажатий Энтеров мы пишем любые буквы. . .	Krabik - рыболовный бот для WoW 3.3.5a AmbA 21.03.2026 без регистрации и смс. Это не торговля, приложение не содержит рекламы. Выполняет свою непосредственную задачу - автоматизацию рыбалки в WoW - и ничего более. Однако если админы будут против -. . .	Программный отбор значений справочника Maks 21.03.2026 Установка программного отбора значений справочника "Сотрудники" из модуля формы документа. В качестве фильтра для отбора служит предопределенное значение перечислений. Процедура. . .	Переходник USB-CAN-GPIO Eddy_Em 20.03.2026 Достаточно давно на работе возникла необходимость в переходнике CAN-USB с гальваноразвязкой, оный и был разработан. Однако, все меня терзала совесть, что аж 48-ногий МК используется так тупо: просто. . .
Оттенки серого Argus19 18.03.2026 Оттенки серого Нашёл в интернете 3 прекрасных модуля: Модуль класса открытия диалога открытия/ сохранения файла на Win32 API; Модуль класса быстрого перекодирования цветного изображения в оттенки. . .	SDL3 для Desktop (MinGW): Рисуем цветные прямоугольники с помощью рисовальщика SDL3 на Си и C++ 8Observer8 17.03.2026 Содержание блога Финальные проекты на Си и на C++: finish-rectangles-sdl3-c. zip finish-rectangles-sdl3-cpp. zip	Символические и жёсткие ссылки в Linux. algri14 15.03.2026 Существует два типа ссылок — символические и жёсткие. Ссылка в Linux — это запись в каталоге, которая может указывать либо на inode «файла-ИСТОЧНИКА», тогда это будет «жёсткая ссылка» (hard link),. . .	[Owen Logic] Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ФедосеевПавел 14.03.2026 Поддержание уровня воды в резервуаре количеством включённых насосов: моделирование и выбор регулятора ВВЕДЕНИЕ Выполняя задание на управление насосной группой заполнения резервуара,. . .

Psilon Master of Orion 6102 / 4958 / 905 Регистрация: 10.07.2011 Сообщений: 14,522 Записей в блоге: 5
	16.05.2014, 20:47 [ТС]
	Catstail, мда, после нормальных IDE, мысль об отладочной печати даже в голову не пришла спасибо, действительно, надо попробовать. 0

Psilon Master of Orion 6102 / 4958 / 905 Регистрация: 10.07.2011 Сообщений: 14,522 Записей в блоге: 5
	17.05.2014, 11:22 [ТС]
	Black Fregat, огромное спасибо, буду знать! ) Добавлено через 1 минуту посмотрел код, понял ошибку. Еще раз спасибо! 0

Интегральный синус

Решение