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