Найти все такие числа N (N<=100), что у числа N! Сумма цифр – простое число

for i:=1 to n do
write(a[i],'   ');

{функция для нахождения факториала}
{Функция проверки на "простоту" IsSimple(s:integer):boolean}
begin
{ввод данных}
for i:=2 to N do
{нахождение суммы цифр s для факториала N}
if IsSimple(s)=true then {какое действие добавить тут или где-либо в программе?}
end.

if IsSimple(s)=true then writeln (s);

if IsSimple(s)=true then writeln (a[i]);

uses crt;
var a,b,c:array[byte] of byte; {размерность массива(количество знаков в числе) можно и больше}
    na,nb,i,j,n,p:byte;
begin
clrscr;
randomize;
write('Kolichestvo znakov  chisla A na=');
readln(na);
write('Kolichestvo znakov chisla B nb=');
readln(nb);
writeln('Chislo A:');
a[1]:=random(9)+1;{первая цифра не ноль}
for i:=2 to na do
a[i]:=random(10);
for i:=1 to na do
write(a[i]);
writeln;
writeln('Chislo B:');
b[1]:=random(9)+1;
for i:=2 to nb do
b[i]:=random(10);
for i:=1 to nb do
write(b[i]);
writeln;
if na>nb then {если числа разной длины, то в начало короткого добавляем ноли, чтоб уравнять длину}
  begin
    n:=na;{запоминаем новую общую длину}
    for i:=1 to na-nb do
      begin
        nb:=nb+1;
        for j:=nb downto 2 do{сдвигаем массив на 1 вправо}
        b[j]:=b[j-1];
        b[1]:=0;{вперед добавляем ноль}
      end;
   end
else if nb>na then
   begin
     n:=nb;
     for i:=1 to nb-na do
       begin
         na:=na+1;
         for j:=na downto 2 do
         a[j]:=a[j-1];
         a[1]:=0;
       end;
    end
else if nb=na then n:=na;{если длина одинаковая}
p:=0; (перенос)
for i:=n downto 1 do {начиная с конца складываем столбиком}
  begin
   c[i]:=a[i]+b[i]+p;{сумма=цифра 1 числа+цифра 2 числа + перенос}
   if (c[i]>9)and(i>1) then{если число больше 9,}
      begin
       c[i]:=c[i] mod 10;{ то записываем последнюю цифру}
       p:=1; {один переносим}
      end
    else p:=0;{иначе ничего не переносим}
  end;
if c[1]>9 then {если первое число больше 9}
  begin
    c[1]:=c[1] mod 10;{записываем вторую цифру}
    n:=n+1; {увеличиваем длину}
    for j:=n downto 2 do{сдвигаем вправо}
    c[j]:=c[j-1];
    c[1]:=1;{добавляем в начало 1}
  end;
writeln('Chislo C:');
for i:=1 to n do
write(c[i]);
readln
end.