Предложены разностные схемы для уравнений Фоккера-Планка-Колмогорова и уравнений Понтрягина, позволяющие численно находить статистические характеристики систем синхронизации в переходном режиме при различных дискриминационных характеристиках фазового детектора и при ненулевой частотной расстройке между частотами управляемого генератора и сигнала, что невозможно выполнить с помощью известных аналитических процедур. Приведена процедура оценки сходимости функционального ряда для указанной плотности распределения вероятностей в стационарном режиме, которая служит для оценки точности численного метода. Показано, что предлагаемые процедуры численного решения не требуют значительных затрат машинного времени и могут быть реализованы на стандартных ПЭВМ.