Help with Euler Method

I've been trying to solve this problem:

dydt =-5.33*10^-6*(-3.69*10^-6*y^4+2.33*10^-4*y^3 + 1.35*10^-3*y^2+5.42*10^-2*y+5.588)*(y+33)

h= 4200 second(stepsize)

y = temperature

boundaries for the temp should between 27 to -33 deg celsius

y(0) = 27.

help me solve this by using euler method using matlab. thanks