Текст програми.
Program mzkt;
uses crt;
label 1, 2;
const E=10; tau1=4; tau2=1; T=6.0; n=20; Ut0=0;
var fajl : Text;
omega, dt, te : real;
C, B, A, fi, Ut, f : array[1..n] of real;
z1, z2, z3, Wser, C0 : real;
i, k : integer;
Begin
clrscr;
Assign(fajl, 'D:\Docs\mzramon.txt');
rewrite(fajl);
omega:=2*pi/T;
for i:=1 to n do f[i]:=i/T;
for k:=1 to n do
begin
z1:=(1.0/(k*omega*tau1))*(cos(k*omega*tau1-1));
z2:=((sqr(tau1)-tau1+tau2)/tau2)*sin(k*omega*tau1);
z3:=(1.0/(k*omega*tau2))*(cos(k*omega*tau1)-cos(k*omega*(tau1+tau2)));
C[k]:=(2*E/(k*omega*T))*(z1+z2+z3)
end;
for k:=1 to n do
begin
z1:=(1.0/tau1)*sin(k*omega*tau1);
z2:=(1.0/tau2)*(sin(k*omega*(tau1+tau2))-sin(k*omega*tau1));
B[k]:=((-1.0)*2*E)/(sqr(k*omega)*T)*(z1+z2)
end;
for k:=1 to n do
begin
A[k]:=sqrt(sqr(C[k])+sqr(B[k]));
fi[k]:=arctan(B[k]/C[k])
end;
Wser:=0;
for k:=1 to n do Wser:=Wser+sqr(A[k]);
Wser:=1.0/2.0*Wser;
for i:=1 to n do writeln(fajl, i, 'f=',f[i]);
writeln(fajl,'');
for k:=1 to n do writeln(fajl,'C',k,'=',C[k]);
writeln(fajl,'');
for k:=1 to n do writeln(fajl,'B',k,'=',B[k]);
writeln(fajl,'');
for k:=1 to n do writeln(fajl, 'A',k,'=',A[k]);
writeln(fajl,'');
for k:=1 to n do writeln(fajl, 'fi',k,'=',fi[k]);
writeln(fajl,'');
writeln(fajl,' Wser=',Wser);
dt:=1.0*T/n;
te:=0;
k:=1;
C0:=(2*E/T)*(3.0/2.0*tau1+tau2-tau1/2.0*(2.0*tau1+tau2));
C0:=C0/2.0;
1: te:=te+dt;
if te>T then goto 2;
Ut[k]:=0;
Ut[k]:=Ut[k]+C0;
for i:=1 to n do Ut[k]:=Ut[k]+C[i]*cos(i*omega*te)+B[i]*sin(i*omega*te);
k:=k+1;
goto 1;
2: writeln(fajl,'');
te:=0;
for i:=1 to n do
begin
te:=te+dt;
writeln(fajl,'U(t=',te:4:2,')= ',Ut[i])
end;
writeln(fajl,'');
writeln(fajl,' END');
close(fajl)
end.
Результати роботи програми.
1f= 1.6666666667E-01
2f= 3.3333333333E-01
3f= 5.0000000000E-01
4f= 6.6666666667E-01
5f= 8.3333333333E-01
6f= 1.0000000000E+00
7f= 1.1666666667E+00
8f= 1.3333333333E+00
9f= 1.5000000000E+00
10f= 1.6666666667E+00
11f= 1.8333333333E+00
12f= 2.0000000000E+00
13f= 2.1666666667E+00
14f= 2.3333333333E+00
15f= 2.5000000000E+00
16f= 2.6666666667E+00
17f= 2.8333333333E+00
18f= 3.0000000000E+00
19f= 3.1666666667E+00
20f= 3.3333333333E+00
C1=-3.9635076356E+01
C2= 1.8005309632E+01
C3= 7.2109461556E-01
C4=-9.0065359660E+00
C5= 7.0596295051E+00
C6= 1.1405014723E-02
C7=-5.1970070526E+00
C8= 4.4849923084E+00
C9= 8.0121624336E-02
C10=-3.5912284467E+00
C11= 3.2356116267E+00
C12= 2.8512539174E-03
C13=-2.7791219812E+00
C14= 2.5615192701E+00
C15= 2.8843784862E-02
C16=-2.2427387263E+00
C17= 2.0987102976E+00
C18= 1.2672242436E-03
C19=-1.8966478782E+00
C20= 1.7926901154E+00
B1= 6.5810039230E-01
B2= 1.1516756866E+00
B3= 1.3162109476E-11
B4=-2.8791892164E-01
B5=-2.6324015695E-02
B6=-9.6500385002E-13
B7= 1.3430620250E-02
B8= 7.1979730410E-02
B9= 7.7611080869E-12
B10=-4.6067027462E-02
B11=-5.4388462196E-03
B12=-4.8250192501E-13
B13= 3.8940851595E-03
B14= 2.3503585440E-02
B15= 3.3837658213E-12
B16=-1.7994932603E-02
B17=-2.2771639890E-03
B18=-1.2866718000E-12
B19= 1.8229927761E-03
B20= 1.1516756866E-02
A1= 3.9640539526E+01
A2= 1.8042104419E+01
A3= 7.2109461556E-01
A4= 9.0111368546E+00
A5= 7.0596785835E+00
A6= 1.1405014723E-02
A7= 5.1970244069E+00
A8= 4.4855698732E+00
A9= 8.0121624336E-02
A10= 3.5915239005E+00
A11= 3.2356161978E+00
A12= 2.8512539174E-03
A13= 2.7791247093E+00
A14= 2.5616270981E+00
A15= 2.8843784862E-02
A16= 2.2428109175E+00
A17= 2.0987115330E+00
A18= 1.2672242436E-03
A19= 1.8966487543E+00
A20= 1.7927271085E+00
fi1=-1.6602463904E-02
fi2= 6.3876097941E-02
fi3= 1.8252957645E-11
fi4= 3.1956892675E-02
fi5=-3.7287925201E-03
fi6=-8.4612240620E-11
fi7=-2.5842932702E-03
fi8= 1.6047642103E-02
fi9= 9.6866584410E-11
fi10= 1.2826948115E-02
fi11=-1.6809313739E-03
fi12=-1.6922446719E-10
fi13=-1.4011916848E-03
fi14= 9.1753851472E-03
fi15= 1.1731351615E-10
fi16= 8.0234698062E-03
fi17=-1.0850297433E-03
fi18=-1.0153465785E-09
fi19=-9.6116534637E-04
fi20= 6.4241992150E-03
Wser= 1.0647429924E+03
U(t=0.30)= -4.1962704665E+01
U(t=0.60)= -4.1316207392E+01
U(t=0.90)= -4.0671204475E+01
U(t=1.20)= -4.0068503037E+01
U(t=1.50)= -3.9666768512E+01
U(t=1.80)= -4.0707754141E+01
U(t=2.10)= 2.9912754693E+01
U(t=2.40)= 2.9744708758E+01
U(t=2.70)= 2.9770486377E+01
U(t=3.00)= 2.9970040161E+01
U(t=3.30)= 3.0275417113E+01
U(t=3.60)= 3.0757429487E+01
U(t=3.90)= 3.1454601734E+01
U(t=4.20)= -4.0216278357E+01
U(t=4.50)= -4.0913925575E+01
U(t=4.80)= -4.3052983615E+01
U(t=5.10)= -4.4474308869E+01
U(t=5.40)= -4.3831525046E+01
U(t=5.70)= -4.3218495534E+01
U(t=6.00)= 0.0000000000E+00
END

Графіки.