from scipy import *
import sys


def calc(x,k):
  x_0 =pi/3 # minimum at 60 deg
  if x<x_0: # Put divergence at 0 deg
    x_m =pi/3
  else:
    x_m =2*pi/3

  # Return force and energy
  return \
    k *(x-x_0)*x_m**2 /(x**2 -2*x*x_0 +x_0**2 -x_m**2),\
    -k *x_m**2 *log(1-((x-x_0)/x_m)**2)



n =100000
k=float(sys.argv[1])

print "#"+str(n),0,pi
for x in linspace(0,pi,n):
 if abs(x)<1E-10:
   print 0.,100000000000,100000000000
   continue
 if abs(x-pi)<1E-10:
   print 0.,-100000000000,100000000000
   continue
 res=calc(x,k)
 print x,res[0],res[1]