Табулирование функции(Prolog)

[Tomogochi]

Добрый день! Задача такова:написать программу,выполняющую табулирование функции.Программу я написала, она работет, но работает неправильно... То есть, у меня есть тесты правильные(от препода), так вот моя прога выдаёт результаты отличные от них. Не мог бы кто-нибудь проверить меня?Где моя ошибка?
Задание:


Моя программа:

domains
r=real
predicates
y(r,r,r)
w
goal
write("X0="),readreal(X0),nl,write( "Xk"),readreal(Xk),nl,write("H="),r eadreal(H),nl,y(X0,Xk,H).
clauses
w:-clearwindow,makewindow(25,26,203,"t abulirovanie",0,0,25,80).
y(X0,Xk,H):-X0>Xk+H/2,write("End").
y(X0,Xk,H):-
X0<Xk+H/2,X0>-0.000000001,X0<0.000000001,write("x =",X0," y=No"),nl,X=X0+H,y(X,Xk,H);
X0<Xk+H/2,X0>1.999999999,X0<2.000000001,wri te("x=",X0," y=No"),nl,X=X0+H,y(X,Xk,H);
X0<Xk+H/2,Z1=X0-2,Z2=X0/Z1,Z3=-X0/Z1,Z4=exp(Z2)-exp(Z3),Z5=exp(X0),Z6=Z5/Z4,write("x=",X0," y=",Z6),nl,X=X0+H,y(X,Xk,H).


Спасибо!

[mirakl]

формула синуса гиперболического расписана не верно... отсутствует деление на 2 "sh(x) = [exp(x) - exp(-x)]/2"

[rrrFer]

write("H="),r eadreal(H),

"readreal" пишется слитно(пробел лишний)

X0<Xk+H/2

эти условия(повторяется 3 раза) можно убрать, они тут лишние

w:-clearwindow,makewindow(25,26,203,"t abulirovanie",0,0,25,80).

вроде бы правило не используется - тоже можно убрать

То есть, у меня есть тесты правильные(от препода)

а у нас нет таких тестов, выкладывайте на форуме, тогда МБ и поможет кто-нибудь

X0>-0.000000001,X0<0.000000001

полагаю, вы тут с нулем сравниваете, но опасаетесь что ноль окажется не равен нулю(ну такое бывает...когда мусор в младших разрядах накапливается), у меня подозрения что в прологе этой проблемы не возникнет(хотя, не проверял). А ближе к теме...лучше не 0.0000(0)1 писать. а нечто типа (H/100) ну и (-H/100) во-первых это красивей чем константы такого чудного вида как у вас, а во вторых представьте что шаг меньше чем 0.000000001... хотя в данном случае, МБ и так сгодится

А еще красивей выделить это хозяйство в функцию сравнения с заданной точностью

Z2=X0/Z1,Z3=-X0/Z1

лучше:

Z2=X0/Z1,Z3=-Z2

А еще лучше вообще не использовать Z3, а просто, вместо него использовать -Z2