NDVI Landsat 5,7 and 8

Hi everyone, I’m studying historical land cover change from 1984 to 2020, I need to calculate NDVI for every image, can anyone help me with a python script to do the task? considering all 3 sensors of Landsat.
Thank you