Ver el código de «Deformaciones de un cuarto de anillo. Grupo 13-B»

{{ TrabajoED | Deformaciones de un cuarto de anillo. Grupo 13-B | [[:Categoría:Teoría de Campos|Teoría de Campos]]|[[:Categoría:TC22/23|2022-23]] | Felipe Delgado Rubio, Carlos Reillo González, Jorge Granadino Aranda, Diego de la Flor Bravo }}

== Enunciado ==

Consideramos la placa plana que ocupa un cuarto de anillo circular centrado en el origen y comprendido entre los radios 1 y 2, en el plano \(y≥|x|\). Dibujar con los ejes en el cuadrado \([-3,3]×[-1,3]\) tomando como función temperatura <math>T(x,y)=log((x-3)^2+2)</math>.
Tomaremos el siguiente campo: <math>\vec u(ρ,θ)=\frac{log(ρ)}{2}sin(2θ-π/2) ~</math>\(\vec e\)<sub>θ</sub>.

== Mallado ==

Dibujar un mallado que represente los puntos interiores del sólido. Tomar los ejes
(comando axis) en el rectángulo \((x, y) ∈ [-3; 3] × [−1; 3]\) y como paso de muestreo
h = 2/10 para las variables x e y.

Para ello, nos valemos del comando meshgrid, que nos define un mallado con los datos dados en el eje a representar (con el comando axis).

== Curvas de nivel de la temperatura==

Dibujar las curvas de nivel de la temperatura (comando contour) y decidir en qué punto la
temperatura es máxima a partir de la gráfica.