















Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity
Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium
Prepara tus exámenes
Prepara tus exámenes y mejora tus resultados gracias a la gran cantidad de recursos disponibles en Docsity
Prepara tus exámenes con los documentos que comparten otros estudiantes como tú en Docsity
Los mejores documentos en venta realizados por estudiantes que han terminado sus estudios
Estudia con lecciones y exámenes resueltos basados en los programas académicos de las mejores universidades
Responde a preguntas de exámenes reales y pon a prueba tu preparación
Consigue puntos base para descargar
Gana puntos ayudando a otros estudiantes o consíguelos activando un Plan Premium
Comunidad
Pide ayuda a la comunidad y resuelve tus dudas de estudio
Descubre las mejores universidades de tu país según los usuarios de Docsity
Ebooks gratuitos
Descarga nuestras guías gratuitas sobre técnicas de estudio, métodos para controlar la ansiedad y consejos para la tesis preparadas por los tutores de Docsity
ejercicios sobre la materia de metodos numéricos
Tipo: Ejercicios
1 / 23
Esta página no es visible en la vista previa
¡No te pierdas las partes importantes!
En oferta
Instituto Tecnológico de Matehuala MÉTODOS NUMÉRICOS
Problemas únicos
ENUNCIADO DEL PROBLEMA
Solución de código en Matlab
Código de Matlab
c=0; g=0; %Este for obtiene y guarda todos los valores de t %También se utiliza para evaluar la ecuación diferencial for p=a:h:b d=1+d; t(d)=p; v(d)=subs(m,p); end %Este for se usa para contabilizar las iteraciones for s=c:1:(d-1) g=1+g; k(g)=(g-1); end %Este for obtiene los valores aproximados de solución for j=a:h:(b-h) i=1+i; w(i+1)=w(i)+((h/4)(subs(f,{x,y},{t(i),w(i)})))+(((3/4)h)(subs(f,{x,y},{(t(i) +((2/3)h)),(w(i)+(((2/3)h)(subs(f,{x,y},{t(i),w(i)}))))}))); end %Presentación de los datos en forma de tabla utilizando matrices [k' t' w' v']
funcion de trabajo F(Ti,Wi)=Dy/Dx formula del metodo: Wi+1=Wi+h/4F(Ti,Wi)+3h/4F(Ti+2h/3, Wi+2h/3*F(Ti,Wi))
fprintf(' ----------------------\n') fprintf(' MÉTODO DEL PUNTO MEDIO\n') fprintf(' ----------------------\n') fprintf('\n'); syms x y disp('Wi+1=Wi+hF(Ti+h/2,Wi+h/2F(Ti,Wi)) la funcion de trabajo Dy/Dx=F(Ti,Wi)') fprintf('\n'); d=input(' - Introduzca la ecuación diferencial encerradas en tildes y que Dy no este x y : '); n=input(' - Introduzca la condición y(a)=b encerradas en tildes : '); l=input(' - Introduzca la variable encerradas en tildes : '); f1=input(' - Introduzca la función de trabajo : '); ya=input(' - Introduzca la condición inicial : '); a=input(' - Introduzca el valor de a : '); b=input(' - Introduzca el valor de b : '); h=input(' - Introduzca el tamaño de paso h : '); fprintf('\n\n'); fprintf(' - La solución de la ecuación diferencial es : \n\n\n'); m = dsolve(d,n,l); pretty(m); fprintf('\n\n\n\n'); fprintf(' i ti wi y(t)'); %Condiciones para el funcionamiento de los lazos FOR f=f1; w(1)=ya; i=0; t(1)=a; v(1)=a; d=0; c=0; g=0; %Este for obtiene y guarda todos los valores de t
%También se utiliza para evaluar la ecuación diferencial for p=a:h:b d=1+d; t(d)=p; v(d)=subs(m,p); end %Este for se usa para contabilizar las iteraciones for s=c:1:(d-1) g=1+g; k(g)=(g-1); end %Este for obtiene los valores aproximados de solución for j=a:h:(b-h) i=1+i; w(i+1)=w(i)+(h(subs(f,{x,y},{(t(i)+h/2),w(i)+((h/2)(((subs(f,{x,y}, {t(i),w(i)})))))}))); end %Presentación de los datos en forma de tabla utilizando matrices [k' t' w' v']