To complete the comparison of the techniques, add the lines
[Wcg1, Ecg1]=pccg(W,50,1); %# 50 iters from W of cg versions 1, 2, 3 ... [Wcg2, Ecg2]=pccg(W,50,2); [Wcg3, Ecg3]=pccg(W,50,3);immediately after the call to pclm(W,50) in your current version of ``sch.m'' and replace the current plotting call to semilogy(...) with
semilogy([1:length(Elm)],Elm-43.3371147782040,'-;lm;', ...
[1:length(Epclm)],Epclm-43.3371147782040,'-;pclm;', ...
[1:length(Ecg1)],Ecg1-43.3371147782040,'-;pccg-FR;', ...
[1:length(Ecg2)],Ecg2-43.3371147782040,'-;pccg-PR;', ...
[1:length(Ecg3)],Ecg3-43.3371147782040,'-;pccg-HS;');
%# Set W to best result for plotting
W=Wcg3;