Before revision:  In the article it is also shown how to valuate numerically wide range of barrier options using the FD approach . The proposed FD method has 2-\alpha order of accuracy with respect to time, where \alpha\in(0,1) is the subdiffusion parameter, and 2 with respect to space .
After revision:  In the article it is also shown how to valuate numerically wide range of barrier options using the FD approach .