Imports System.Math
Imports System.IO
Module Module1
 Sub Main()
  Dim a(14) As Long
  Dim dmean As Double = 0.0025
  Dim dose As Double = 0.01
  Dim freq(1001) As Double
  Dim g(1001) As Double
  Dim i, j As Long
  Dim mz(1001) As Double
  Dim Xso(1001) As Double
  Dim n As Double = 101
  Dim nz(1001) As Double
  Dim p As Double = 1300
  Dim pot(1000000) As Double
  Dim pmean As Double
  Dim potency As Double = 0
  Dim pt As Double = 0
  Dim r(1001) As Double
  Dim rinc As Double
  Dim rmax As Double
  Dim rmean As Double
  Dim rmin As Double
  Dim rs As Double
  Dim rsd As Double
  Dim rt As Double = 0
  Dim s As Double
  Dim sd As Double = 2
  Dim t As Long = 100000
  Dim tfreq As Double = 0
  Dim vz(1001) As Double
  For i = 1 To 14
   a(i) = 0
  Next i
  rmean = dmean / 2
  rmax = rmean * sd ^ 3
  rmin = rmean / sd ^ 3
  rinc = Log(rmax / rmin) / (n - 1)
  Dim file As StreamWriter
  file = My.Computer.FileSystem.OpenTextFileWriter(
   "C:\Users\Kevin C. Johnson\cu.txt", False)
  For i = 1 To n
   r(i) = Exp(rinc * (i - 1) + Log(rmin))
   freq(i) = 1 / Log(sd) / Sqrt(2 * PI) * Exp(-((Log(r(i)) - Log(rmean)) / Log(sd)) ^ 2 / 2)
   tfreq = tfreq + freq(i)
   vz(i) = 4 / 3 * PI * r(i) ^ 3
   mz(i) = vz(i) * p
   g(i) = 0
  Next i
  For i = 1 To n
   Xso(i) = freq(i) / tfreq * dose
   nz(i) = Xso(i) / mz(i)
  Next i