    {%- if default_stimulus_equations["equations"]|length > 0 %}
    boost::shared_ptr<RegularStimulus> {{class_name}}::UseCellMLDefaultStimulus()
    {
        // Use the default stimulus specified by CellML metadata
        {% for eq in default_stimulus_equations["equations"] %}{%- if not eq.lhs_modifiable %}const double {{ eq.lhs }} = {{ eq.rhs }}; // {{eq.units}}
        {% endif %}{% endfor %}boost::shared_ptr<RegularStimulus> p_cellml_stim(new RegularStimulus(
                -fabs({{default_stimulus_equations["membrane_stimulus_current_amplitude"]}}),
                {{default_stimulus_equations["membrane_stimulus_current_duration"]}},
                {{default_stimulus_equations["membrane_stimulus_current_period"]}},
                {% if default_stimulus_equations["membrane_stimulus_current_offset"] is defined %}{{default_stimulus_equations["membrane_stimulus_current_offset"]}}{%- else %}0.0{%- endif %}{% if default_stimulus_equations["membrane_stimulus_current_end"] is defined %}, {{default_stimulus_equations["membrane_stimulus_current_end"]}}{%- endif %}
                ));
        mpIntracellularStimulus = p_cellml_stim;
        return p_cellml_stim;
    }
    {%- endif %}