Товарищи! Помогите написать программу с применением формулы корреляции в питоне.ПОЖАЛУЙСТА

import matplotlib.pyplot as plt
%matplotlib inline 
import numpy
import numpy as np
from numpy import ndarray
y1=ndarray((21,),int)
y2=ndarray((21,),int)
%Кroc
corf=ndarray(15,int)
 for j in range (0,15,1):
    sum(y1y2)=0
    sum(y2)=0
    for i in range (0,len(y1)-j,1):
        sum(y1y2)=sum(y1y2)+y1[i]*y2[j+i]
        sum(y2)=sum(y2)+y2[i+j]
    sum(y1y2)=sum(y1y2)/(len(y1-j))
    s2=0
    for i in range (0,len(y1)-j,1):
        s2=s2+sum(y2)*y1[i]
    s2=s2/((len(y1)-j)**2)
         corf[j]=(sum(y1*y2)-s2)/((cv1*cv2)**.5)
     print (corf[j])
 %поменять потом местами y1 и y2

Пишет ошибки не могу разобраться . нужно чтобы выводилось каждое значение кореляции(пишу на платформе анаконда в юпитере)