I'm trying to speed up the following code:

```
from math import log
from random import random

def logtest1(N):
    tr=0
    for i in range(1,N):
        T= 40 + 10*random()
        tr += -log(random())/T

```

I'm fairly new to python (coming from matlab)... and this same code runs 5x slower in python than matlab (and Julia), which got my attention.

I tried using a numba and a parakeet wrapper, and numpy functions instead of python ones, but didn't get any improvement at all.

I'd appreciate any help.
Thanks.

edit: the whole thing is a Monte Carlo simulation, so N is very large... **10e6** for testing purposes