1) Что это y1:=xmin+(xmin/0.001) ?
2) Почему dx := 0.1 ?
3) В ваших переменных sqrt(x) = sqrt(z)
4) Зачем переименовали переменные?? Проверять тяжело.