C5: inverse transform

Here we demonstrate the inverse transform defined in (T2).

%%% prelim

x=linspace(1,100,25); %equally spaced values bet. 1 and 100

xprime=1./(x); %transformed coordinates

p=ones([length(xprime) 1])./length(xprime); %uniform distribution

%%% plotting

figure; subplot(2,1,1); %create figure and sub-figure for original plot

stem(x,p,'k-','MarkerFaceColor','k','MarkerSize',11); %plot black (colorcode ‘k’) stems

axis([-1 102 0 1.05*p(1)]); box off

subplot(2,1,2); stem(xprime,p,'k-','MarkerFaceColor',.6*[1 1 1],'MarkerSize',11);

axis([-.01 1.02 0 1.05*p(1)]); box off