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;