programacion

Ingeniería Computacional
Juan Carlos Gutiérrez Arredondo
Proyecto #6

1.
!INTEGRACION POR EL METODO MONTE CARLO

!definicion de variables
double precision sum, f, r1, prom
integer i, wcharacter name, ig

!Archivo
open(unit=18,file='RESULTADOS MONTECARLO')

sum=0
w=10
do i=1,10000
call random_number(r1)
f=1/(1+(r1)**2)
sum=sum+f
prom=sum/i
res=mod(i,w)
if (res.eq.0)then

write(18,*)i,prom,f
end if

end do
end

El valor de la integral según el método Monte Carlo es aproximadamente π/4,y se observa que el método converge aproximadamente a las 250iteraciones












2.
!definicion de variables
double precision sum, f, r1, prom, N
integer i,imp


!Archivo
open(unit=18,file='Resultados')

N=10000
imp=10
sum=0

doi=1,N
call random_number(r1)
r1=((100-20)*r1)+20
f=(1/(1+(r1*r1)))*(100-20)

sum=sum+f
prom=sum/i
res=mod(i,imp)
if(res.eq.0.0)then

write(18,*)i,prom,f
end if
end do
end3.











Se puede observar que para los tres diferentes tamaños del material se presenta la misma tendencia de distribución de energía, pues las curvas están prácticamentesuperpuestas, y también se observa que al aumentar la temperatura aumenta proporcionalmente la energía, y esto se debe a que los átomos se alteran adquiriendo mayor energía, debido a las vibraciones.4. De la gráfica del calor específico, estime la temperatura de transición.

Para los tres tamaños del material, el calor especifico aumenta en forma exponencial a una temperatura de 2.2, dondeadquiere el valor máximo, y luego disminuye de igual manera en forma exponencial, y esto se debe a que mientras aumenta un poco la temperatura los electrones saltan a la banda de valencia donde puedentransportar energía, pero al pasar de la temperatura de transición los electrones vibran de manera desordenada formando desordene impidiendo el transporte de energía.

5. Grafique el magnetismo y...