Еще не увидела, где вы находите значение функции на концах отрезка.
Например, для у=0 хє[-1; 1], поэтому в этом случае еще надо найти z(-1), z(1). Аналогично и для дргих участков.
Посмотрите примеры