Запуск Python кода под Windows

import math
from  decimal import Decimal
 
import matplotlib.pyplot as plt
 
def f (d):
    return Decimal(str(d))
 
 
 
def constTicks(array):
    array = array[1:]
    for i in range(1,len(array)):
        if abs(array[i-1] - array[i])>1:
            ma = int(round(max(array)))
            mi = int(round(min(array)))
            step = (ma-mi)//6
            ma += step
            return [int(i) for i in range(mi, ma+1, step)]
    x = int(round(array[1]))
    return [x-4, x, x+4]
 
 
 
def defirinV(coordinate, time, func):
    t = 0.000001
    coordinateDif = [[coordinate[0], eval(func)]]
    for i in range(1, len(coordinate)):
        t = time[i]+0.000001
        coordinateDif.append([coordinate[i], eval(func)])
    result = []
    for i in coordinateDif:
        result.append((i[1]-i[0])/0.000001)
    return result
 
def AAA(coordinate, time, func):
    t = time[0]
    coordinateDif = [[coordinate[0], eval(func)]]
    for i in range(1, len(coordinate)):
        t = time[i]+0.000001
        coordinateDif.append([coordinate[i], eval(func)])
    result = []
    for i in coordinateDif:
        result.append((i[1]-i[0])/0.000001)
    return result
 
def Crutch(func):
    last = func.find("f(t**2)")
    if last < 0: return 0
    first = last-3
 
    answer = 0
    k = 0
    while func[first] in [str(i) for i in range(10)]:
        answer = 10**k*func[first]+answer
        k += 1
        first -= 1
 
    if first -2 >= 0:
        if func[first-2] == "-": return -answer*2
 
    return answer*2
 
 
def defirinA(coordinate, time, func):
    speed1 = AAA(coordinate, time, func)
    speed2 = AAA([eval(func) for t in [i+0.000005 for i in time]], [i+0.000005 for i in time], func)
    print(speed1)
    print(speed2)
    result = []
    for i in range(len(speed1)):
        result.append((speed2[i]-speed1[i])/0.00001*2)
    result[0] = Crutch(func)
    return result
 
def parser(inp):
    row = ["cos", "sin","log","ln", "e^","^"]
    pretty = ["math.cos","math.sin","math.log", "math.log", "math.exp","**"]
    for i in range(len(row)):
        if row[i] in inp:
            print(i)
            inp = inp.replace(row[i], pretty[i])
 
    return inp
 
 
print("Введите уравнение S(t): ")
func = input()
print("Укажите кол-во секунд(целые числа):")
t1 = int(input())
func = parser(func)
print(func)
n =256*t1+1
a0 = 0
 
time = [i/128 for i in range(0,n) if ("log" not in func and "ln" not in func) or (i!=0)]
coordinate = [eval(func) for t in time]
movement = [i-coordinate[0] for i in coordinate]
speed = defirinV(coordinate, time, func)
axilir = defirinA(coordinate, time, func)
k = 0
if "log" in func or "ln" in func:
    k = 1
 
#print ("\n\n\n\ncoordinate  ", [str(i).ljust(20," ") for i in coordinate])
#print ("movement    ", [str(i).ljust(20," ") for i in movement])
#print ("time        ", [str(i).ljust(20," ") for i in time])
#print ("speed       ", [str(i).ljust(20," ") for i in speed])
#print ("axilir      ", [str(i).ljust(20," ") for i in axilir])
 
# sub 1
plt.subplot(221)
 
plt.plot(time[1+k:], coordinate[1+k:])
 
plt.grid(True)
plt.title(r'$Координаты\ от\ времени$')
plt.ylabel(u'Координата [м]', fontsize=12)
 
# sub 2
plt.subplot(222)
 
plt.plot(time[1+k:], movement[1+k:])
 
plt.grid(True)
plt.title(r'$Перемещения\ от\ времени$')
plt.ylabel(u'Перемещение [м]', fontsize=12)
 
 
# sub 3
plt.subplot(223)
plt.plot(time[k:], speed[k:])
 
plt.grid(True)
plt.title(r'$Скорости\ от\ времени$')
plt.xlabel(u'Время [c]', fontsize=12)
plt.ylabel(u'Скорости [м/c]', fontsize=12)
plt.yticks(constTicks(speed))
 
# sub 4
plt.subplot(224)
 
plt.plot(time[k+1:], axilir[k+1:])
 
plt.grid(True)
plt.title(r'$Полного\ ускорения\ от\ времени$')
plt.xlabel(u'Время [c]', fontsize=12)
plt.ylabel(u'Полного ускорения [м/c^2]', fontsize=12)
plt.yticks(constTicks(axilir))
 
 
plt.show()
 
#fout = open("calc.txt", "w")
#fout.write("N".ljust(5, " ")+"".join([i.ljust(15, " ") for i in ["coordinate", "movement", "time", "speed", "axilir"]])+"\n")
#for i in range(len(coordinate)):
#    fout.write(str(i+1).ljust(5, " ")+"".join([str(round(i, 4)).ljust(15, " ") for i in [coordinate[i], movement[i], time[i], speed[i], axilir[i]]])+"\n")
#fout.close()