exponential-distribution

I have some points that look like a logarithmic curve. The curve that I'm trying to obtain look like: y = a * exp(-b*x) + c My code: x <- c(1.564379666,1.924250092,2.041559879,2.198696382,2.541267447,2.666400433,2.922534874,2.965726615,3.009969443,3.248480245,3.32927682,3.371404563,3.423759668,3.713001284,3.841419166,3.847632349,3.947993339,4.024541136,4.030779671,4.118849343,4.154008445,4.284232251,4.491359108,4.585182188,4.643299476,4.643299476,4.643299476,4.684369939,4.84424144,4.867973977,5.144490521,5.324298915,5.324298915,5.988637637,6.146599422,6.674937463) y <- c(25600,23800,11990

I would like to generate some pseudorandom numbers and up until now I've been very content with the .Net library's Random.Next(int min, int max) function. PRNGs of this variety are supposed to be using a Uniform distribution , but I would very much like to generate some numbers using an Exponential Distribution . I'm programming in C#, although I'll accept pseudocode or C++, Java or the like. Any suggestions / code snippets / algorithms / thoughts? Since you have access to a uniform random number generator, generating a random number distributed with other distribution whose CDF you know is