obsline=csvread('obsline.csv');
nobl=length(obsline(:,1));

for i=1:nobl
    nfrom=obsline(i,1);
    nto=obsline(i,2);
    ni=locLN(nfrom, nto, line);
    mui=PLineX(ni,1);
    smgi=sqrt(PLineX(ni,2));
    xi=-150;
    for j=1:10
        PCL(i,j)=xi+(j-1)*30;
        PCDFL(i,j)=NCalCdf(PCL(i,j), mui, smgi, PLineU(ni, 3), PLineU(ni,4), PLineU(ni,5));
        PPDFL(i,j)=NCalPdf(PCL(i,j), mui, smgi, PLineU(ni, 3), PLineU(ni,4), PLineU(ni,5));
    end    
    figure(i);
    plot(PCL(i,:), PCDFL(i,:))
    figure(2*i);
    plot(PCL(i,:), PPDFL(i,:))
end   

PCL1=PCL(1,2:12);
for i=1:11
    PL(1,i)=PPL(1,i+1)-PPL(1,i);
end
figure(2)
plot(PCL1(1,:),PL(1,:));