Вывести значение функции Y=ARCTG(X)+(x^2) на промежутке [a..b] c шагом 0.5, используя цикл с предусловием — QBasic(Бейсик)

INPUT A, b
DO
        X = A
        DO WHILE X < b
        PRINT "Y ="; 1 / ATN(X) + (X ^ 2)
        X = X + .5
        LOOP
        K = K + 1        
LOOP UNTIL K = 2