Численное интегрирование(метод Симпсона)

[km77]

Есть программа,которая должнна вычислять интеграл ф-ии x^2/x-3 на отрезке [1,4].Но есть одна проблема в то что при задании пределов интегрирование от 1 до 3 или от 2 до 3,программа не работает,потому что если посмотреть на ф-ию и поставить вместо х=3,то прогр. не будет считать интеграл,така как знаменатель будет =0,то есть нужно ввести условие при котором она будет работать на всем отрезке от 1 до 4.

program test;
uses crt;
var a,b,r,r2,eps:real;
n:integer;
function f(x:real):real;
begin
f:=sqr(x)/(x-3);
end;
function Simpson(a,b:real; n:integer):real;
var s,h:real;
m,mn:integer;
begin
h:=(b-a)/(2*n);
s:=f(a)+f(b);
mn:=4;
for m:=1 to 2*n-1 do
begin
s:=s+mn*f(a+h*m);
if mn=4 then mn:=2
else mn:=4;
end;
Simpson:=s*h/3;
end;
begin
writeln ('введите нижнюю границу');
Readln (a);
writeln('введите верхнюю границу');
readln(b);
writeln('введите погрешность');
readln(eps);
n:=2;
r:=simpson(a,b,n);
repeat
r2:=r;
n:=n*2;
r:=simpson(a,b,n);
until abs(r-r2)<eps;
Writeln;
Writeln ('Значение интеграла = ' ,r:6:6);
writeln ('количество разбиений = ',n);
readkey;
end.

[KobolD]

function f(x:real):real;
 begin
 if(x=3) return 0; // Если делитель равен нулю то функция равна 0 (пропуск точки)
 f:=sqr(x)/(x-3);
 end;