Abstract. This paper discusses the properties and the numerical discretizations of the fractional substantial integraland the fractional substantial derivativewhere Ds = ∂ ∂x + σ = D + σ, σ can be a constant or a function without related to x, say σ(y); and m is the smallest integer that exceeds µ. The Fourier transform method and fractional linear multistep method are used to analyze the properties or derive the discretized schemes. And the convergences of the presented discretized schemes with the global truncation error O(h p ) (p = 1, 2, 3, 4, 5) are theoretically proved and numerically verified.