Matlab - RESOLVER 2 ECUACIONES DIFERENCIALES ORDINARIAS SIMULTÁNEAMENTE POR EL MÉTODO DE RUNGE KUTTA 3 Y 4TO

function dydt=odet(t,y,param)

    %Descripción breve del contenido de la función...

    dydt=-0.002395882181697;

end

function ynew=runkuta2(t,y,h,param)

%% Algoritmo de Runge-Kutta de 2° Orden 

% Este programa es una implementación del método explícito de R-K 2er orden

% ...

%

% Llamada inicial de la EDO en los argumentos iniciales (t0,y0);

%

f0=feval('odet',t,y,param);

neq=length(y);  %Número de EDOs de 1er orden...

%

%% Se inicializan los parámetros del método:

%

A=[1;1];

B=zeros(3,2);

B(1,1)=1; B(1:2,2)=[1/2;1/2];

f=zeros(neq,3);

f(:,1)=f0;

hA=h*A;

hB=h*B;

f(:,2)=feval('odet',t+hA(1),y+f*hB(:,1),param);

%f(:,3)=...

tnew=t+hA(2);

ynew=y+f*hB(:,2);

f(:,3)=feval('odet',tnew,ynew,param);

end

% Simulación de un sistema de EDOs

%Inicializando el programa

%Inicializando:

clc,clear

%Datos y parámetros del sistema: 

t0=0;

yinic=15; %Hi, iniciales, en mts.

y=yinic;

%Parámetros de la simulación:

t=t0;

%tinic=t;

tmax=2520; %Tiempo total de simulación, en seg.

dt=360;  %Incremento de tiempo, en seg.

ndt=ceil((tmax-t0)/dt);

int=0;

tsim=zeros(1,ndt+1);

%Solución numérica de las EDOs: 

while int<ndt

    yy=runkuta2(t,y,dt,param); %Se aplica método de R-K 2do orden...

    tsim(int+2)=t+dt;

    ysim(:,int+2)=yy;

    t=tsim(int+2);

    y=yy;

    int=int+1;

end

tsim(1,1)=t0;

ysim(:,1)=yinic;

res=[tsim', ysim'];

%Resultados de la simulación:

format short g

disp(res)

%Salida de resultados a una tabla con descripción de contenido o enviar los

%resultados a una hoja de excel...

T = table(tsim',ysim', 'VariableNames', {'Tiempo','Alturadelnivel'});

T.Properties.Description = 'Variación del nivel del tanque cónico a través del tiempo';

T.Properties.VariableUnits =  {'segundos'  'metros'};

disp(T)

disp(T.Properties)

dH1/dt= - ((H1-H2)/423.624024)^0.5*(1/78.539846)

dH2/dt=((H1-H2)/423.624024)^0.5*(1/78.539846)