Issue in calculating NDVI using rasterio python

Question

I am trying to calculate NDVI using two raster images (NIR & Red Bands) using the following code:

import rasterio as rio
import numpy as np
import matplotlib.