BW [p_]  :=  3 (1 - p) p^2  +  p^3

g1 =  Plot [Log [10, BW[10^(-x)]],  {x,  .2, 6}]

BC11 [p_]  :=  
      1  -  (1 - p)^11  -  11 (1 - p)^10 *p  -  
        55 (1 - p)^9 *p^2

BC11 [p]

g2  =  
      Plot [Log [10, BC11[10^(-x)]],  {x,  .2, 5}]

Show [{g1,  g2}]

120 !

GCD [64,  13]

ExtendedGCD [64, 13]

e  =  ExtendedGCD [11121, 34343]

l =  {a, 1, 2}

l2  =  Append [l, t]

l

l3  =  Append [l2,  s]

{a, 1, 2, t, s} [[3]]

e  =  
      ExtendedGCD [34354, 2222]

e [[2]] [[1]]

<<  poliesModP.m

zweierlog [
        n_]  := 
Module[{l,  
          n2}, 
l  =  0; 
n2  =  
          n; 
While [
          n2  >  1, 
n2  =  
            Floor [n2/2]; 
l  =  
            l  +  1]; 
l]

zweierlog [8]

fac [n_]  :=  If [n  ≤ 0, 1, n*fac[n - 1]]

fac[3]

y  =  -3

Plot[Sin [x] *y,  {x, 1, 16}]

ExtendedGCD [2^1000 - 1,  2^1030  -  1]

Mod [7, 4]

Mod [555, 7]

PowerMod [555, 1000000000000000000000010, 7]

PowerMod [5, -1, 7]