import random
import numpy as np

def probtrials(target, num_psywaves, num_trials):

    successful_trials = 0

    psywave_rolls = np.zeros(num_psywaves)

    for i in range(num_trials):
        for j in range(num_psywaves):
            psywave_rolls[j] = random.randint(50,150)

        average_roll = np.mean(psywave_rolls)

        if average_roll >= target:
            successful_trials += 1

    return successful_trials/num_trials