MSC: 68Q25 65C30Keywords: Stochastic Itô integrals Singular problems Optimal algorithm Standard information r-fold integrated Brownian motion a b s t r a c tWe deal with numerical approximation of stochastic Itô integrals of singular functions. We first consider the regular case of integrands belonging to the Hölder class with parameters r and . We show that in this case the classical Itô-Taylor algorithm has the optimal error Θ(n −(r+ ) ). In the singular case, we consider a class of piecewise regular functions that have continuous derivatives, except for a finite number of unknown singular points. We show that any nonadaptive algorithm cannot efficiently handle such a problem, even in the case of a single singularity. The error of such algorithm is no less than n − min{1/2,r+ } . Therefore, we must turn to adaptive algorithms. We construct the adaptive Itô-Taylor algorithm that, in the case of at most one singularity, has the optimal error O(n −(r+ ) ). The best speed of convergence, known for regular functions, is thus preserved. For multiple singularities, we show that any adaptive algorithm has the error Ω(n − min{1/2,r+ } ), and this bound is sharp.