function Explicit(fname,t0,y0,h,n)

tvals=[];yvals=[];

tc = t0;
yc = y0;
tvals = tc;
yvals = yc;
fc = fname(tc,yc);
for j = 1:n
  yc = yc + h*fc;
  tc= tc+h;
  fc = fname(tc, yc);
  yvals = [yvals,yc];
  tvals = [tvals,tc];
end
plot2d(tvals, yvals)


function Implicit(fname,t0,y0,h,n)

tvals=[];yvals=[];

tc = t0;
yc = y0;
tvals = tc;
yvals = yc;
fc = fname(tc,yc);
for j = 1:n
  deff('[y]=fsol1(x)', 'y=-x+yc+h*fname(tc,x)');
  yc = fsolve(yc+h*fc,fsol1);
  tc= tc+h;
  fc = fname(tc, yc);
  yvals = [yvals,yc];
  tvals = [tvals,tc];
end
plot2d(tvals, yvals)


// function [yp]=f2(t,y)
// yp=[];
// yp = -50*(y-cos(t));