[pva]

Если применить чуток теории и сделать предварительные рассчёты, то задача упрощается:

Код:

(* это скрипт для пакета Wolfram Research Mathematica 4.0 *)
In[1]: tab1 = Prime[Range[2000]];
(* перебираем все 4-значные простые числа, для которых сумма разрядов равна произведению *)
Part[tab1, 
In[2]:  Flatten@Position[tab1, 
      x_ /; Block[{dig1 = RealDigits[x]}, 
          dig1[[2]] == 4 && (Plus @@ dig1[[1]] == Times @@ dig1[[1]])]]]
(* получаем *)
Out[1]: {2141, 2411, 4211}

то есть нужно просто проверить на совпадение с любым из этих чисел ;-)