\r jsloglib_v0.5.gp
\p 48
\ps 64
{
  InitFixedEstimator(48, 128);
  for (k=1,10,
    lnb = k/5. + max(0, (k-5)/5.)*I;
    print("log(B)= ",lnb);
    print("Base  = ",exp(lnb));
    print("L1(B) = ", FindFixed_lnB(lnb, -1));
    print("L2(B) = ", FindFixed_lnB(lnb,  1));
    print();
  );
    
    b = sqrt(2.);
    lnb = log(b);
    print("log(B)= ",lnb);
    print("Base  = ",b);
    print("L1(B) = ", FindFixed_lnB(lnb, -1));
    print("L2(B) = ", FindFixed_lnB(lnb,  1));
    print();
    
    b = 1.001;
    lnb = log(b);
    print("log(B)= ",lnb);
    print("Base  = ",b);
    print("L1(B) = ", FindFixed_lnB(lnb, -1));
    print("L2(B) = ", FindFixed_lnB(lnb,  1));
    print();
    
    b = -1.;
    lnb = log(b);
    print("log(B)= ",lnb);
    print("Base  = ",b);
    print("L1(B) = ", FindFixed_lnB(lnb, -1));
    print("L2(B) = ", FindFixed_lnB(lnb,  1));
    print();

    b = 0.001;
    lnb = log(b);
    print("log(B)= ",lnb);
    print("Base  = ",b);
    print("L1(B) = ", FindFixed_lnB(lnb, -1));
    print("L2(B) = ", FindFixed_lnB(lnb,  1));
    print();
    
    b = 0.318131505204764135312654251587664517203517613871 + 1.33723570143068940890116214319371061253950213846*I;
    print("Base  = ",b);
    print("L1(B) = ", FindFixed_B(b, -1));
    print("L2(B) = ", FindFixed_B(b,  1));
    print();
}