system of ode plotting solution

1 view (last 30 days)
Shubham
Shubham on 19 Feb 2021
Commented: Shubham on 20 Feb 2021
function qdot = vdPav(t,q)
Delta= -2;
F= 1;
epsilon = 0.1;
qdot = [epsilon((1/2)*q(1)-(1/8)*q(1)^3+(1/2)*Delta*q(2)-(1/8)*q(1)*q(2)^2-(1/2)*F);epsilon((1/2)*q(2)-(1/8)*q(2)^3+(1/2)*Delta*q(1)-(1/8)*q(2)*q(1)^2)];
op = odeset('relTol',1e-7,'abstol',1e-7);
[t,q] = ode15s('vdPav',[0 50],[0.1 0.1],op);
plot(q(:,1),q(:,2));
I am trying to plot q(1) v/s q(2) but I am getting error
  5 Comments
Shubham
Shubham on 20 Feb 2021
Yes I rectified my mistake thanks a lot

Sign in to comment.

Answers (0)

Tags

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!