Generate truth tables in Python
pip install truth-table-generator
>>> import ttg >>> print(ttg.Truths([‘a’, ‘b’, ‘c’, ‘d’], [‘(a and b) or not (c xor d)’]))
Avoiding loss of precision with multiple partial sums
Days between two dates
import sys
from datetime import date
a = sys.argv[1]
b = sys.argv[2]
[m1,d1,y1] = a.split('/')
[m2,d2,y2] = b.split('/')
y1 = int(y1)
m1 = int(m1)
d1 = int(d1)
y2 = int(y2)
m2 = int(m2)
d2 = int(d2)
t1=date(y1,m1,d1)
t2=date(y2,m2,d2)
print((t2-t1).days)
Get constellation name from given celestial coordinates
#!/usr/bin/python3 import sys from astropy import units as u from astropy.coordinates import SkyCoord, get_constellation coord=sys.argv[1:] co=" ".join(coord) c=SkyCoord(co, unit=(u.hourangle,u.deg)) print(get_constellation(c,short_name=False))
Print all primitive Pythagorean triples up to n
#!/usr/bin/python3 import sys,math n = int(sys.argv[1]) print([(x,y,z) for x in range(1,n+1) for y in range(x,n+1,1) for z in range(y,n+1,1) if x**2 + y**2 == z**2 and math.gcd(x,y) == 1])
Express a rational number as the sum of unit fractions
#!/usr/bin/python3
import sys
from fractions import Fraction as frac
from sympy import Rational
from sympy.ntheory.egyptian_fraction import egyptian_fraction
n = int(sys.argv[1])
d = int(sys.argv[2])
rat=frac(n,d)
egypt=egyptian_fraction(Rational(rat))
print(str(rat) + ' = 1/' + str(egypt[0]), end = )
for i in range(1,len(egypt)):
print(' + 1/',str(egypt[i]), sep ="", end = )
print ()
Approximate pi from increasingly better fractions
#!/usr/bin/python3
from fractions import Fraction
import sys
import numpy as np
x = np.pi
d = 10
while d < 1000000000:
f=x, Fraction.from_float(x).limit_denominator(d)
a=f[1]
b=str(a)
c=eval(b)
print ("Denom. limit=",d)
print (a,"=",c)
d = d * 10
Normalized (red) and unnormalized (blue) sinc function
#!/usr/bin/python3 import numpy as np import matplotlib.pyplot as plt x=np.linspace(-np.pi*6,np.pi*6,1000) y1=np.sinc(x) y2=np.sinc(x/np.pi) plt.plot(x,y1,color="blue") plt.plot(x,y2,color="red") plt.show()
Plotnorm
Plot the normal distribution for a mean of zero and a standard deviation of one (blue) and two (red):
#!/usr/bin/python3 import numpy as np import matplotlib.pyplot as plt from scipy.stats import norm x1 = np.arange(-3, 3, 0.001) x2 = np.arange(-3, 3, 0.001) plt.plot(x1, norm.pdf(x1, 0, 1)) plt.plot(x2, norm.pdf(x2, 0, 2), color='red', linewidth=3) plt.show()
Capone
Make the first field of standard input uppercase
#!/usr/bin/python3
import sys,re
for line in sys.stdin:
word=line.split()
word1=re.split("\s",line)[0]
word1=word1.upper()
print(str(word1)," ".join(word[1:]),end=" ")
print()
Print a range of binary numbers in #Python
- python -c “for i in range(0,32): print i, (format(i,’b’))”
#!/usr/bin/python
import sys
a = int(sys.argv[1])
b = int(sys.argv[2])
for i in range (a, b+1):
print(format (i, 'b'))
Find anagrams of a word in a list
#!/usr/bin/python
from collections import Counter
import sys
s=sys.argv[1]
def anagrams(word, words):
result = []
counter_word = Counter(word)
for other_word in words:
if Counter(other_word) == counter_word:
result.append(other_word)
return result
with open("words.txt") as wordfile:
lines=wordfile.read().splitlines()
print(anagrams(s,lines))
#!/usr/bin/python
from collections import Counter
def anagrams(word, words):
result = []
counter_word = Counter(word)
for other_word in words:
if Counter(other_word) == counter_word:
result.append(other_word)
return result
with open("words.txt") as wordfile:
lines=wordfile.read().splitlines()
for x in lines:
ana=anagrams(x,lines)
if len(ana) >1:
print(str(ana))
Select mean, standard deviation and count and plot the spread
#!/usr/bin/python3
#USAGE: ./plotrand mean standard-deviation count
import sys,pylab,random
mean=float(sys.argv[1])
dev=float(sys.argv[2])
count=int(sys.argv[3])
x = [random.gauss(mean,dev) for i in range (count)]
y = [random.gauss(mean,dev) for i in range (count)]
pylab.plot(x,y,'bo')
pylab.savefig('points.png')
Lotto
#!/usr/bin/python3
import random
random.seed()
items = [*range(1,50)]
for i in range (0,10):
print (random.sample(items,6))
Plot the Lorentz factor against velocity
#!/usr/bin/python3
import matplotlib.pyplot as plt
import numpy as np
fig = plt.figure()
ax = fig.add_subplot(1, 1, 1)
major_ticks = np.arange(0, 101, 5)
minor_ticks = np.arange(0, 101, 1)
ax.set_xticks(minor_ticks, minor=True)
ax.set_yticks(minor_ticks, minor=True)
ax.grid(which='minor', alpha=0.3)
ax.grid(which='major', alpha=0.8)
x = np.arange(0.7, 1, 0.001)
y = 1/(np.sqrt((1-x**2)))
plt.plot(x*100, y)
plt.title("Lorentz Factor")
plt.xlabel("Real velocity (percentage of c)")
plt.ylabel("Apparent velocity (multiple of c)")
plt.show()
Pairpal
Find pairs of words which form palindromes and write to standard output
#!/usr/bin/python
import sys
file = open("words.txt","r")
words = file.read()
data = words.split("\n")
for y in data:
for x in data:
z = "".join([y,x])
w = z[:][::-1]
if w == z:
v = "".join([y," ",x])
sys.stdout.write('\n'+v)
Plot some nearby stars
#!/usr/bin/python3
import matplotlib.pyplot as plt
types = ['Sol', 'A Centauri', 'Epsilon Eridani', '61 Cygni', 'Procyon', 'Epsilon Indi', 'Tau Ceti', '40 Eridani']
x_coords = [0,-1.6,6.2,6.4,-4.8,5.7,10.3,7.2]
y_coords = [0,-1.4,8.2,-6.1,10.3,-3.2,5.0,14.6]
for i,type in enumerate(types):
x = x_coords[i]
y = y_coords[i]
plt.scatter(x, y, marker='x', color='red')
plt.text(x+0.3, y+0.3, type, fontsize=9)
plt.show()
Convert float to fractions
#!/usr/bin/python3
from fractions import Fraction
import sys
d = float(sys.argv[1])
f=d, Fraction.from_float(d).limit_denominator(1000)
a=f[1]
b=str(a)
c=eval(b)
print ("Decimal=",d)
print ("Fraction:",a,"=",c)
print ("Error=",c-d)
Abundant, deficient, and perfect numbers
#!/usr/bin/python3
import sys
beg=int(sys.argv[1])
end=int(sys.argv[2])
tt = []
for n in range(beg, end+1):
for x in range(1,1+n//2):
if n%x == 0:
tt.append(x)
if sum(tt) == n:
print(str(n) + " is perfect")
elif sum(tt) > n:
print(str(n) + " is abundant")
elif sum(tt) < n:
print(str(n) + " is deficient")
tt = []
Euler’s constant:
#!/usr/bin/python3
import math
import numpy as np
e0 = 0
e = 2
n = 0
fact = 1
while(e-e0>np.nextafter(0,1)):
e0=e
n+=1
fact*=2*n*(2*n+1)
e+=(2.*n+2)/fact
print(e),
print ("Computed e = "+str(e))
print ("Real e = "+str(math.e))
print ("Error = "+str(math.e-e))
Convert Arabic numerals to Roman
#!/usr/bin/python3
import sys
x=int(sys.argv[1])
anums = [1000, 900, 500, 400, 100, 90, 50, 40, 10, 9, 5, 4, 1]
rnums = "M CM D CD C XC L XL X IX V IV I".split()
ret = []
for a,r in zip(anums, rnums):
n,x = divmod(x,a)
ret.append(r*n)
print (.join(ret))
maze
#!/usr/bin/python3
from random import shuffle, randrange
def make_maze(w = 20, h = 10):
vis = [[0] * w + [1] for _ in range(h)] + [[1] * (w + 1)]
ver = [["| "] * w + ['|'] for _ in range(h)] + [[]]
hor = [["+--"] * w + ['+'] for _ in range(h + 1)]
def walk(x, y):
vis[y][x] = 1
d = [(x - 1, y), (x, y + 1), (x + 1, y), (x, y - 1)]
shuffle(d)
for (xx, yy) in d:
if vis[yy][xx]: continue
if xx == x: hor[max(y, yy)][x] = "+ "
if yy == y: ver[y][max(x, xx)] = " "
walk(xx, yy)
walk(randrange(w), randrange(h))
s = ""
for (a, b) in zip(hor, ver):
s += .join(a + ['\n'] + b + ['\n'])
return s
if __name__ == '__main__':
print(make_maze())
Fourier analysis
#!/usr/bin/python3
import numpy as np
import matplotlib.pyplot as plot
t=np.arange(0, 6.28, .01);
a=np.rint(np.sin(t))
plot.subplot(211)
plot.plot(t,a,'bo')
plot.xlabel('Time')
plot.ylabel('Amplitude')
plot.subplot(212)
plot.magnitude_spectrum(a,Fs=4)
plot.show()
fibonacci
fibseq=[1,1,]
fiblen=40
for x in range(1,fiblen-1):
xcount=fibseq[x-1]+fibseq[x]
fibseq.append(xcount)
print(xcount)
Plot a set of data points in #Python
#!/usr/bin/python3
x = [60,65,47,88,62,12,32,72,38,100,87,99,35,93,58]
y = [40,41,45,11,63,67,83,30,57,23,3,96,25,48,54]
import pylab
pylab.plot(x,y,'bo')
pylab.savefig('xypoints.png')
Print a calendar
Plot the sinc function
#!/usr/bin/python3
import numpy as np
import matplotlib.pyplot as plt
in_array=np.linspace(-np.pi,np.pi,100)
out_array=np.sinc(in_array)
plt.plot(in_array,out_array,color="blue", marker="o")
plt.title("Sinc function")
plt.xlabel("X")
plt.ylabel("Y")
plt.show()
Calculate days until Christmas
#!/usr/bin/python3
import datetime
today=datetime.date.today()
xmas=datetime.date(2022,12,25)
diff=xmas-today
d=diff.days
print ("Only",d,"shopping days until Christmas.")
Return all primes over an interval:
#!/usr/bin/python3 import sys m=int(sys.argv[1]) n=int(sys.argv[2]) primes = [i for i in range(m,n) if all(i%j !=0 for j in range(2,int(i**0.5) + 1))] print(primes)
Use list comprehension to make a #Python one-liner that prints the prime numbers from 1 to 5000:
python -c “print([i for i in range(1,5000) if all(i%j !=0 for j in range(2,int(i**0.5) + 1))])”
Plot a function in three dimensions
#!/usr/bin/python3
from sympy import symbols,cos,sin
z = symbols('z')
from sympy.plotting import plot3d_parametric_line
plot3d_parametric_line(cos(z)*sin(z/20),sin(z)*sin(z/20),z,(z,-31.41,31.41))
Max, min, mean, and standard deviation
#!/usr/bin/python3
import numpy
file = open("sample.txt","r")
data = file.read().splitlines()
elements = list(numpy.array(data, dtype = 'float'))
max = max(elements)
min = min(elements)
mean = round((numpy.mean(elements)),5)
std = round(numpy.std(elements),5)
print(f"Max = {max}")
print(f"Min = {min}")
print(f"Mean = {mean}")
print(f"Standard Deviation = {std}")
ABSMAG
Calculate absolute magnitude from stellar apparent visual magnitude and radial distance (LY)
#!/usr/bin/python
import math
import sys
Mvis = float(sys.argv[1])
Dist = float(sys.argv[2])
Vabs = Mvis - 5 * (math.log10(Dist / 3.2616) -1)
print round(Vabs,3)
AVERAGE
Return the arithmetic mean of the arguments:
Excel:
Python:
data =[3546,3027,22,99,3208,348,77,1398,8567,2866,99,66,4272,2784,4429,88,3286,86,55,1766]
from numpy import mean
elements = [float(elements) for elements in data]
print (mean(elements))
PYPRIME
In Python return all primes over an interval:
#!/usr/bin/python3 import sys m=int(sys.argv[1]) n=int(sys.argv[2]) primes = [i for i in range(m,n) if all(i%j !=0 for j in range(2,int(i**0.5) + 1))] print(primes)
TRIG
Why the moon and the sun seem to be the same size in the sky:
import math
Diameter_moon=2159
Distance_moon=238855
Diameter_sun=864938
Distance_sun=92955807
print (math.degrees(math.atan2(Diameter_moon,Distance_moon)))
print (math.degrees(math.atan2(Diameter_sun,Distance_sun)))
With Python find the shortest routes in a network of nearby stars:
#!/usr/bin/python
import math
a=["Sol","A Centauri","E Eridani","61 Cygni","E Indi","Tau Ceti","40 Eridani"]
x=[0,-1.643,6.197,6.489,5.658,10.28,7.195]
y=[0,1.374,8.294,-6.109,-3.157,5.018,14.63]
z=[0,-3.838,-1.725,7.152,-9.896,-3.267,-2.191]
for i in range(0,6):
for w in range(i+1,6):
print round(math.sqrt((x[i]-x[w])**2 + (y[i]-y[w])**2 + (z[i]-z[w])**2),3),"light-years from ", a[i],"to", a[w]
Calculate the date of Easter, accurate from 1900 to 2368 inclusive.
Easter is celebrated on the first Sunday after the first full moon on or after 21 March.
#!/usr/bin/python3 import sys y=int(sys.argv[1]) c=y//100;n=y-19*(y//19);k=(c-17)//25;i=c-c//4-(c-k)//3+19*n+15 i=i-30*(i//30);i=i-(i//28)*(1-(i//28)*(29//(i+1))*((21-n)//11)) j=y+y//4+i+2-c+c//4;j=j-7*(j//7);l=i-j;m=3+(l+40)//44;d=l+28-31*(m//4) print (m,d,y,sep="/")
Print a range of Unicode characters
#!/usr/bin/python
import sys
a = int(sys.argv[1])
b = int(sys.argv[2])
for i in range(a,b + 1):
print(" ".join([str(i)," ",unichr(i)]))
Spell out numbers
pip install inflect
import inflect inflector=inflect.engine() print(inflector.number_to_words(8675309))
Find the point of intersection of two co-planar lines defined by four points
import numpy as np A = np.array([[4,0],[4,-3]]) B = np.array([[6,2],[10,2]]) t,s = np.linalg.solve(np.array([A[1]-A[0],B[0]-B[1]]).T, B[0]-A[0]) print((1-t)*A[0]+t*A[1])
Calculate the gamma function over a range
#!/usr/bin/python
import sys
import math
from numpy import arange
start = float(sys.argv[1])
stop = float(sys.argv[2])
step = float(sys.argv[3])
for x in arange(start, stop, step):
print (round(x,2),round(math.gamma(x),8))
Terminal Cruise 