Ecuaciones Diferenciales De Primer Y Segundo Orden En Aplicaciones, Resueltas Con Matlab Usando Los Métodos De Runge Kutta, Euler Y Euler Modificado

3604 palabras 15 páginas
Universidad Nacional Experimental del Táchira.
Vicerrectorado Académico.
Decanato de Docencia.
Departamento de MATEMáTICA.
MATEMÁTICA ESPECIAL.
PROFESORA: Olga Roa.

Ecuaciones Diferenciales Ordinarias: Ejercicios

Realizado por:
Christian Santos C.I. 21221747
Ingeniería civil

San Cristóbal, 20 de abril de 2012.
Ejercicio 28.8: El compuesto A se difunde a través de un tubo de 4cm de largo y reacciona conforme se difunde. La ecuación que gobierna la difusión de la reacción es:
D∙d2Adx2-K∙A=0
En un extremo se encuentra una fuente grande de A con concentración de 0.1 M. En el otro extremo del tubo está un material que absorbe con rapidez cualquier A y hace la concentración sea 0 M. Si D=1.5x10-6 cm2/s y K=5x10-6 s-1
…ver más…
function EDO2O() ecu1=input('Ingrese la primera ecuación entre comillas en término de f(x,A,y): '); ecu2=input('Ingrese la segunda ecuación entre comillas en término de F(x,A,y): '); a=input('Ingrese el límite inferio del intervalo: '); b=input('Ingrese el límite superior del intervalo: '); in1=input('Ingrese la condición inicial de la primera ecuación diferencial A(0): '); in2=input('Ingrese la condición inicial de la segunda ecuación diferencial y(0): '); h=input('Ingrese el valor del tamaño de paso h: '); n=(b-a)/h; x=zeros(1,n+1); ye=zeros(1,n+1); Ae=zeros(1,n+1);
Aem=zeros(1,n+1);
yem=zeros(1,n+1);
Ark=zeros(1,n+1);
yrk=zeros(1,n+1); x(1)=a; Ae(1)=in1; ye(1)=in2; Aem(1)=in1; yem(1)=in2; Ark(1)=in1; yrk(1)=in2; for i=1:n w=x(i);A=Ark(i);y=yrk(i);a=eval(ecu1); ka=h*a; w=x(i);A=Ark(i);y=yrk(i);a=eval(ecu2); qa=h*a; w=x(i)+(h/2); A=Ark(i)+(ka/2); y=yrk(i)+(qa/2); b=eval(ecu1); kb=h*b; w=x(i)+(h/2); A=Ark(i)+(ka/2); y=yrk(i)+(qa/2); b=eval(ecu2); qb=h*b; w=x(i)+(h/2); A=Ark(i)+(kb/2); y=yrk(i)+(qb/2); c=eval(ecu1); kc=h*c; w=x(i)+(h/2); A=Ark(i)+(kb/2); y=yrk(i)+(qb/2); c=eval(ecu2); qc=h*c; w=x(i)+(h); A=Ark(i)+(kc); y=yrk(i)+(qc); d=eval(ecu1); kd=h*d; w=x(i)+(h); A=Ark(i)+(kc); y=yrk(i)+(qc); d=eval(ecu2); qd=h*d; x(i+1)=x(i)+h; %---------------------R.K-----------------------% l=(ka+(2*kb)+(2*kc)+kd)/6; Ark(i+1)=Ark(i)+l;

Documentos relacionados

  • Ecuaciones de primer y segundo grado
    1052 palabras | 5 páginas
  • Ecuaciones diferenciales ordinarias
    2977 palabras | 12 páginas
  • APLICACIONES DE LA ECUACIONES DIFERENCIALES DE 1ER ORDEN
    682 palabras | 3 páginas
  • Péndulo simple con ecuaciones diferenciales y aplicaciones
    1741 palabras | 7 páginas
  • Unidad 1: Ecuaciones Diferenciales De Primer Orden
    9524 palabras | 39 páginas
  • Metodo De Euler
    1403 palabras | 6 páginas
  • Rutina en c++
    672 palabras | 3 páginas
  • Diagramas de euler y venn
    1026 palabras | 5 páginas
  • Resonancia con ecuaciones diferenciales de segundo orden
    4592 palabras | 19 páginas
  • Ecuaciones diferenciales ejercicios de primer orden
    2551 palabras | 11 páginas