4806.sagews

def timeE(t):
    f(a,x,b)=cos((a)^2*pi*b/1281)*sin(a*x*pi)/a
    g(a,x,b)=sin((a)^2*pi*b/1281)*sin(a*x*pi)     
    h(x,b)=sum(f(2*z+1,x,b) for z in range(0,90))        
    p1=plot(h(x,t),-1,1,ymin=-4,ymax=4)       
    return p1    
 
a=animate(timeE(i) for i in range(0,25))
a.show()