Starting diary file "diary-2016-02-22.txt" on Mon Feb 22 2016 12:26 PM.

edit weightf

tspan = [0 10*pi]
tspan =
         0   31.4159

ode23tx (@weightf, tspan, [1/2;0]);
axify

ode23tx (@weightf, tspan, [1 + 1/2;0]);
axify

ode23tx (@weightf, tspan, [1 ;0]); axify;
ode23tx (@weightf, tspan, [1 ; -1/4]); axify;
ode23tx (@weightf, tspan, [1 ; 10]); axify;
ode23tx (@weightf, tspan, [.5 ; 0]); axify;

[T,Y] = ode23tx (@weightf, tspan, [.5 ; 0]);

whos
  Name         Size            Bytes  Class     Attributes

  T          125x1              1000  double              
  Y          125x2              2000  double              
  tspan        1x2                16  double              

T
T =
         0
    0.0002
    0.0010
    0.0050
   ...
   31.0668
   31.3120
   31.4159
Y
Y =
    0.5000         0
    0.5000    0.0001
    0.5000    0.0005
    0.5000    0.0025
    ...
    0.5455   -0.1602
    0.5203   -0.0451
    0.5182    0.0049

plot(T,Y)

h = diff(T)
h =
    0.0002
    0.0008
    0.0040
    0.0200
    ...
    0.2881
    0.2451
    0.1039
whos
  Name         Size            Bytes  Class     Attributes

  T          125x1              1000  double              
  Y          125x2              2000  double              
  h          124x1               992  double              
  tspan        1x2                16  double              

plot(h),shg

clc

edit ode1

ode1 (@weightf, tspan, [.5 ; 0], .1);
Took 315 steps

ode1 (@weightf, tspan, [.5 ; 0], .01); 
Took 3142 steps
axify;


ode1 (@weightf, tspan, [.5 ; 0], .001); axify;
Took 31416 steps

diary off