function [y, x] = myeuler(xstart, ystart, ydotstart, stepsize, xstop)
steps = (xstop-xstart)/stepsize;
cury = ystart;
curx = xstart;
ct = 0;
ydot = ydotstart;
for s = 1:steps
	ct = ct+1;
	y(:,ct) = cury;
	x(:,ct) = curx;
	curx = curx+stepsize;
	newydot = ydot - stepsize.*cury;
	cury = cury + (ydot.*stepsize);
	ydot = newydot;
end
end
