Система дифференциальных уравнений

Restart:
k = 1   # номер коэффициента ряда Фурье
RR = 100  # см    радиус сферической оболочки
h = 1  #  см   толщина сферической оболочки
EU = 21000   #  кН/см2  модуль упругости материала
ν = 0.5  #   коэффициент Пуассона
θ0 = 180   #  градусы
δ  = 3   #   градусы
θ1 = 90   #  градусы
qq = 0.1  #    кн/см2
qrk = 0.0  #    
 qzk=qq/EU    #   безразмерная величина
qφk = 0.0  #   
 
qzk = piecewise( 0,7 ≤ x ≤ 0.8) 
 
θ1/θ0≤x≤(1-δ/θ0)
θi = x∙θ0
ri = RR∙sin(θ0(1− x))
p0=(1+  (h^2∙〖sin⁡θi〗^2)/(3∙〖ri〗^2 ))
f1 = 
f2 = 
f3 = 
Y1(),Y2(
граничные условия 
Y5(s0) = 0, Y6(s0) = 0, Y7(s0) = 0, Y8(s0) =0,
Y5(sn) = 0, Y6(sn) = 0, Y7(sn) = 0, Y8(sn) =0,
si=RR∙(θi-θ1)
si/(RR∙θ0)=  RR/RR  ∙(θi/θ0-θ1/θ0)
 
si/(RR∙θ0 )=ss,      si=ss∙θ0∙RR,         ss=(θi/θ0- θ1/θ0)=(x-θ1/θ0)
ss=(x-θ1/θ0)
ss0=(θ1/θ0-θ1/θ0)=0                       ssn=(1-δ/θ0-θ1/θ0)
 
0≤ss≤(1-δ/θ0-θ1/θ0)                            dT/(dsi )=dT/(d(ss∙θ0∙RR))=1/(θ0∙RR)  ∙dT/dss
dY1/dss=θ0∙RR∙[-(1-ν)]∙(cos⁡θi  )/ri∙Y1(ss)+ν∙sin⁡θi/ri∙Y2(ss)-
-k∙(cos⁡θi  )/ri∙Y3(ss)+(ν∙k^2∙h∙sin⁡θi)/〖ri〗^2 ∙Y4(ss)+
+h/〖ri〗^2 ∙(1+(k^4∙h^2∙sin^2⁡θi)/(12∙〖ri〗^2 ))∙Y5(ss)-(k^4∙h^3∙sin⁡θi∙cos⁡θi)/(12∙〖ri〗^4 )∙Y6(ss)+
+(k∙h)/〖ri〗^2 ∙(1+(k^2∙h^2∙sin^2⁡θi)/(12∙〖ri〗^2 ))∙Y7(ss)+(k^2∙h^2∙sin⁡θi∙cos⁡θi)/(12∙〖ri〗^2 )∙Y8(ss)-qrk/h]
 
dY2/dss=θ0∙RR∙[…………..Y8(ss)-qzk/h]
dY3/dss=θ0∙RR∙[…………..Y8(ss)-qφk/h]
dY4/dss=θ0∙RR∙[…………..Y8(ss)]
dY5/dss=θ0∙RR∙[…………..Y8(ss)]
dY6/dss=θ0∙RR∙[…………..Y8(ss)]
dY7/dss=θ0∙RR∙[…………..Y8(ss)]
dY8/dss=θ0∙RR∙[…………..Y8(ss)]