http://www.chem.mtu.edu/~tbco/cm416/RKTutorial.html Study TI-89T CellSheet first To define plot characteristics, do this before cell program Home Diamond F1 F1:9:Format Enter Coordinates Rect Graph Order Seq Grid On Axes On Leading Cursor Off Labels On Mode Current Folder acel Exact/Approx Approximate Graph Function Split 1 App CellSheet Bring up CellSheet Name it x(arbitrary) and start the process below Define column size with and numbers size as desired with F3:Column format. I used Float 6 and cell width 6. is the url for this example. x | A | B | C | D | E | F | G | H | ---------------------------------------------------------------------------- 1 | tk | yk | dt | d1 | d2 | d3 | d4 | y(k+1) | ---------------------------------------------------------------------------- 2 | 0 | 1.0 | 0.1 | -0.1 | -0.0944 | -0.0955 | -0.0895 | 0.9051 | ---------------------------------------------------------------------------- 3 | 0.1 | 0.9051 | 0.1 | -0.0898 | -0.0842 | -0.0852 | -0.0796 | 0.82041 | ---------------------------------------------------------------------------- 4 | 0.2 | 0.82041 | 0.1 | -0.0797 | -0.0746 | -0.0754 | -0.0703 | 0.74539 | ---------------------------------------------------------------------------- 5 | 0.3 | 0.74539 | 0.1 | -0.0705 | -0.0659 | -0.0665 | -0.0621 | 0.67917 | ---------------------------------------------------------------------------- 6 | 0.4 | 0.67917 | 0.1 | -0.0621 | -0.0581 | -0.0586 | -0.0547 | 0.6208 | ---------------------------------------------------------------------------- 7 | 0.5 | 0.6208 | 0.1 | -0.0548 | -0.0512 | -0.0516 | -0.0483 | 0.56934 | ---------------------------------------------------------------------------- 8 | 0.6 | 0.56934 | 0.1 | -0.0483 | -0.0452 | -0.0456 | -0.0426 | 0.52393 | ---------------------------------------------------------------------------- ******************************************************************************** Note cell variables are not available in Home mode Put titles on row 1 a2 =0 b2=1 c2=0.1 d2 =c2*(-2*b2^2+e^(-3*a2)) e2 =c2*(-2*(b2+.5*d2)^2+e^(-3*(a2+.5*c2))) f2 =c2*(-2*(b2+.5*e2)^2+e^(-3*(a2+.5*c2))) g2 =c2*(-2*(b2+.5*e2)^2+e^(-3*(a2+.5*c2))) h2 =b2+1/6*d2+1/3*e2+1/3*f2+1/6*g2 A3 =a2+.1 B3 =h2 c3 =.1 d3 =c3*(-2*b3^2+e^(-3*a3)) e3 =c3*(-2*(b3+.5*d3)^2+e^(-3*(a3+.5*c3))) g3 =c3*(-2*(b3+.5*e3)^2+e^(-3*(a3+.5*c3))) h2 =b3+1/6*d3+1/3*e3+1/3*f3+1/6*g3 Use F3:3 Fill Range to extend column down a level one level at a time. Example Select B3 Enter Copy F3:3 Fill Range Formula =h2 Range b3:b4 Enter Enter b4 is now accurate. Repeat with rest of rows When all cells are accurate, plot the results. Plot F2:1:Plot Setup See plot menu F1 Define Plot 1 Plot Type Scatter Mark Box xRange a2:a8 yRange h2:h8 No Enter F5 ZoomData See graph See graph, scatter plot with squares, grid and labels x y Hit APPS Cell Sheet Current to return to cell layout