NameError: name 'S' is not defined. Did you mean: 'Sa'?

#! /usr/bin/env python3
import math
 
def FloatInput(message):
    while True:
        try:
            a=float(input(message).replace(',', '.'))
            return a
        except:
            print("Ошибка ввода вещественного числа! Попробуйте заново")
 
print("______________________ПАРАМЕТРЫ СНАРЯДА__________________________")
 
print("______________________Введите полукалибр,м_______________________")
D=FloatInput("D=")
print("____________________Введите высоту головной части,м______________")
H=FloatInput("H=")
print("_____________________Введите радиус оживала,м____________________")
R=FloatInput("R=")
print("_____________________Введите скорость при контакте,м/с___________")
V=FloatInput("V=")
print("_____________________Введите массу снаряда,кг____________________")
m=FloatInput("m=")
 
print("______________________ПАРАМЕТРЫ ПРЕГРАДЫ_________________________")
 
print("_________________Введите плотность преграды,кг/м**3______________")
ro=FloatInput("Ro=")
print("___________________Введите предел текучести,Па___________________")
sigmat=FloatInput("sigmat=")
print("________________Введите модуль упругости преграды,Па_____________")
E=FloatInput("E=")
print("В______________ведите кожффмцмент разупрочнения преграды_________")
limbda=FloatInput("limbda=")
print("_______________Введите коэффициент сопростивления преграды_______")
mi=FloatInput("mi=")
 
print("______________________ВВОД ГЕОМЕТРИИ_____________________________")
print("введите радиус головной части в соответствии с сечениием ")
r=str (input("r=").replace(',', '.'))
rr=[float(i) for i in r.split()]
 
r_sr = []
for i in range(1, len(rr)):
    print(f"r[i-1]={rr[i-1]} r[i]={r[i]}")
    r_sr.append((rr[i-1] + rr[i])*0.5)
    
print(rr)
print(f"{r_sr=}")
 
S_j = []
for i in range(1, len(rr)):
    print(f"r[i-1]={rr[i-1]} r[i]={r[i]}")
    S_j.append((rr[i-1]**2 - rr[i])**2)
    
print(f"{S_j=}")
 
print("______________введите координату сечения по оси x________________")
x=str (input("x=").replace(',', '.'))
xx=[float(i) for i in x.split()]
print(xx)
 
print("Введите угол наклона касательной к оживалу в соответствующем сечении")
ggr = input("grad=").replace(',', '.')
grr = [float(i) for i in ggr.split()]
 
grad = list(map(math.radians, grr))
 
print(grr)
print(f"grad = {grad}")
 
print("___________РАСЧЕТ КОЭФФИЦИЕНТОВ ОПИСЫВАЮЩИЕ СВОЙСТВА ПРЕГРАДЫ____")
 
print("___________________Прочностная составляющая______________________")
 
alfa=sigmat/(2*E*limbda)
print("Alfa=",sigmat/(2*E*limbda))
 
gamma=(1/(1-math.e**(-3*alfa)))**(1/3)
print("gamma=",(1/(1-math.e**(-3*alfa)))**(1/3))
 
A=2*(sigmat*math.log(gamma)+(2/3)*E*(((math.pi**2)/18)*(1-limbda)*((1-gamma**-3)**0.62)+alfa))
print("A=",(2*(sigmat*math.log(gamma)+(2/3)*E*(((math.pi**2)/18)*(1-limbda)*((1-gamma**-3)**0.62)+alfa))))
 
print("_________________Динамическая стационарная состовляющая__________")
 
С=ro/2
print("ro=",ro/2)
 
print("_________________Динамическая нестационарная состовляющая________")
D_j = []
for i in range(1, len(r_sr)):
    print(f"r_sr[i-1]={r_sr[i-1]} r_sr[i]={r_sr[i]}")
    D_j.append((ro*r_sr[i]*((gamma-1)/(gamma+1))))
    
print(f"{D_j=}")
 
print("_________________Учет промежуточных сечений________")
 
def Sa(S, x, h, k):
    if x[k] <= h:
        return S[k]
    else:
        if x[k]>h or x[k-1]<h:
            return S[k]*(h - x[k-1])/(x[k] - x[k-1])
        else:
            return 0
print(f'{S}{x}')

def Sa(S, x, h, k):
    if x[k] <= h:
        return S[k]
    else:
        if x[k]>h or x[k-1]<h:
            return S[k]*(h - x[k-1])/(x[k] - x[k-1])
        else:
            return 0
    print(f'{S}{x}')