Решение во вложении ниже

небольшие пометки:
1. писал в ABC.net паскале команда power(x,2) - возведение Х во 2 степень, если писать в других оболочках паскаля, возможно можно будет использовать запись x^2

2. Формула сферы (x-x0)²+(y-y0)²+(z-z0)² = R²
где x0, y0, z0 - координаты начала (т.к. у нас начало в центре получаются эти координаты 0 0 0) и получаем формулу сферы x²+y²+z²=R²
причем эта формула для поверхности сферы, по условию задачи нужно найти точки внутри сферы, поэтому ставим не =, а < (если же понадобятся еще точки лежащие на поверхности сферы, нужно из < сделать <=)