想做一个简单的拟合，条件：给定一组数据，它们服从对数正态分布，或者说，在区间[15，49]服从对数正态分布，现在想用对数正态分布的函数将其拟合出来，就是要估计出分布函数中的那两个参数。书中给出了方法，应该是要在MATLAB的M文件里进行编程，但是看了半天不懂，希望高手帮忙看一下：
step1：由已知数据得到一组散点图b(r)，这组散点图大概服从对数正态分布,
       由b(r)得到分布函数F(r)=F(r-1)+b(r);(注：r是从15到49的整数)
step2:令a=Max{F(r)},b = Min{F(r)}，data0 为1*10000 矩阵。在区间[b,a] 上产生    10000个服从均匀分布的数据,用矩阵data 记录。若F(r -1) < data(i) <F(r)， data0(i) = r 。那么， data0 包含的10000个数据服从对数正态分布。
step3 ：利用data0 的数据，调用matlab 统计工具箱中lognfit 函数，估计参数b(r)的参数u、d 。（这两个参数是对数正态分布的参数）
希望高手帮帮忙，给出这三个步骤下的matlab程序，谢谢！！