% mydsolve

% "old school" method:
ys = dsolve('D2y - 2*Dy -8*y = 0, y(0)=3, Dy(0)=6')

% "new school" alternative:
% syms y(t);
% Dy = diff(y, t);
% D2y = diff(Dy, t);
% ys = dsolve(D2y - 2*Dy -8*y == 0, y(0)==3, Dy(0)==6)

% ezplot (ys)

% Alternative plotting method (with more control)
x = 0:0.05:.5;
syms t;
y = subs(ys,t,x);
plot(x,y)
