Jacek Kierzenka

There is a bug in example 4, which cause it to not follow the boundary values,
the correct res. function should be:
res = [ ya(1)
          yb(1)
          (ya(2)-yb(2)) ];
which is not dependent on T, instead of
res = [ya(1) - yb(1)
       ya(2) - yb(2)
       T*(-1/3)*(ya(1) - 0.7 + 0.8*ya(2)) - 1];

I can't understand how can I chose the right guess... and for examples the two guesses in example 1. Can somebody help me?

I was wondering to solve:

y1'=A1*y2'+B1*y3+C1*y4
y2'=A2*y1'+B2*y3+C2*y4
y'3=0
y'4=0

but I could not use the tutorial,