I would do something like this:
endtc = inputdlg('Maximum time for the simulation ');
endt = str2double(endtc{:});
t = linspace(0, endt, endt*10);
y = t.^2 .* exp(-t/2) .* sin(4*t); % Use Element-Wise Operations
Experiment to get the result you want.
