Help me resolve an equation using the Fourier series

This is my Fourier analisis about single phase rectifier resitive load
The voltage load and current load
(1)
There is a concept related to nonsinusoidal power definition
(2)
I got in matlab the next using the discrete fourier block wich gives me the ability to read the magnitud an angle phase of every harmonic order...

DC COMPONENT = 1.151467E+01

HARMONIC FREQUENCY FOURIER NORMALIZED PHASE NORMALIZED
NO (HZ) COMPONENT COMPONENT (DEG) PHASE (DEG)

1 5.000E+01 1.821E+01 1.000E+00 7.078E-03 0.000E+00
2 1.000E+02 7.965E+00 4.373E-01 -8.999E+01 -9.001E+01
3 1.500E+02 1.937E-01 1.064E-02 1.798E+02 1.797E+02
4 2.000E+02 1.583E+00 8.689E-02 -9.001E+01 -9.004E+01
5 2.500E+02 1.149E-01 6.310E-03 1.789E+02 1.788E+02
6 3.000E+02 6.708E-01 3.683E-02 -8.990E+01 -8.994E+01
7 3.500E+02 8.038E-02 4.413E-03 1.790E+02 1.790E+02
8 4.000E+02 3.670E-01 2.015E-02 -8.987E+01 -8.993E+01
9 4.500E+02 5.977E-02 3.281E-03 1.785E+02 1.784E+02
10 5.000E+02 2.313E-01 1.270E-02 -8.972E+01 -8.979E+01
11 5.500E+02 4.860E-02 2.668E-03 1.782E+02 1.781E+02
12 6.000E+02 1.576E-01 8.654E-03 -8.948E+01 -8.956E+01
13 6.500E+02 4.077E-02 2.238E-03 1.782E+02 1.781E+02
14 7.000E+02 1.132E-01 6.215E-03 -8.924E+01 -8.934E+01


TOTAL HARMONIC DISTORTION = 4.483567E+01 PERCENT

well indeed tis a report about orcad but the fourier component an phase angle is what I got in matlab simulink

My question is how do I make or resolve the equation (2) using the series of fourier or the report got from matlab simulink

thanks for your replies