Tableau_type_symbol

import sympy as sp
import numpy as np
from matplotlib import pyplot as plt

def w(t,i,j,x):
    if t[i]==t[i+j]:
        return 0
    else:
        return (x-t[i])/(t[i+j]-t[i])
    
    
def inter_vides(t):
    
    ind_vides=[]
    for i in range(len(t)-1):
        if t[i]==t[i+1]:
            ind_vides.append(i) #indices des intervalles vides 
    return ind_vides
        

def genere_Bik(t,K):
    
        
   
    n = len(t)-1
    
    ind_vides=inter_vides(t)

        
        
    c=n-len(ind_vides)
    
    B=np.zeros((n,c))
    
    #initialisation de B avec les Bi0
    
    for i in range(n):
        a=0
        if not(i in ind_vides):
            B[i][a]==1
        a+=1
            
     #generation des Bik       
       
    x = sp.symbols('x') 
    
    for k in range (1,K,1):
        for i in range (n-k):
            B[i]=w(t,i,k,x)*B[i]+w(t,i+1,k,x)*B[i+1]
            
    return B