UNIVERSIDAD AUTÓNOMA DE NUEVO LEÓN
FACULTAD DE CIENCIAS FÍSICO - MATEMÁTICAS
FÍSICA COMPUTACIONAL
Tarea 5
DOCENTE: DR. ALFREDO TLAHUICE FLORES
INTEGRANTE GPO: 022
EDUARDO HIRAM PUEBLA ESCOBAR
1926361
02/03/2023
5. Calcule el volumen generado por función f(x) al girarlo con respecto al eje
x. Elija un intervalo en x de manera que se tenga un sólido hueco y un sólido
“macizo”. Haga el programa que calcule el volumen usando el método de
arandelas y el de cilindros. Use al menos 5 ejemplos donde compruebe que su
programa funciona para ambos tipos de sólidos de revolución. F(x) = 1/x
Objetivo: Hacer el programa usando la menor cantidad de funciones y subrutinas para a la
hora de querer usarlo no tener que hacer muchas modificaciones en el programa.
Primero se busca las fórmulas de los métodos que se piden en páginas web.
Formula del Método de Cilindro: FIG1
Formula del Método de Arandela: FIG2
Las FIG1 y FIG2 fueron sacadas de la pagina web: [Link]
%C3%A1lculo-2.0/section/6.5/primary/lesson/s%C3%B3lidos-de-revoluci%C3%B3n
%3A-vol%C3%BAmenes-mediante-capas-cil%C3%ADndricas-calc-spn/,
[Link]
%C3%B3lidos-de-revoluci%C3%B3n%3A-vol%C3%BAmenes-por-arandelas-calc-spn/.
PROGRAMA
En el programa se utilizó los tres métodos de integración para resolver las integrales:
Simpson, Doble Integral y se compararon con Romberg; Se declara la dimensión para el
programa de Romberg y una variable que se utilizara para el Método de Arandela (RECTA)
que se explicara más adelante.
También se intento declarar otra variable que se utiliza para el Meto de Arandela
(SIGNOS) pero al hacerlo las operaciones salían incorrectas por los decímales, nunca se
supo el porqué del error.
Después se piden los limites y N, aunque se vera que para la subrutina de Romberg se
reevalúa con N=20.
Ahora se empieza con el Método de Cilindro siguiendo las formulas antes mencionadas, se
mandan a llamar las dos subrutinas con las que se harán este método y así comparar los
resultados. (nótese que se llama “RombergIntegral1” porque no se logró hacer que se
utilizara la misma subrutina de Romberg para el método de cilindros y arandelas)
Ignórese la F2(x) que no se utilizara en el Método de Cilindros; En el código de Simpson se
menciona que las variables serán de primera precisión porque al poner REAL*8 los
números salían muy pequeños y no se llegaba al resultado, en el código no se modifica
nada mas excepto que la F(x) (nótese que la parte que se cambiara al querer usar otra
función será la parte “exiliada”) donde se utiliza la fórmula de Cilindros.
En el código de Romberg1 se reevalúa la N ya que al querer usar una N mayor los
resultados se iban alejando más y más, el código utiliza la misma función que usa el
método de Simpson y no se hizo más cambios al código, sin incluir los WRITE para hacer
el resultado más vistoso y fácil de entender.
Ahora toca el turno del método de Arandela donde antes de llamar a los CALL se decide
por un IF cual será el radio menor, es decir, la segunda curva a utilizar para la fórmula de
Arandela. (Aquí se utiliza la F2(x))
Después se evalúan las variables antes mencionadas (RECTA y SIGNOS) y se manda a
llamar una subrutina “Doble Integral” 1 o 2 dependiendo de cuál será la segunda curva para
utilizar, no se logró simplificar el código de Doble Integral para solo utilizar una sola
subrutina y que dependa de las variables del IF, como se hizo con el código de Romberg2.
En el código de Romberg2 se evalúa una tercera función F3(x) utilizando la formula del
método de Arandela y que depende de la función F2(x) (se logró que este código
dependiera de esta función al igual que lo hace la parte donde se evalúa cual será la
segunda curva en el IF), además aquí también se reevalúa la N.
Para el código de Doble Integral se intento hacer solo una subrutina y hacer un IF antes de
las funciones C(XZ), D(XZ) y dependiendo del resultado se utilizaría unas funciones para
el código u otras y así simplificar las líneas de código que se modificarían si se quisiera
utilizar otra función, pero el resultado seria el mismo asique ni se intento y solo se usaron
dos subrutinas donde la única diferencia es en cual seria la segunda curva y cual seria la
curva que estaría arriba de la otra.
En la ultima parte del codigo solo se multiplica por PI la integral y se cambia el formato de
los WRITE.
Para terminar se reutiliza unas lineas del codigo de Lagrange para repetir el proceso desde
el comienzo. (donde pide los limites)
EJEMPLO NO.1 = 1/X
Volumen con Hueco (Método de Arandela):
Ahora se grafica con Gnuplot la zona que girara en torno al eje x.
Ahora para poder visualizar mejor el volumen girado se utiliza la ayuda de un programa
externo en GeoGebra (Hecho por Laura del Rio):
Volumen sin Hueco (Método de Cilindro):
Ahora se grafica con Gnuplot la zona que girara en torno al eje x.
Ahora para poder visualizar mejor el volumen girado se utiliza la ayuda de un programa
externo en GeoGebra (Hecho por Laura del Rio):
EJEMPLO NO.2 =
Volumen con Hueco (Método de Arandela):
Ahora se grafica con Gnuplot la zona que girara en torno al eje x.
-
Ahora para poder visualizar mejor el volumen girado se utiliza la ayuda de un programa
externo en GeoGebra (Hecho por Laura del Rio):
-
Volumen sin Hueco (Método de Cilindro):
Ahora se grafica con Gnuplot la zona que girara en torno al eje x.
-
Ahora para poder visualizar mejor el volumen girado se utiliza la ayuda de un programa
externo en GeoGebra (Hecho por Laura del Rio):
-
EJEMPLO NO.3 =
Volumen con Hueco (Método de Arandela):
Ahora se grafica con Gnuplot la zona que girara en torno al eje x.
-
Ahora para poder visualizar mejor el volumen girado se utiliza la ayuda de un programa
externo en GeoGebra (Hecho por Laura del Rio):
-
Volumen sin Hueco (Método de Cilindro):
Ahora se grafica con Gnuplot la zona que girara en torno al eje x.
-
Ahora para poder visualizar mejor el volumen girado se utiliza la ayuda de un programa
externo en GeoGebra (Hecho por Laura del Rio):
-
EJEMPLO NO.4 =
Volumen con Hueco (Método de Arandela):
Ahora se grafica con Gnuplot la zona que girara en torno al eje x.
-
Ahora para poder visualizar mejor el volumen girado se utiliza la ayuda de un programa
externo en GeoGebra (Hecho por Laura del Rio):
-
Volumen sin Hueco (Método de Cilindro):
Ahora se grafica con Gnuplot la zona que girara en torno al eje x.
-
Ahora para poder visualizar mejor el volumen girado se utiliza la ayuda de un programa
externo en GeoGebra (Hecho por Laura del Rio):
-
EJEMPLO NO.5 =
Volumen con Hueco (Método de Arandela):
Ahora se grafica con Gnuplot la zona que girara en torno al eje x.
-
Ahora para poder visualizar mejor el volumen girado se utiliza la ayuda de un programa
externo en GeoGebra (Hecho por Laura del Rio):
-
Volumen sin Hueco (Método de Cilindro):
Ahora se grafica con Gnuplot la zona que girara en torno al eje x.
-
Ahora para poder visualizar mejor el volumen girado se utiliza la ayuda de un programa
externo en GeoGebra (Hecho por Laura del Rio):
-
LIBRETA
Se hizo la planificación del programa en una libreta para cumplir con el objetivo:
-
Se vera que al principio se intento usar una sola función y así solo modificar solo una línea
del código, pero no se logró, pero se pudo simplificar mucho las subrutinas y hacer que
dependieran de la misma función todas las posibles.
