prob60(N)=
{
	\\ ab, ba
	\\ ac, ca, bc, cb,
	\\ ad, da, db, bd, dc, cd
	\\ ae, ea, eb, be, ec, ce, de, ed
	default(primelimit,10000000);
	forprime(a=11,N,
		print(a);
		counta = length(Strtex(a));
		forprime(b=7,a-1,
			countb = length(Strtex(b));
			if (isprime(a*10^countb + b) && isprime(b*10^counta+a),
				forprime(c=5,b-1,
					countc = length(Strtex(c));
					if (isprime(a*10^countc + c) && isprime(c*10^counta+a) &&
						isprime(c*10^countb + b) && isprime(b*10^countc+c),
						forprime(d=3,c-1,
							countd = length(Strtex(d));
							if (isprime(a*10^countd + d) && isprime(d*10^counta+a) &&
							isprime(d*10^countb + b) && isprime(b*10^countd+d) &&
							isprime(d*10^countc + c) && isprime(c*10^countd+d),
								forprime(e=2,d-1,
									counte = length(Strtex(e));
									if (isprime(a*10^counte + e) && isprime(e*10^counta+a) &&
									isprime(e*10^countb + b) && isprime(b*10^counte+e) &&
									isprime(e*10^countc + c) && isprime(c*10^counte+e) &&
									isprime(e*10^countd + d) && isprime(d*10^counte+e),
										print1(a);
										print1(", ");
										print1(b);
										print1(", ");
										print1(c);
										print1(", ");
										print1(d);
										print1(", ");
										print(e);
										
										print (a + b + c + d + e);
										quit;
									);
								);
							);
						);
					);
				);
			);
		);
	);
}