Integrals - Exemple 2 (2016)

Ejercicio Catalán
Universidad Universidad de Barcelona (UB)
Grado Física - 3º curso
Asignatura Eines Informàtiques
Año del apunte 2016
Páginas 1
Fecha de subida 19/03/2016
Descargas 1
Subido por

Vista previa del texto

Exercici 2: from __future__ import division from pylab import * close('All') interactive(True) from scipy import integrate #funcio a integrar def pendol(theta,theta_0): return 1/(sqrt((sin(theta_0/2.)**2)-(sin(theta/2.)**2))) N=10 #pasos theta0_array=linspace(0+0.001,pi/4.,N) #valors thtea_0 integral=zeros(N) for i, th0 in enumerate(theta0_array): #guarda els valors de l'integral en un vector integral[i],error=integrate.quad(pendol,0,th0,args=th0) #integral #grafic figure() plot(theta0_array,integral) xlabel(r'$\theta_o$') ylabel('Integral') title(r"Integral oscil.lacions pendol per a diferents $\theta_0$" ) ...