% difeq1.m  simple ODE
function dydt=difeq1(t,y)
dydt = t*t + 2.0*t + 3.0;
