Khi nghiên cứu giải quyết các bài toán thực tế trong các môi trường liên tục, thông qua phương pháp mô hình hóa thì đại đa số các bài toán đều đưa đến mô hình được mô tả bởi các phương trình vi phân đạo hàm riêng tức là các mô hình có chứa toán tử vi phân. Một lớp rất nhỏ các bài toán ứng với mô hình và điều kiện biên đơn giản, ta có thể thu được lời giải trực tiếp của bài toán thông qua các phương pháp giải tích, còn đại đa số các bài toán phức tạp đều thông qua các phương pháp rời rạc hóa các toán tử vi phân để chuyển về các hệ phương trình sai phân. Khi đó nghiệm xấp xỉ sẽ thu được thông qua việc giải các hệ phương trình sai phân dựa trên công cụ của máy tính điện tử. Với yêu cầu cần thu được lời giải với độ chính xác cao thì vấn đề nghiên cứu các phương pháp rời rạc hóa các toán tử vi phân với độ chính xác cao là một hướng nghiên cứu được các nhà toán học đặc biệt quan tâm. Trong bài báo này, chúng tôi đề xuất một thuật toán rời rạc hóa đạo hàm các cấp với độ chính xác bậc cao. Các kết quả lý thuyết và tính toán thực nghiệm đã khẳng định độ chính xác của thuật toán.When studying and solving practical problems in continuous environments, through modeling methods, the vast majority of problems lead to models described by partial differential equations, i.e. models that contain differential operators. For a very small class of problems corresponding to simple models and boundary conditions, we can obtain a direct solution of the problem through analytical methods, while the vast majority of complex problems can be obtained through analytical methods. Methods of discretizing differential operators to convert to systems of difference equations. Then the approximate solution will be obtained through solving the system of difference equations based on the tools of an electronic computer. With the need to obtain solutions with high accuracy, the issue of researching methods to discretize differential operators with high precision is a research area of special interest to mathematicians. In this paper, we propose an algorithm to discretize the n-th order derivative with high-order accuracy. Theoretical results and experimental calculations have confirmed the accuracy of the algorithm.