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