This question is closed. Reopen it to edit or answer.
can anyone help me with "ode45" ?
1 view (last 30 days)
i do have a function like;
and a "ode45" like
I would like to take "NOD" out given in my function;
How can I get it in a ".mat" ? is that possible? Or I need to use the Runge Kutta method manually?
Please a bit urgent, thanks in advance. If you want a clear understanding, I can upload a draft code.
Star Strider on 25 Nov 2018
Probably the easiest way would be to integrate your function first, the using the solved values for ‘t’ and ‘x’, run it again, recover the ‘NOD’ output, and save it.
xdot(1,:) = t-x(1); % Create Equation
xdot(2,:) = t+x(2); % Create Equation
NOD = rand; % Create ‘NOD’
tspan = [0 1];
x0 = [1;1];
[T,X] = ode45(@random_func, tspan, x0);
NODv = zeros(numel(T),1); % Preallocate
for k1 = 1:numel(T)
[~,NOD(k1,:)] = random_func(T(k1),X(k1,:)); % Calculate ‘NOD’, Save As Vector
fprintf(1, '%.4f,%.4f\n', [T NOD]') % Save Data
I simply printed out the ‘NOD’ vector here (alternatively saving it as a ‘.csv’ file with the same statement although with an appropriate file identification number created by fopen), since I did not want to create the file and then have to delete it. Use the appropriate save syntax to create your ‘.mat’ file.
This illustrates the idea. You will likely need to experiment with it and tweak it to make it work with your functions.