// This MAGMA script makes a list of values of
// a particular linear combination of Bessel functions.

numvalues := 200000;
a := 0;
b := 15;
RR := RealField(50);

vals := [RR|];
for m in [1..numvalues] do
  if (m mod 500) eq 0 then
    printf "Working on value %o.\n",m;
  end if;
  input := RR!((m/numvalues)*(b-a));
  funcval := (-6/Pi(RR))*Sqrt(input)*KBessel2(1,4*Pi(RR)*Sqrt(input)) + 24*input*KBessel2(0,4*Pi(RR)*Sqrt(input));
  Append(~vals,funcval);
end for;
printf "Writing to file...";
PrintFile("beslist.out",Sprintf("vals := "));
PrintFileMagma("beslist.out",vals);
PrintFile("beslist.out",Sprintf(";"));
printf "Done!\n";