# wave.jl using Winston # A_n = 32/(pi^3 n^3) if n is odd # wave parameters cc = 10.0 # x parameters dx = 0.1 LL = 5.0 trunc = 20 # t parameters period = 2*LL/cc nsteps = 10 dt = 0.5*period/nsteps # Plots cmap = Colors.colormap("RdBu",nsteps) p1 = plot() p2 = plot() # Math xx = collect(0:dx:LL) uu = 0.*xx tm = 0. for nn=1:nsteps uu = 0.*xx for kk=1:trunc ak = 32/(pi^3*kk^3) * isodd(kk) uu += ak*sin(pi*kk*xx/LL) * cos(pi*kk*cc*tm/LL) end if true p1 = oplot(xx,uu,color=cmap[nn]) else p2 = plot(xx,uu) xlim(0,LL); ylim(-1,1) figname = string("./figs/fig",string(nn),".pdf") savefig(p2, figname, width=600, height=400) println(nn) end tm += dt end figure(width=500,height=400); display(p1) xlabel("x"); ylabel("u(x,t)")