metodos numericos ejercicios
Publicado por josefina (6 intervenciones) el 08/01/2021 04:16:03
necesito ayuda con este ejercicio, necesito encontrar las ecuaciones diferenciales para luego llevarlas a un sistema de ecuaciones diferenciales de primer grado y aplicar un metodo numerico
uso el siguiente codigo, donde antes debo definir los siguientes parametros, mi funcion f, un intervalo, la condicion inicial y mi L.
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
function [t,y]=rk4(f, inter, y0, L)
%function [t,y]=rk4(f, [T0 Tf], y0, L)% Metodo de Runge-Kutta de orden 4% para resolver% y' = f(t,y) en [t0,TF]
% y(t0) = y0% Usando L pasos% y0 puede ser vectorial o escalart = linspace(inter(1),inter(2),L+1)';h = (inter(2)-inter(1))/L;
% reservamos lugar en memoria para yy = zeros( length(y0), L+1 );
y(:,1) = y0;
for n=1:L
k1 = h * f(t(n) , y(:,n));
k2 = h * f(t(n)+h/2, y(:,n)+k1/2);
k3 = h * f(t(n)+h/2, y(:,n)+k2/2);
k4 = h * f(t(n+1), y(:,n)+k3);
y(:,n+1) = y(:,n) + (k1+2*k2+2*k3+k4)/6;
end
y=y';
Valora esta pregunta
0