Python function for computing rms integral from high frequency to low frequency

rms.py

import numpy as np

'''
Find the rms value of y(x). Here y(f) could be a spectral density or y(t) could be a time series

usage:
    rms_of_y = rms(X,Y);
'''


def rms(x, y):
    diff = np.diff(np.squeeze(x))
    dx = np.concatenate(([diff[0]], diff), axis=0)

    # Return rms intergrated from high to low
    rms = np.flipud(np.sqrt(np.cumsum(np.flipud(np.squeeze(y)**2 * dx))))
    return rms