420.sagews

omg = 2*pi/5
def wave0(t): return -10*cos(pi*1)/(pi)*sin(1*omg*t)
    
def wave1(t): return -10*cos(pi*2)/(2*pi)*sin(2*omg*t)
    
def wave2(t): return -10*cos(pi*3)/(3*pi)*sin(3*omg*t)
    
def wavesum(t): return wave0(t) + wave1(t) + wave2(t)

plotRange = (-5/2,5/2)
plot0 = plot(wave0,plotRange,color='black')
plot1 = plot(wave1,plotRange,color='red')
plot2 = plot(wave2,plotRange,color='green')
plotsum = plot(wavesum,plotRange,color='blue')
show(plot0+plot1+plot2+plotsum)

def waveTerm(t, i): return -10*cos(pi*i)/(pi*i)*sin(i*omg*t)

def waveS(t): return sum((waveTerm(t,i) for i in xrange(1,40,1)))

plot(waveS, plotRange)

def waveS(t): return sum((waveTerm(t,i) for i in xrange(1,25,1)))

plot(waveS, plotRange)