def simulateAndPlot(r, **kwargs):
    """ Simulate with r.simulate with given arguments and plot with tellurium.

    :returns: simulation result
    :rtype: NamedArray
    """
    s = r.simulate(**kwargs)
    r.plot(s)
    return s