RosettaCodeData/Task/Bernoulli-numbers/Phix/bernoulli-numbers-2.phix
2016-12-05 23:44:36 +01:00

10 lines
319 B
Text

sequence a = {}
for m=0 to 60 do
a = append(a,{ba_new(1),ba_new(m+1)})
for j=m to 1 by -1 do
a[j] = ba_frac_mul({ba_new(j),ba_new(1)},ba_frac_sub(a[j+1],a[j]))
end for
if a[1][1]!=BA_ZERO then
printf(1,"B(%2d) = %44s / %s\n",{m,ba_sprint(a[1][1]),ba_sprint(a[1][2])})
end if
end for