Распишите пожалуйста эту формулу для паскаля

Question

заранее спасибо большое!

Программа выдаёт ошибку - вещественное деление на 0 =(((

program simps;
var
int, int1, y1, yn, s4, s2, h, del, x, U, V, eps: real;
i, j, k, km, ERR: integer;
function F (x : real) : real;
begin
F := Arctan((exp(power(3*x,-1/3)))/Sqrt(1-sqr((exp(power(3*x,-1/3))))));
end;
begin
int1 := 0;
y1 := F(u);
yn := F(v);
k := 2;
ERR := 1;
for j := 1 to Km do
begin
int := 0;
h := (v - u) / k;
s4 := 0;
s2 := 0;
i := 1;
x := u + h;
repeat
if odd(i) then s4 := s4 + F(x)
else s2 := s2 + F(x);
int := h * (y1 + yn + 4 * s4 + 2 * s2) / 3;
inc(i);
x := x + h;
until i >= k;
del := abs(int1 - int);
if del < eps then
begin
ERR := ERR - 1;
exit
end
else
begin
int1 := int;
k := k * 2;
end;
end;
end.

Юлия Князева · Accepted Answer

Сначала пишем свою функцию на арксинус: 
function Arcsin (x : real) : real; 
begin 
Arcsin := Arctan(x/Sqrt(1-sqr(x))); 
end; 
 
Затем уже работаем с ней: 
otvet:=Arcsin(exp(power(3*x,-1/3)));