A088460型:=proc(n)局部i,p;i:=1;p:=0;while true do while ithprime(i+1)-ithprime(i)<>2 do i:=i+1;od;p:=p+1;如果p=n,则返回(ithprime(i));fi;i:=i+1;od;结束时间:A082674号:=proc(n)局部rhs,co,row,col;rhs:=直线[矢量](n+1);co:=linalg[矩阵](n+1,n+1);对于从1到n+1的行,执行rhs[row]:=A088460型(行);对于从1到n+1的列,做co[row,col]:=行^(col-1);od;od;linalg[linsolve](co,rhs)[1];结束:对于从1到30的n,执行printf(“%d,”,A082674号(n) );od#R.J.马塔尔2006年10月31日