% Calculate the unit impulse response of the gain. 
simType('symbolic');                        % Symbolic simulation
dataType('impulse');                        % Calculate the unit impulse response of the gain (gain type has already been set)
gainImpulse = execute();                    % Assign the execution result to the variable 'gainImpulse'               
gainImpulseExpr = gainImpulse.results(1);   % Unpack the execution result and assign the expression to 'gainImpulseExpr'
htmlPage('Unit impulse response');          % Start a new html page
head2html('Symbolic expression');           % Put a heading on it
syms('h(t)');                               % Define the left hand side of the equation
eqn2html(h(t), gainImpulseExpr);            % Write the equation to the html page
% generate the plot of h(t)
simType('numeric');                         % This requires numeric simulation                  
figImpulse = plotTime('Impulse response', execute(), 0, 1e-3, 100); % Make the time-domain plot
head2html('Plot');                          % Put a heading on the page
fig2html(figImpulse, 'figImpulse.svg', 800);% Put the plot on the page