log[((Subscript[v(r), x])/(Subscript[v, o]))(x/D)]=0.79-33((r^2)/(x^2))