% Parámetros del circuito RLC
R = 100; % Resistencia en ohmios
L = 0.1; % Inductancia en henrios
C = 0.001; % Capacitancia en faradios
V0 = 10; % Voltaje inicial en voltios
% Definir la ecuación diferencial del circuito RLC
equation = @(t, y) [(V0 - y(1)/C - R*y(2))/L; y(1)];
% Condiciones iniciales
initialConditions = [0; 0];
% Tiempo de simulación
tspan = [0 0.5];
% Resolver la ecuación diferencial
options = odeset('Events', @(t, y) y(2), 'RelTol', 1e-6, 'AbsTol', 1e-9);
[t, y, te, ye, ie] = ode45(equation, tspan, initialConditions, options);
% Graficar la carga del condensador
figure;
plot(t, y(:, 1));
hold on;
% Marcar los valores máximos
[maxValues, maxIndices] = findpeaks(y(:, 1));
plot(t(maxIndices), maxValues, 'ro', 'MarkerFaceColor', 'r');
% Añadir etiquetas y título
title('Carga del Condensador en un Circuito RLC');
xlabel('Tiempo (s)');
ylabel('Voltaje en el Condensador (V)');
legend('Carga del Condensador', 'Valores Máximos');
% Mostrar los valores máximos
disp('Valores Máximos:');
disp([t(maxIndices), maxValues]);