Matlab:Evaluating ODE at Points

  • MATLAB
  • Thread starter qspeechc
  • Start date
  • Tags
    Ode Points
  • #1
qspeechc
843
15
Hello everyone.

I have a second order ODE, which I reduced to two first order ODEs, and solved with ode45. Now I want to evaluate the solution and the derivative at certain points (specifically, I want to plot the Poincare section).

So, what I did is this:

Code:
>> sol=ode45(@...);
>> xint=a:b;    
>> [Sxint, Spxint]=deval(sol, xint);
>> A=[Sxint, Spxint];
>> plot(Sxint, Spxint,'o')

Where:
xint is the points at which I want to evaluate the solution.
Sxint is the value of the solution at the points xint.
Spxint is the value of the derivative at those points.

But this is clearly wrong. I know the solution should be periodic, but I get something random. I think it's the plotting I'm doing incorrectly. Any help appreciated.
 
Physics news on Phys.org
  • #2
Ok, instead, I created a separate M-file:

Code:
hold on,
for i=1:length(xint)
    plot(Sxint(:,i), Spxint(:,1),'.')
end

Is this correct?
 
  • #3
Nevermind, I've figured it out, you can close the thread.
 

Similar threads

Replies
5
Views
2K
Replies
3
Views
2K
Replies
0
Views
2K
Replies
7
Views
3K
Replies
8
Views
2K
Replies
0
Views
1K
Replies
7
Views
3K
Back
Top