title: Implement Radau-IIA-method
body: <p>I´m looking for an implementation of the Radau-IIA-method with order 3 in MATLAB to solve a DAE system. This is my system (Rober problem):  </p>
  
  <pre><code>y1'=-0.04*y1+10^4*y2*y3
  y2'= 0.04*y1-10^4*y2*y3-3*10^7*y2^2
  0=y1+y2+y3-1
  </code></pre>
  
  <p>Initial values are </p>
  
  <pre><code>y1(0)=1 and y2(0)=y3(0)=0.
  </code></pre>
  
  <p>The Runge-Kutta tableau looks like this</p>
  
  <pre><code>1/3  |  5/12  -1/12
  1    |  3/4    1/4
  _____|_______________
       |  3/4    1/4
  </code></pre>
  
  <p>(Sorry, I don´t know if there is a possibility to write LaTeX or something in StackOverflow).</p>
  
  <p>I have no idea how to do that..</p>