Ring object in STL file

import math

out = "ring.stl"

f = open(out, 'w')

n = 10
m = 20

r1 = 5.0
r2 = 1.0

dtheta = 2.0 * math.pi / m
dphi = 1.0 * math.pi / n

f.write("solid Something (Meshed)\n")
for i in range(0, m):
    theta = 2.0 * math.pi * i * 1.0 / m
    ct = math.cos(theta)
    st = math.sin(theta)
    ct1 = math.cos(theta + dtheta)
    st1 = math.sin(theta + dtheta)
    for j in range(0, n):
        phi = math.pi * j * 2.0 / n
        cp = math.cos(phi)
        sp = math.sin(phi)
        cp1 = math.cos(phi + dphi)
        sp1 = math.sin(phi + dphi)
        f.write("\tfacet normal %.6f %.6f %.6f\n" % (ct * cp, st * cp, sp))
        f.write("\t\touter loop\n")
        f.write("\t\t\tvertex %.6f %.6f %.6f\n" % ((r1 + r2*cp) * ct, (r1 + r2*cp) * st, r2*sp))
        f.write("\t\t\tvertex %.6f %.6f %.6f\n" % ((r1 + r2*cp1) * ct, (r1 + r2*cp1) * st, r2*sp1))
        f.write("\t\t\tvertex %.6f %.6f %.6f\n" % ((r1 + r2*cp) * ct1, (r1 + r2*cp) * st1, r2*sp))
        f.write("\t\tendloop\n")
        f.write("\tendfacet\n")
        f.write("\tfacet normal %.6f %.6f %.6f\n" % (ct * cp, st * cp, sp))
        f.write("\t\touter loop\n")
        f.write("\t\t\tvertex %.6f %.6f %.6f\n" % ((r1 + r2*cp1) * ct1, (r1 + r2*cp1) * st1, r2*sp1))
        f.write("\t\t\tvertex %.6f %.6f %.6f\n" % ((r1 + r2*cp1) * ct, (r1 + r2*cp1) * st, r2*sp1))
        f.write("\t\t\tvertex %.6f %.6f %.6f\n" % ((r1 + r2*cp) * ct1, (r1 + r2*cp) * st1, r2*sp))
        f.write("\t\tendloop\n")
        f.write("\tendfacet\n")

f.write("endsolid Mesh\n")

f.close()