%famous numbers

format long
   %display 15 digits
f='y';
   % because we are solving y'(x)=f(x,y) for f(x,y)=y,
   % and estimating y(1), since the solution to y'=y
   % with y(0)=1 is exp(x) and exp(1)=e.
x0 = 0; 
y0 = 1;
b=1;
n = 40;

 % Use Runge Kutta method to solve:
	[xR yR]=Runge_Kutta(x0,y0,b,n,f);
 % compare our estimate to actual value of e, as a column vector.
 ans1=[n;yR(n+1);exp(1)]
     % everything in Matlab is a matrix.  Here ans1 is a 3 by 1 matrix,
     % because I used semicolons between the quantities.  If I'd have used
     % commas Matlab would've made a 1 by 3 matrix.
 