Generating random variable from inverse gaussian dist.

Hi. I want to generate random variable from inverse gaussian dist with accepted-rejection method but I didnt do it. help me please..

Answer 1

[Jeff Miller]

Do you really need to write the random generator yourself? This one is already available:

>> pd = makedist('inversegaussian')
%pd = 
%  InverseGaussianDistribution
%
%  Inverse Gaussian distribution
%        mu = 1
%    lambda = 1
>> pd.random(1,10)
% ans =
%      0.82544      0.35561      0.46273       6.1449       2.5297      0.22435      0.20619        2.594       0.2875      0.89231