Мне нужно переделать этот код под мой вариант:
```
restart:
with(inttrans):
Rn_:=1700: Um:=140: T:=130/1000: f1:=1/(2*T): f2:=5/(2*T): n:=15: N:=500:
sigma:=t->piecewise(t>=0,1):
Uvh_1 := t-> Um/T * t * sigma(t):
Uvh_2 := t-> -Um * sigma(t-T):
Uvh:=t->Uvh_1(t)+Uvh_2(t):
plot(Uvh_1(t),t=-T..2*T);

plot(Uvh_2(t),t=-T..2*T);

plot(Uvh_1(t)+Uvh_2(t),t=0..T);

omega:=array(0..n):
 for k from 0 to n do
  omega[k]:=k*2*Pi/T od:
A:=array(0..n):
 for k from 0 to n do
  A[k]:=2/T*evalf(int(Uvh(t)*cos(omega[k]*t),t=0..T)) od;
                      A[0] := 140.0000000

                           A[1] := 0.

                           A[2] := 0.

                           A[3] := 0.

                           A[4] := 0.

                           A[5] := 0.

                           A[6] := 0.

                           A[7] := 0.

                           A[8] := 0.

                           A[9] := 0.

                          A[10] := 0.

                          A[11] := 0.

                          A[12] := 0.

                          A[13] := 0.

                          A[14] := 0.

                          A[15] := 0.

B:=array(1..n):
 for k to n do
  B[k]:=2/T*evalf(int(Uvh(t)*sin(omega[k]*t),t=0..T)) od;
                      B[1] := -44.56338406

                      B[2] := -22.28169204

                      B[3] := -14.85446135

                      B[4] := -11.14084601

                      B[5] := -8.912676810

                      B[6] := -7.427230676

                      B[7] := -6.366197722

                      B[8] := -5.570423006

                      B[9] := -4.951487116

                     B[10] := -4.456338406

                     B[11] := -4.051216732

                     B[12] := -3.713615336

                     B[13] := -3.427952620

                     B[14] := -3.183098862

                     B[15] := -2.970892270

Uvh_f:=t-> A[0]/2+sum('A[k]*cos(omega[k]*t)+B[k]*sin(omega[k]*t)',' k'=1..n);
Uvh_f := proc (t) options operator, arrow; (1/2)*A[0]+sum('A[k]*\

  cos(omega[k]*t)+B[k]*sin(omega[k]*t)', 'k' = 1 .. n) end proc


simplify(Uvh_f(t)); 
              70. - 44.56338406 sin(48.33219466 t)

                 - 22.28169204 sin(96.66438936 t)

                 - 14.85446135 sin(144.9965840 t)

                 - 11.14084601 sin(193.3287787 t)

                 - 8.912676810 sin(241.6609734 t)

                 - 7.427230676 sin(289.9931681 t)

                 - 6.366197722 sin(338.3253628 t)

                 - 5.570423006 sin(386.6575575 t)

                 - 4.951487116 sin(434.9897522 t)

                 - 4.456338406 sin(483.3219466 t)

                 - 4.051216732 sin(531.6541413 t)

                 - 3.713615336 sin(579.9863361 t)

                 - 3.427952620 sin(628.3185308 t)

                 - 3.183098862 sin(676.6507255 t)

                 - 2.970892270 sin(724.9829203 t)


plot({Uvh(t),Uvh_f(t)},t=0..2*T); 

Fchh:=proc(omega_) 
global P,Q,Res: 
local F: 
if P(omega_)>0 then F:=arctan(Q(omega_)/P(omega_)) end if: 
if P(omega_)<0 and Q(omega_)>=0 then 
F:=evalf(arctan(Q(omega_)/P(omega_))+Pi) end if: 
if P(omega_)<0 and Q(omega_)<0 then 
F:=evalf(arctan(Q(omega_)/P(omega_))-Pi) end if: 
if P(omega_)=0 and Q(omega_)>0 then F:=evalf(Pi/2) end if: 
if P(omega_)=0 and Q(omega_)<0 then F:=evalf(-Pi/2) end if: 
Res:=F: 
end proc:
omega1:=evalf(2*Pi*f1);
omega2:=evalf(2*Pi*f2);
                     omega1 := 24.16609734

                     omega2 := 120.8304867

omega0:=sqrt(omega1*omega2);
                     omega0 := 54.03703640

L1_:=evalf(Rn_/(2*Pi*(f2-f1)));
                       L1_ := 17.58662121

L2_:=evalf(Rn_*(f2-f1)/(2*Pi*f1*f2));
                       L2_ := 56.27718786

C1_:=evalf((f2-f1)/(2*Pi*f1*f2*Rn_)); 
                    C1_ := 0.00001947307539

C2_:=evalf(1/(2*Pi*(f2-f1)*Rn_));
                    C2_ := 0.000006085336058

ZL1:=L1*s: 
ZC1:=1/(C1*s):
ZL2:=L2*s:
ZC2:=1/(C2*s):
Z2:=ZL2*ZC2/(ZL2+ZC2):
Z1:=ZL1+ZC1:
Ku:=simplify(Z2*Rn/(Z1*Z2+(Z1+Z2)*Rn)):
Ku_n:=collect(numer(Ku),s):
Ku_d:=collect(denom(Ku),s):
Ku:=Ku_n/Ku_d;
           /       2   \//       4                3      
     Ku := \Rn L2 s  C1/ \Rn L1 s  C1 L2 C2 + L2 s  L1 C1

                                            2            \
        + (Rn L2 C2 + Rn L2 C1 + Rn L1 C1) s  + L2 s + Rn/


L1:=L1_: L2:=L2_: C1:=C1_: C2:=C2_: Rn:=Rn_:
Ku;
    /             2\//                 4                  3
    \1.863012867 s / \0.0001993805689 s  + 0.01927300095 s 

                      2                       \
       + 3.027395909 s  + 56.27718786 s + 1700/


s:=I*omega_;
                         s := I omega_

Ku_c:=evalc(Ku):
P:=coeff(Ku_c,I,0): 
P:=unapply(P,omega_):
Q:=coeff(Ku_c,I,1):
Q:=unapply(Q,omega_):
Achh:=sqrt(P(omega_)^2+Q(omega_)^2):
plot(Achh,omega_=0..N);

h:=1: 
omega_:=array(1..N): 
Fh:=array(1..N): 
for i from 1 to N do 
  omega_[i]:=h*i: 
  Fh[i]:=Fchh(omega_[i]):  
end:
Fhs:=[seq([omega_[i],Fh[i]],i=1..N)]:
PLOT(CURVES(Fhs,THICKNESS(2),LINESTYLE(1),COLOR(RGB,1.0,0.0,0.0)),AXESSTYLE(BOX)); 

Uvih_m_c:=array(0..n): 
Uvih_m_s:=array(1..n): 
phi:=array(1..n): 
omega_:=0: 
Uvih_m_c[0]:=A[0]*Achh; 
for k to n do 
  omega_:=evalf(omega[k]); 
  Uvih_m_c[k]:=evalf(A[k]*Achh); 
  Uvih_m_s[k]:=evalf(B[k]*Achh); 
  phi[k]:=evalf(Fchh(omega_)) od:
                       Uvih_m_c[0] := 0.

U_vih:=t->Uvih_m_c[0]/2+sum('Uvih_m_c[k]*cos(omega[k]*t+phi[k])+Uvih_m_s[k]*sin(omega[k]*t+phi[k])','k'=1..n);
U_vih := proc (t) options operator, arrow; (1/2)*Uvih_m_c[0]+sum\

  ('Uvih_m_c[k]*cos(omega[k]*t+phi[k])+Uvih_m_s[k]*sin(omega[k]*\

  t+phi[k])', 'k' = 1 .. n) end proc


 simplify(U_vih(t));
      -44.91010038 sin(48.33219466 t + 0.1263081175)

         - 25.71586681 sin(96.66438936 t - 0.9164750817)

         - 10.21373950 sin(144.9965840 t - 2.048324518)

         - 3.681722104 sin(193.3287787 t - 2.486408627)

         - 1.709681406 sin(241.6609734 t - 2.668561942)

         - 0.9362140006 sin(289.9931681 t - 2.767933745)

         - 0.5701104136 sin(338.3253628 t - 2.831193447)

         - 0.3736963395 sin(386.6575575 t - 2.875356406)

         - 0.2585711893 sin(434.9897522 t - 2.908109394)

         - 0.1864997706 sin(483.3219466 t - 2.933456541)

         - 0.1390207231 sin(531.6541413 t - 2.953701978)

         - 0.1064425037 sin(579.9863361 t - 2.970272049)

         - 0.08333131370 sin(628.3185308 t - 2.984100310)

         - 0.06647406288 sin(676.6507255 t - 2.995825184)

         - 0.05388539111 sin(724.9829203 t - 3.005899010)


plot(U_vih(t),t=0..2*T);
```
Исходные данные моего варианта:
```
Тип фильтра: Г-образный ПФ типа «k» с Т-образным входом
Граничные частоты: $f_1=\frac{1}{2T};f_2=\frac{5}{2T}$
Сигнал: смотри описание ниже
$R_н$=500 Ом
Параметры сигнала:
$U_m$=36 В
T=10 мс
```
Задания выполняются в Maple версии 6.
Описание сигнала:
```md
Это периодический кусочно-линейный сигнал с периодом $T$, амплитудой $\pm U_m$ и нулевым средним, который можно описать так:

1. **Форма**

   * На интервале $0 \le t < \tfrac T2$ идёт линейный спад от нуля до $-U_m$.
   * В точке $t=\tfrac T2$ мгновенный скачок вверх с $-U_m$ до $+U_m$.
   * На интервале $\tfrac T2 < t < T$ — линейное убывание от $+U_m$ до нуля.

2. **Математически**

$$
u(t)=
\begin{cases}
-\dfrac{2U_m}{T}\,t, & 0 \le t < \tfrac T2,\\[1em]
U_m - \dfrac{2U_m}{T}\,(t - \tfrac T2), & \tfrac T2 \le t < T,\\[0.3em]
u(t+T)=u(t).
\end{cases}
$$

3. **Ключевые особенности**

* **Разрывы:** скачок напряжения в середине периода.
* **Линейность участков:** каждая половина — обычная «треугольная» линейная функция.
* **Среднее значение:** ноль (положительная и отрицательная области равны по площади).

Такой сигнал иногда называют **асимметричным пилообразным/треугольным** с центральным разрывом.
```