AQUI SE LOS DEJO ... ESPERO QUE USTEDES TAMBIEN AYUDEN A LA COMUNIDAD function f f=input('ingrese la ecuacion diferencial:','s'); x0=input('ingrese el valor inicial de x:'); x1=input('ingrese el valor final de x:'); y0=input('ingrese el valor de y:'); n=input('ingrese el numero de iteraciones:'); h=(x1-x0)/n; xs=x0:h:x1; fprintf(' ''´Pto x0 y(x1)'); for i=1:n Pto=i; x0=xs(i+1); x=x0; y=y0; k1=h*eval(f); x=xs(i); y=y0+k1; k2=h*eval(f); y0=y0+(k1+k2)/2; fprintf(' %2.0f%10.6%.10.6f ',Pto,x0,y0); end fprintf(' El punto aproximado Y(x1) es= %8.6f ',y0);
Una pregunta, que método o cómo puedo encontrar la solución analítica de esa ecuación diferencial? Cuando "y" esta dentro de la raiz? Eso hace que sea una ecuación diferencial ordinaria de primer orden no lineal. En Matlab me tira que no se puede obtener la solución explícita, es decir en forma de ecuación :/
que tal, para hacer la grafica usando plot, cuales deberian ser los valores que deverian ir dentro de la funcion plot? tanto para que de la grafica de las iteraciones junto con la grafica original de la funcion que se ingreso
muchas gracias, amigo! Comprendí mejor el método gracias a tu video. Una consulta. Realizo los pasos en el orden tal cual, y bajo el mismo código, pero sólo me arroja un resultado; es como si el - if - no avanzara. Gracias!
Ricardo Fuentes gracias por comentar, el código debes copiarlo en el orden que está, y la verdad no utilice en este script la condicional “if”.. podrías revisarlo nuevamente y observar cuál es el detalle que quieres aclarar?
Gabriel Zurita gracias por tu comentario, te sugiero que veas el vídeo detenidamente y compares tu script para que esta igual en las funciones, en los “;” y sobre todo en los colores que es una buena guía para que el script no tenga errores. Además puedes revisar los comentarios para que verifiques que si funciona el script
anderson mora termine haciéndolo a mano jajaja estaba estudiando para mi parcial de métodos numéricos pero lo bueno que gane esa clase sin tanto problema
AQUI SE LOS DEJO ... ESPERO QUE USTEDES TAMBIEN AYUDEN A LA COMUNIDAD
function f
f=input('ingrese la ecuacion diferencial:','s');
x0=input('ingrese el valor inicial de x:');
x1=input('ingrese el valor final de x:');
y0=input('ingrese el valor de y:');
n=input('ingrese el numero de iteraciones:');
h=(x1-x0)/n;
xs=x0:h:x1;
fprintf('
''´Pto x0 y(x1)');
for i=1:n
Pto=i;
x0=xs(i+1);
x=x0;
y=y0;
k1=h*eval(f);
x=xs(i);
y=y0+k1;
k2=h*eval(f);
y0=y0+(k1+k2)/2;
fprintf('
%2.0f%10.6%.10.6f
',Pto,x0,y0);
end
fprintf('
El punto aproximado Y(x1) es= %8.6f
',y0);
tienes un error en la línea antes del "end" la línea correcta es: fprintf('
%2.0f%10.6f%10.6f
',Pto,x0,y0);
Muchas gracias!! Te mereces mas seguidores por todo tu trabajo ... 100000 likes
por que solo llegas hasta k2 si es de 4 etapas
Muy buen video amigo
Disculpa y para sacar la gráfica como debería hacerlo?
porqué solo usas k1 y k2?
se supone que es de 4to orden?
veo que ahi tienes el euler mejorado, lo compartes ?
grandeeee
Una pregunta, que método o cómo puedo encontrar la solución analítica de esa ecuación diferencial? Cuando "y" esta dentro de la raiz? Eso hace que sea una ecuación diferencial ordinaria de primer orden no lineal. En Matlab me tira que no se puede obtener la solución explícita, es decir en forma de ecuación :/
a mi me da un error en la linea 15 y nose porque??? si lo escribi tal cual
genial¡¡¡¡
que tal, para hacer la grafica usando plot, cuales deberian ser los valores que deverian ir dentro de la funcion plot?
tanto para que de la grafica de las iteraciones junto con la grafica original de la funcion que se ingreso
Ayuda, cómo puedo graficar los valores con la funcion?
bro estas haciendo el codigo de segundo orden
Sería bueno que hagas el código de RK de orden 2, sería mas sencillo (?)
ALDAIR WILLIAN BIGMAN OLAYA CHIRA gracias por comentar, el script de 4to orden es más complicado que el de 2do orden, por eso lo hice
tienes el archivo en math lab?
muchas gracias, amigo! Comprendí mejor el método gracias a tu video. Una consulta. Realizo los pasos en el orden tal cual, y bajo el mismo código, pero sólo me arroja un resultado; es como si el - if - no avanzara.
Gracias!
Ricardo Fuentes gracias por comentar, el código debes copiarlo en el orden que está, y la verdad no utilice en este script la condicional “if”.. podrías revisarlo nuevamente y observar cuál es el detalle que quieres aclarar?
Me pasa lo mismo , es como si nunca entrara el for
sale error en la linea 16 .... es la parte donde se k1=h*eval(f);
tengo el mismo error, lo pudiste solucionar?
@@eduardohuchimmena1012 pon f entre comas eval('f')
Hi,I need rungekutta to solve delay differential equation
Please helpe me I need code
me sale error en la line 16 .....alguien pudo hacerlo correctamente?
Gabriel Zurita gracias por tu comentario, te sugiero que veas el vídeo detenidamente y compares tu script para que esta igual en las funciones, en los “;” y sobre todo en los colores que es una buena guía para que el script no tenga errores. Además puedes revisar los comentarios para que verifiques que si funciona el script
Puedes subir el programa de Runge-Kutta a mega o a algún servidor para solamente descargarlo y no tener que copiar todo eso?
ExtremE GameR has tu tarea, no esperen que te lo hagan, para que quieres que lo suban, para copiar y pegar. Jajaja
entra a mi pagina lo tengo en .exe ya listo para usar
y la grafica papu? y la comparación con el método analítico papu?
bajale el volumen a tu intro so gei
Como le podria hacer para graficar el resultado de la ecuacion diferencial?
Con ese mismo codigo
me da error en la funcion, ni para copiar sirvo jaja...
ayuda.
fredy peralta te sugiero que pauses el video y copies el script con tranquilidad :D
Tutoriales Básicos Ec lo pause, time screenshoot y lo copié pero no me corrió jajajajaja ayuda!!!!
fredy peralta fijate las filas 23 y 21, yo cometí error a poner comas por puntos. Cuando decides las cantidades de cifras décimales
anderson mora termine haciéndolo a mano jajaja estaba estudiando para mi parcial de métodos numéricos pero lo bueno que gane esa clase sin tanto problema