Write a matlab script along with required functions, that solve the following system of differential equations on time horizon [0,1]

y'''(t)+ay(t)=-by''(t)+u(t)

Where u(t) is any function of time

Please explain your solution to be a valid answer