,
        point_size    = 3,
        points(makelist([k, res[2][k]], k, p)),
        xlabel = &quot;Variances&quot;,
        ylabel = &quot;Percentages&quot;,
        xtics  = setify(makelist([concat(&quot;PC&quot;,k),k], k, p))) $
</pre>
<p>En caso de que la matriz de covarianzas sea conocida, se le puede pasar
a la función, pero debe utilizarse la opción <code>data=false</code>.
</p>
<pre class="example">(%i1) load (descriptive)$
(%i2) S: matrix([1,-2,0],[-2,5,0],[0,0,2]);
                              