IPN

IPN
IPN-ESIQIE

jueves, 11 de noviembre de 2010

program newton_r
implicit none
real::xo,tol,E,x,df,f
integer::max,I
f(xo)=xo**2+10*xo-1
df(xo)=2*xo+10
write(*,*)'dame el valor de xo,tol,max'
read(*,*)xo
read(*,*)tol
read(*,*)max
write(*,18)"Ite","xo","x","error"
do I=1,max,1
x=xo-(f(xo)/df(xo))
write(*,19)I,xo,x,abs(x-xo)
if(abs(x-xo)<=tol)then write(*,20)"x=",x stop else xo=x end if end do write(*,21)'No converge' 20 format(3/,tr10,a2,f10.3,3/) 21 format(a11) 19 format(I, 3e11.4) 18 format(4/,t11,a3,tr5,a2,tr9,a
,tr9,a5,2/)
end program

domingo, 10 de octubre de 2010

calcular presiones de saturacion de antoine de una sustancia con un itervalo dado

program antoine
!tarea 1
!=================
!carmona rodriguez andrea elizabeth 3Iv5
implicit none
character::sustancia
real::t,p,a,b,c,tmin,tmax
write(*,*)'dame el nombre de la sustancia'
read (*,*)sustancia
write(*,*)'dame el valor de tmin'
read (*,*)tmin
write(*,*)'dame el valor de tmax'
read (*,*)tmax
write(*,*)'dame el valor de a'
read (*,*)a
write(*,*)'dame el valor de b'
read (*,*)b
write(*,*)'dame el valor de c'
read (*,*)c
write(*,*)'dame el valor de t'
read (*,*)t
write(*,*)'A continuacion dse calculan las presion'
if(tmin<=t.and.t<=tmax)then
p=exp(a-(b/(t+c)))
write(*,*)'el valor de las presiones es',p
else
write(*,*)'la temperatura dada no pertenece al intervalo de temperaturas minimas y maximas'
end if
end program

calcular el valor de y por medio de programacion en fortran

program valor_de_y2
!=================
!carmona rodriguez andrea elizabeth 3Iv5
implicit none
real::x,y
write(*,*)'dame el valor de x'
read(*,*)x
if(x>=0)then
y=5*x**2+sqrt(x)
write(*,*)y,'es el valor solucion de y'
else if (-10<=x<0)then
y=24*x**3+exp(x)
write(*,*)y,'es el valor solucion de y'
else
y=exp(x)
write(*,*)y,'es el valor solucion de y'
end if
end program

por medio de un programa en fortran clacular presiones de saturacion de una sustancia x

program antoine
!tarea 1
!=================
!carmona rodriguez andrea elizabeth 3Iv5
implicit none
character::sustancia
real::t,p,a,b,c,tmin,tmax
write(*,*)'dame el nombre de la sustancia'
read (*,*)sustancia
write(*,*)'dame el valor de tmin'
read (*,*)tmin
write(*,*)'dame el valor de tmax'
read (*,*)tmax
write(*,*)'dame el valor de a'
read (*,*)a
write(*,*)'dame el valor de b'
read (*,*)b
write(*,*)'dame el valor de c'
read (*,*)c
write(*,*)'dame el valor de t'
read (*,*)t
write(*,*)'A continuacion dse calculan las presion'
if(tmin<=t.and.t<=tmax)then
p=exp(a-(b/(t+c)))
write(*,*)'el valor de las presiones es',p
else
write(*,*)'la temperatura dada no pertenece al intervalo de temperaturas minimas y maximas'
end if
end program

el de arriba es newton rhapson con do y este con do while

program newton_r
implicit none
real::xo,tol,E,x,df,f
integer::max,I
f(xo)=xo**2+10*xo-1
df(xo)=2*xo+10
write(*,*)'dame el valor de xo,tol,max'
read(*,*)xo
read(*,*)tol
read(*,*)max
I=1
E=1
do while(I<=max.and.E>tol)
x=xo-(f(xo)/df(xo))
E=abs(x-xo)
xo=x
I=I+1
end do
if(E<=tol)then
write(*,*)x
else
write(*,*)'no converge'
end if
end program

programa q calcula el valor de una ecuciaon de newton raphson

program n_rdo
implicit none
real::xo,tol,x,df,f
integer::max,I
f(xo)=xo**2+10*xo-1
df(xo)=2*xo+10
write(*,*)'dame el valor de xo,tol,max'
read(*,*)xo
read(*,*)tol
read(*,*)max
I=1
do I=1,max,1
x=xo-(f(xo)/df(xo))
xo=x
if(abs(x-xo)<=tol)then
write(*,*)x
stop
else
xo=x
end if
end do
end program

por medio del metodo del trapecio calcular el valor solucion a una integral de un intervalo numero dado

program metodo_trapecio
!carmona rodriguez andrea elizabeth
!3IV5
!6/10/10
IMPLICIT NONE
INTEGER::N,K
REAL::A,B,H
REAL::SUMA,INT,F,X0
REAL,ALLOCATABLE::X(:)
f(X0)=X0**2
WRITE(*,*)"DAME EL No Total de intervalos"
read(*,*)N
WRITE(*,*)"DAME LOS VALORES DE A Y B"
READ(*,*)A,B
H=(B-A)/N
ALLOCATE (X(0:N))
X(0)=A
X(N)=B
DO K=1,N-1,1
X(K)=X(K-1)+H
END DO
SUMA=0
DO K=1,N-1,1
SUMA=SUMA+F(X(K))
END DO
INT=H/2*(F(X(0))+2*SUMA+F(X(N)))
WRITE(*,*)"LA INTEGRAL ES=",INT
DEALLOCATE(X)
END PROGRAM