Rafael Costa wrote: > Why Sage not execute the python code? > > from scipy import * > from pylab import * > > sample_rate=1000.00 > t=r_[0:0.6:1/sample_rate] > N=len(t) > s=sin(2*pi*50*t)+sin(2*pi*70*t+pi/4) > S=fft(s) > f=sample_rate*r_[0:(N/2)]/N > n=len(f) > plot(f,abs(S[0:n])/N) >
See http://sagenb.org/home/pub/676/ for a working version: from scipy import * from pylab import * sample_rate = 1000.0 t = r_[0:0.6:1/sample_rate].astype(float) N = len(t) s = sin(2*pi*50*t)+sin(2*pi*70*t+pi/4) S = fft(s) f=sample_rate*r_[0:(N/2)]/N n=len(f) line(zip(f,abs(S[0:n])/N)) As Minh points out, the need for the "astype" above will go away once the three tickets (that are up for review!) are reviewed and merged; hopefully in the next few days. Thanks, Jason --~--~---------~--~----~------------~-------~--~----~ To post to this group, send email to sage-devel@googlegroups.com To unsubscribe from this group, send email to sage-devel-unsubscr...@googlegroups.com For more options, visit this group at http://groups.google.com/group/sage-devel URLs: http://www.sagemath.org -~----------~----~----~----~------~----~------~--~---
