Archivo discusión:Ejerciciopolo.jpeg

% Paso de muestreo h = 0.1;

% Valor de u (radios) u = 1:h:2;

% Calculamos cuántos puntos caben aproximadamente puntos = round(pi/h) + 1;

% Vector v (ángulos) de 0 a pi v = linspace(0, pi, puntos);

% Matrices de coordenadas polares [rho, theta] = meshgrid(u, v);

% Conversión a Cartesianas xx = rho .* cos(theta); yy = rho .* sin(theta);

% Crear figura figure;

% Dibujar líneas radiales y circunferencias plot(xx, yy, 'r'); hold on plot(xx', yy', 'r');

% Ajustes visuales axis equal axis([-3, 3, -0.5, 2.5]); title('Sección de arco con campo vectorial radial'); xlabel('X'); ylabel('Y'); grid on

% --- Campo vectorial radial --- r1 = 1; r2 = 2; theta_vec = linspace(0, pi, 20); % menos puntos para flechas rho_vec = linspace(r1, r2, 20); [R, THETA] = meshgrid(rho_vec, theta_vec); X = R .* cos(THETA); Y = R .* sin(THETA);

% Componente radial U_rho = (1/5) * (R - 1) .* R; U_theta = zeros(size(U_rho));

% Convertir a cartesianas U = U_rho .* cos(THETA) - U_theta .* sin(THETA); V = U_rho .* sin(THETA) + U_theta .* cos(THETA);

% Graficar flechas quiver(X, Y, U, V, 0.5, 'b', 'LineWidth', 1.2); % azul, grosor 1.2

hold off