MATLAB Ode45: Solving the Initial Value Problem dv/dt = -y

Using ode45 in MATLAB to Solve the Initial Value Problem dv/dt = -y

This code demonstrates how to use the ode45 function in MATLAB to solve the initial value problem dv/dt = -y with the initial condition y(0) = 1 over the time interval 0 to 5.

Code:

function dydt = myode(t,y)
  dydt = -y;
end

[t,y] = ode45(@myode,[0 5],1);
plot(t,y)
xlabel('t')
ylabel('y')
title('Solution of dv/dt = -y, y(0) = 1')

Explanation:

1. Function Definition: The first part defines a function myode that represents the differential equation. It takes the time t and the current state y as input and outputs the derivative of the state dydt:

   function dydt = myode(t,y)
     dydt = -y;
   end
   

2. Solving the ODE: The ode45 function is used to solve the differential equation numerically. It takes the following arguments:

   • @myode: The function handle of the differential equation (myode in this case).
   • [0 5]: The time interval over which the solution is calculated (from 0 to 5).
   • 1: The initial value of the state variable y at time 0.

   The function returns two arrays: t containing the time points and y containing the corresponding state values.

   [t,y] = ode45(@myode,[0 5],1);
   

3. Plotting the Solution: Finally, the code plots the solution using the plot function. The xlabel, ylabel, and title functions are used to label the axes and provide a descriptive title for the plot.

   plot(t,y)
   

xlabel('t') ylabel('y') title('Solution of dv/dt = -y, y(0) = 1')

This code demonstrates how to use MATLAB's `ode45` function to solve simple initial value problems. You can modify the `myode` function to represent different differential equations and adjust the time interval and initial condition to solve a variety of problems.