import pyfits
import scipy, pylab

dcube = raw_input('Please select the datacube: ')

dcube = scipy.load('%(dcube)s'%vars())

i = input('Which row: ')
i = int(i)
                                          #getting the row and column
j = input('Which column: ')
j = int(j)

y1 = dcube[:,i,j]                             #extracts a vector row i column j
y1axis = scipy.linspace(0, len(y1)-1,num=len(y1))

##########Plotting
pylab.subplot(111)
pylab.plot(y1axis,y1)
pylab.title('Datacube, row = %(i)s, column = %(j)s'%vars())
pylab.ylabel('Number of counts')
pylab.show()