function TiO2Rutile_PeakIntensity

Lambda = 1.54056;
u = [0.3; 1.0/3.0];
LC = [[4.5936,2.9587];[4.5845,2.9533]];

for i = 1:2
    for j = 1:2
        fprintf(1,'\n');
        fprintf(1,'Lambda=%f Angstroms a=%f c=%f Angstroms u=%f\n',Lambda,LC(i,1),LC(i,2),u(j));
        fprintf(1,'(h,k,l) M d(Ang)  T(Deg)  fTi    fO  nfTi   nfO  S(hkl) LP(T) Intensity\n');
        TiO2Rutile_PeakIntensity2(1,0,0,u(j),Lambda,LC(i,1),LC(i,2));
        TiO2Rutile_PeakIntensity2(1,1,0,u(j),Lambda,LC(i,1),LC(i,2));
        TiO2Rutile_PeakIntensity2(1,0,1,u(j),Lambda,LC(i,1),LC(i,2));
        TiO2Rutile_PeakIntensity2(2,0,0,u(j),Lambda,LC(i,1),LC(i,2));
        TiO2Rutile_PeakIntensity2(2,0,1,u(j),Lambda,LC(i,1),LC(i,2));
        TiO2Rutile_PeakIntensity2(1,1,1,u(j),Lambda,LC(i,1),LC(i,2));
        TiO2Rutile_PeakIntensity2(2,1,0,u(j),Lambda,LC(i,1),LC(i,2));

        TiO2Rutile_PeakIntensity2(1,1,0,u(j),Lambda,LC(i,1),LC(i,2));
        TiO2Rutile_PeakIntensity2(1,0,1,u(j),Lambda,LC(i,1),LC(i,2));
        TiO2Rutile_PeakIntensity2(1,1,1,u(j),Lambda,LC(i,1),LC(i,2));
        
        
        %TiO2Rutile_PeakIntensity2(1,2,0,u(j),Lambda,LC(i,1),LC(i,2));
        %TiO2Rutile_PeakIntensity2(1,1,1,u(j),Lambda,LC(i,1),LC(i,2));
        %TiO2Rutile_PeakIntensity2(0,2,0,u(j),Lambda,LC(i,1),LC(i,2));
        %TiO2Rutile_PeakIntensity2(0,1,1,u(j),Lambda,LC(i,1),LC(i,2));
        %TiO2Rutile_PeakIntensity2(1,1,0,u(j),Lambda,LC(i,1),LC(i,2));
        
    end
end