def integrand(k): first = 6*(np.sin(k*r)-k*r*np.cos(k*r))*(((k*r)**2 -3)*np.sin(k*r)\ +3*k*r*np.cos(k*r))/(k*r)**7 *linpwer.deltasq(k,1.) second = 6*(np.sin(r/k)-r/k*np.cos(r/k))*(((r/k)**2 -3)*np.sin(r/k)\ +3*r/k*np.cos(r/k))/(r/k)**7 *linpwer.deltasq(1/k,1.) return first + second
def integrand(k): first = 9.*(np.sin(k*r)/(k*r)**3 -np.cos(k*r)/(k*r)**2)**2\ *linpwer.deltasq(k,1)/k second =9*(np.sin(r/k)/(r/k)**3 -np.cos(r/k)/(r/k)**2)**2\ *linpwer.deltasq(1/k,1)/k return math.sqrt(first + second)
#!/usr/bin/env python import linpwer import numpy as np import math import matplotlib.pyplot as plt #This program calls the deltasq module to calculate the linear power spectrum k=1000 x=np.linspace(.0001,k,10000) delta=[] for i in x[:]: delta.append(linpwer.deltasq(i,1)*10) fig=plt.figure() ax=fig.add_subplot(2,1,1) ax.plot(x,delta) ax.set_xscale('log') #plt.plot(x,delta) #plt.axes.set_xscale('log') plt.show()