from numpy import *
from matplotlib.pyplot import *
pi=3.14159265359

Tmin = -0.5
Tmax = 0.5

def f(t):
    if(t**2<0.25):
        return 1
    else:
        return 0
def g(t):
    if(t**2<0.25):
        return 1
    else:
        return 0

Int='(0'
i=0
dt=0.01
for tau in linspace(Tmin,Tmax,100):
    #Int=Int+'+(((f(t-'+str(tau)+')*g('+str(tau)+')+f(t-'+str(tau)+'+0.1)*g('+str(tau)+'+0.1))/2)*(0.1)/(Tmax-Tmin))'
    i=i+1
    Int=Int+'(Tmin+'+str(i)+'*dt)*(f(t)*g(t-'+str(tau)+')'
#Int=Int+')'

lt=[]
lv=[]
for t in linspace(-2,2,100):
    lt.append(t)
    lv.append(eval(Int))
plot(lt,lv)
grid()
show()