Hi i've looked for help with this problem all over but found little to no help.
I'm simply trying to model a forced Van der Pol oscillator (ie instead of making diff equation = 0 now make it =cos(wt)). I need to update my function "vdp" 's 't' term as i run the code. I don't see how i can do this (ie i want the t term in the cos function to increment for each execution). Please any help would be greatly appreciated.
function dydt = vdp(t,y,Mu)
% End solvevdp.m
t = [1:100:1]
title('van der Pol Equation, mu = 1000')
axis([0 3000 -3 3])