MATLAB: Matlab coding for the below

Question

r(t) =1 − e−γ2t ,t ≤ Ta
      4e−γ1(t−Tb)(sin(2πωat) + 1), otherwise
t=500
Ta=250,Tb=500/2.2,γ1=0.07,γ2=0.004,ωa = 0.06,

Answer 1

Ta=250 ;
Tb=500/2.2 ;
t=1:500 ;
g1=0.07 ;
g2=0.004 ;
Wa = 0.06 ;
r = zeros(length(t),1) ;
r(t<=Ta) = 1-exp(1)-g2*t(t<=Ta) ;
r(t>Ta) = 4*exp(1)-g1*(t(t>Ta)-Tb).*(sin(2*pi*Wa*t(t>Ta))+1) ;
plot(t,r)