Hi guys,
I want to integrate the function q_dim function over z wich gives a vector as result
z=np.linspace(0,4,100)
q_dim=1+10*((1*62)**(-2/3)*(z/10.)**(-2/3))*1./2.*1.0**2*(z/10)**2*1
Q_dim_int==-integrate.cumtrapz(q_dim,z)
Q_dim_Analytically=................
The result always looks like this (Black is analytic solution) - What I am doing wrong?
[IMG]http://www.pixoload.de/image.php?dt=MZIR