In this paper we constructed optimal with respect to order quadrature formulas for evaluating one- and multidimensional hypersingular integrals on classes of functions Ωur,γ(Ω,M), Ω¯ur,γ(Ω,M), Ω=[−1,1]l, l=1,2,…,M=Const, and γ is a real positive number. The functions that belong to classes Ωur,γ(Ω,M) and Ω¯ur,γ(Ω,M) have bounded derivatives up to the rth order in domain Ω and derivatives up to the sth order (s=r+⌈γ⌉) in domain Ω∖Γ, Γ=∂Ω. Moduli of derivatives of the vth order (r<v≤s) are power functions of d(x,Γ)−1(1+|lnd(x,Γ)|), where d(x,Γ) is a distance between point x and Γ. The interest in these classes of functions is due to the fact that solutions of singular and hypersingular integral equations are their members. Moreover various physical fields, in particular gravitational and electromagnetic fields belong to these classes as well. We give definitions of optimal with respect to accuracy methods for solving hypersingular integrals. We constructed optimal with respect to order of accuracy quadrature formulas for evaluating one- and multidimensional hypersingular integrals on classes of functions Ωur,γ(Ω,M) and Ω¯ur,γ(Ω,M).