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