Eigenvalue problem associated with free surface waves

Hello,

I need to solve non homogeneous first order differential equation related to Rayleigh Waves.

How can I perform numerical integration in ode23 or ode45? I need recommendations. Thank you.

[Ritankar Sahu]

Hi,

I don't know if you are familier with the 'ode45' or not. So I am sharing a part of my one project I think you will get some idea from it. If afterwards you need some more, I can find for you.

Good Luck!!

% Initial Conditions
A=300;
B=1;
C=500;

[TOUT, YOUT] = ode45(@model, [0 1000], [A;B;C]); % that [0 1000] is the limit of simulation

figure("name", "A")
plot(TOUT, YOUT(:,3));
grid on;
title('Solution of State Space model with ODE45 for A');

function myFunc = model(TOUT, YOUT)
Fa=1.5;
Fb=0.5;
Fd=2;
Ta=300;
Tb=Ta;
H=10.^6;
c=5000;
h=500;

k=0.2*(1-exp(1-(X(3)/173)));

myFunc = [Fa+Fb-Fd; (1 - X(2)) / X(1) * Fa - X(2) / X(1) * Fb - X(2) * k; 
    (Ta - X(3)) / X(1) * Fa + (Tb - X(3)) / X(1) * Fb + H / (c * X(1)) - h / c * X(2) * k]; % Equations Matrix
end