Accumulation, Exposure, and Flow-to-Stock Reasoning

Last Updated June 15, 2026

Accumulation, exposure, and flow-to-stock reasoning explain how local flows become cumulative system consequences. A flow is a quantity moving, changing, emitted, absorbed, consumed, stored, loaded, or transferred per unit time. A stock is the accumulated result of past flows. Exposure is the cumulative contact, burden, intensity, or dose experienced by a person, population, infrastructure system, ecosystem, or institution over an interval.

In systems modeling, this relationship is foundational. Emissions accumulate into atmospheric burden. Inflows and outflows determine resource stocks. Traffic volumes accumulate into infrastructure load. Concentrations and time combine into exposure. Financial contributions accumulate into balances. Infection rates accumulate into case counts. Wear, stress, and usage accumulate into degradation.

This article develops accumulation, exposure, and flow-to-stock reasoning as a modeling principle. It examines net flows, stocks, cumulative exposure, dose, residence time, inflow-outflow balance, units, lags, measurement windows, numerical approximation, audit trails, and responsible interpretation of cumulative claims.

Accumulation is one of the central reasons calculus matters for systems modeling. A system’s present state often reflects not only what is happening now, but what has been flowing into, out of, or through the system over time. Exposure extends the same logic to contact and burden: what matters is not just intensity at one moment, but intensity accumulated across duration, population, location, vulnerability, or pathway.

Why Flow-to-Stock Reasoning Matters

Flow-to-stock reasoning matters because many system outcomes are cumulative. A stock is not explained by a single moment of inflow or outflow. It is explained by the history of net movement into and out of the stock. A cumulative exposure is not explained by a single peak intensity. It is explained by intensity over time, sometimes weighted by vulnerability, pathway, population, location, or duration.

The core relationship is:

If the system has separate inflows and outflows, the net-flow form is:

This relationship appears across systems modeling. Atmospheric concentration reflects emissions, removals, and residence time. Reservoir storage reflects inflows, releases, evaporation, and withdrawals. Inventory reflects production, shipments, losses, and returns. Public health burden reflects exposure intensity, duration, population vulnerability, and pathway. Infrastructure load reflects traffic, weight, frequency, maintenance, and degradation.

Flow-to-stock reasoning prevents a common error: focusing on rates while ignoring accumulation. A low annual rate can create a large stock over a long interval. A high peak exposure can matter less than moderate exposure sustained for a long duration. A system can look stable for a time while accumulated stress moves it toward a threshold.

The key question is not simply “What is the current rate?” It is “What has accumulated, what is still accumulating, what is leaving, and what stock or exposure does that history produce?”

Back to top ↑

Flows, Stocks, and Accumulation

A flow measures change per unit time. A stock measures an amount at a point in time. The relationship between them is one of the most important bridges between calculus and systems thinking.

Integrating the rate recovers the stock trajectory:

In systems modeling, this relationship is often more important than the formula alone. It says that stocks contain memory. The current amount of a stock is shaped by previous flows. This is why stock variables can create delays, persistence, inertia, and path dependence. Even if a flow changes immediately, the stock may respond gradually because it reflects accumulated history.

Flows and stocks also differ in interpretation:

A model that reports flows without stocks may miss cumulative consequence. A model that reports stocks without flows may obscure the processes driving change. Responsible flow-to-stock reasoning connects both.

Back to top ↑

Net Flow and Stock Change

Many stocks are governed by both inflows and outflows. The stock changes according to the net flow:

The stock equation becomes:

Integrating over an interval gives:

This accounting identity is simple but powerful. It reveals why reducing inflow may not immediately reduce the stock if the inherited stock is large. It explains why outflow capacity matters as much as inflow control. It also clarifies why a temporary imbalance can produce long-lasting consequences when the stock persists.

For example, emissions reductions may slow the growth of atmospheric burden without immediately lowering accumulated concentration. A resource stock may continue declining if extraction remains above regeneration. A backlog may keep growing even if processing rates improve but remain below arrival rates. A debt stock may increase if interest and new borrowing exceed repayment.

Net-flow reasoning should also include sign conventions. A positive flow may mean inflow, outflow, emissions, extraction, or depletion depending on the model. Ambiguous signs are a major source of modeling error. Every stock-flow model should specify whether positive values increase or decrease the stock.

Back to top ↑

Exposure as Accumulated Contact

Exposure is a form of accumulation. It combines intensity with duration and often with pathway, location, population, or vulnerability. A simple exposure measure is:

If vulnerability, population, or pathway weighting matters, a weighted exposure measure may be:

Exposure reasoning is central in public health, environmental justice, occupational safety, infrastructure stress, climate risk, and social burden analysis. It shifts attention from momentary intensity to cumulative contact. A short high-intensity event and a long moderate-intensity exposure can produce different burdens even if their peak values look very different.

Exposure models require careful interpretation. What is being accumulated? Concentration? Dose? Time above threshold? Population-weighted burden? Vulnerability-weighted burden? Area under a curve? A simple integral may not capture nonlinear dose-response, recovery, adaptation, saturation, or threshold effects. Still, accumulation provides the first disciplined structure for making exposure auditable.

Exposure is not just intensity. It is accumulated contact under a stated interpretation.

Back to top ↑

Dose, Intensity, and Duration

Dose-like quantities often arise when exposure intensity accumulates over time. If \(r(t)\) is an intake, contact, loading, or burden rate, then cumulative dose is:

Intensity and duration both matter. Constant exposure at intensity \(C\) over duration \(T\) produces:

When intensity varies, the integral generalizes this product:

In applied modeling, dose is rarely just arithmetic. Biological uptake, infrastructure fatigue, ecological response, human vulnerability, and institutional capacity may be nonlinear. Accumulated exposure may matter differently depending on timing, spacing, recovery, thresholds, and sequence. Two exposure histories with the same total area under the curve may not produce the same consequence if the system responds nonlinearly or has memory.

Flow-to-stock reasoning therefore distinguishes several layers:

Accumulation clarifies exposure, but interpretation still depends on the response model.

Back to top ↑

Inflows, Outflows, and Residence Time

Stocks are shaped not only by how much enters but also by how long it stays. Residence time describes how long material, energy, burden, information, or stress remains in a system before leaving, decaying, or being transformed.

A common stock-flow model with proportional removal is:

The parameter \(k\) is a removal, decay, or turnover rate. Larger \(k\) means faster removal. Smaller \(k\) means longer persistence. A rough residence time is:

This structure appears in carbon-cycle models, pollutant decay, reservoir turnover, inventory depletion, queue processing, debt repayment, biological elimination, and infrastructure recovery. It explains why systems with long residence times respond slowly to changes in inflow. Even if inflow falls, the accumulated stock may persist.

Residence time also helps interpret exposure. A short-lived substance may require sustained inflow to maintain a high stock. A persistent substance may accumulate even under modest inflow. A slow recovery process may turn repeated small stresses into substantial cumulative burden.

The balance between inflow and outflow determines whether a stock stabilizes, grows, declines, or oscillates. A steady state occurs when:

Balance is not the same as safety. A system can be in balance at an undesirable stock level. Flow-to-stock reasoning therefore asks both whether flows balance and what stock level the balance sustains.

Back to top ↑

Cumulative Burden and Lagged Consequence

Stocks create lags. A change in flow may not immediately appear as a proportional change in stock. Exposure may accumulate before consequences become visible. Infrastructure wear may build silently before failure. Health burden may emerge after latency. Ecological systems may absorb stress until thresholds are crossed.

A lagged consequence can be represented by a response function \(R(t)\) depending on accumulated exposure:

More generally, recent exposure may matter more than distant exposure. A memory kernel \(K(t-\tau)\) can weight past exposure:

This structure is important when burden decays, systems recover, or past stress remains partially active. It prevents the false assumption that all past exposure matters equally or that only present exposure matters. It also introduces interpretive responsibility: the memory kernel must be justified.

Cumulative burden is especially important in sustainability, public health, environmental justice, infrastructure management, and institutional risk. Systems can accumulate disadvantage, pollutant load, maintenance backlog, ecological stress, financial liability, or governance debt. The visible state may be the result of long histories of uneven inflow, uneven removal, and uneven exposure.

Flow-to-stock reasoning therefore supports a deeper form of accountability. It asks how past flows produced present stocks, who experienced cumulative exposure, and whether current measurements hide inherited burden.

Back to top ↑

Units and Dimensional Accounting

Flow-to-stock reasoning depends on units. If \(F(t)\) is measured in tons per year, then integrating over years produces tons. If exposure concentration is measured in micrograms per cubic meter and accumulated over hours, the result is concentration-time. If traffic load is measured as axle-load cycles per day, accumulation over days produces total load cycles.

For exposure:

Unit checks help distinguish related but different quantities:

Dimensional accounting prevents a common mistake: treating related quantities as interchangeable. A concentration is not an exposure. A flow is not a stock. A rate reduction is not necessarily a stock reduction. A cumulative burden is not necessarily a direct measure of harm unless the response model justifies that interpretation.

Back to top ↑

Measurement Windows and Baselines

Accumulation depends on the interval. A cumulative total is always tied to a start time, end time, baseline, and inclusion rule. Changing the measurement window changes the accumulated result.

Baselines matter because a stock may have a large inherited value before the modeled interval begins. Reporting only accumulation during a recent window can hide prior burden. Reporting only the inherited stock can hide current net flows. Responsible modeling should state both the initial condition and the interval accumulation when both matter.

Measurement windows also affect exposure interpretation. A daily average, annual cumulative exposure, lifetime exposure, and event-specific exposure may answer different questions. A short window may emphasize peaks. A long window may emphasize persistent burden. A moving window may reveal recent trends. A cumulative lifetime window may reveal historical inequality.

Modelers should make interval choices explicit:

Accumulation is not neutral with respect to time. The chosen window is part of the model’s claim.

Back to top ↑

Numerical Approximation and Data Resolution

Real flow and exposure data are often discrete. A model may observe daily emissions, hourly concentrations, monthly withdrawals, annual costs, or irregular sensor readings. Continuous accumulation must then be approximated from sampled data.

A simple discrete approximation is:

For irregular intervals, each \(\Delta t_i\) must be handled explicitly. For smoothly varying data, trapezoidal approximation may be more appropriate:

Data resolution matters. A daily average may miss short high-intensity peaks. Annual reporting may hide seasonal accumulation. Sparse sampling can underestimate or overestimate exposure if the measured points miss important variation. Interpolation choices can change totals.

A responsible computational workflow should record:

Numerical accumulation is not merely a calculation. It is a data interpretation pipeline.

Back to top ↑

Mathematical Deepening

This section adds a more formal layer for mathematically advanced readers. Flow-to-stock reasoning is a direct application of the Fundamental Theorem of Calculus, but systems interpretation requires attention to sign conventions, initial conditions, units, memory, discrete approximation, and the distinction between net change and total activity.

Back to top ↑

Back to top ↑

Computation and Reproducible Workflows

Computational flow-to-stock workflows should treat accumulation as an auditable process. They should record initial stock, inflows, outflows, time steps, unit conversions, integration method, baseline, measurement window, missing data rules, and interpretation of stock or exposure. They should also distinguish net change from gross activity.

A strong workflow checks whether the ending stock equals the initial stock plus accumulated net flow. It also reports cumulative inflow and cumulative outflow separately, because two systems can have the same net stock change while experiencing very different levels of activity, stress, or burden.

For exposure, computational workflows should record intensity, duration, threshold rules, population weighting, vulnerability weighting, and any memory or decay assumptions. A cumulative exposure value should never appear without its measurement window and units.

Back to top ↑

Python Workflow: Flow-to-Stock Audit

The Python workflow below calculates cumulative inflow, cumulative outflow, net accumulation, ending stock, and exposure from simple time-step records.

from __future__ import annotations

from dataclasses import dataclass, asdict
import csv
from pathlib import Path


@dataclass(frozen=True)
class FlowRecord:
    step: int
    duration: float
    inflow: float
    outflow: float
    exposure_intensity: float
    population_weight: float


@dataclass(frozen=True)
class StockExposureAudit:
    initial_stock: float
    cumulative_inflow: float
    cumulative_outflow: float
    net_accumulation: float
    ending_stock: float
    cumulative_exposure: float
    population_weighted_exposure: float
    method: str
    unit_check: str
    warning: str


def sample_records() -> list[FlowRecord]:
    return [
        FlowRecord(1, 1.0, 12.0, 6.0, 20.0, 1000.0),
        FlowRecord(2, 1.0, 10.0, 7.0, 18.0, 1100.0),
        FlowRecord(3, 1.0, 9.0, 8.0, 15.0, 1050.0),
        FlowRecord(4, 1.0, 8.0, 9.0, 13.0, 980.0),
        FlowRecord(5, 1.0, 7.0, 9.0, 11.0, 960.0),
    ]


def audit_flow_to_stock(initial_stock: float, records: list[FlowRecord]) -> StockExposureAudit:
    cumulative_inflow = sum(row.inflow * row.duration for row in records)
    cumulative_outflow = sum(row.outflow * row.duration for row in records)
    net_accumulation = cumulative_inflow - cumulative_outflow
    ending_stock = initial_stock + net_accumulation

    cumulative_exposure = sum(row.exposure_intensity * row.duration for row in records)
    population_weighted_exposure = sum(
        row.exposure_intensity * row.population_weight * row.duration
        for row in records
    )

    warnings = []
    if ending_stock < 0:
        warnings.append("ending stock is negative; check constraints or sign conventions")
    if cumulative_inflow > 0 and cumulative_outflow > 0 and abs(net_accumulation) < 0.1:
        warnings.append("large gross flows nearly cancel; report gross activity separately")

    return StockExposureAudit(
        initial_stock=initial_stock,
        cumulative_inflow=cumulative_inflow,
        cumulative_outflow=cumulative_outflow,
        net_accumulation=net_accumulation,
        ending_stock=ending_stock,
        cumulative_exposure=cumulative_exposure,
        population_weighted_exposure=population_weighted_exposure,
        method="discrete time-step accumulation",
        unit_check="flow multiplied by duration gives stock units; intensity multiplied by duration gives exposure units",
        warning="; ".join(warnings)
    )


records = sample_records()
audit = audit_flow_to_stock(50.0, records)

output_dir = Path("outputs/tables")
output_dir.mkdir(parents=True, exist_ok=True)

with (output_dir / "flow_to_stock_audit.csv").open("w", newline="", encoding="utf-8") as handle:
    writer = csv.DictWriter(handle, fieldnames=asdict(audit).keys())
    writer.writeheader()
    writer.writerow(asdict(audit))

print("Wrote flow-to-stock audit.")

This workflow separates cumulative inflow, cumulative outflow, net accumulation, ending stock, and exposure so that cumulative claims can be reviewed.

Back to top ↑

R Workflow: Exposure and Stock Diagnostics

The R workflow below performs the same stock and exposure accounting with explicit duration, flow, and population-weighted exposure records.

# Accumulation, Exposure, and Flow-to-Stock Reasoning
# Base R workflow for stock and exposure diagnostics.

records <- data.frame(
  step = 1:5,
  duration = c(1, 1, 1, 1, 1),
  inflow = c(12, 10, 9, 8, 7),
  outflow = c(6, 7, 8, 9, 9),
  exposure_intensity = c(20, 18, 15, 13, 11),
  population_weight = c(1000, 1100, 1050, 980, 960)
)

initial_stock <- 50

cumulative_inflow <- sum(records$inflow * records$duration)
cumulative_outflow <- sum(records$outflow * records$duration)
net_accumulation <- cumulative_inflow - cumulative_outflow
ending_stock <- initial_stock + net_accumulation

cumulative_exposure <- sum(records$exposure_intensity * records$duration)
population_weighted_exposure <- sum(
  records$exposure_intensity *
    records$population_weight *
    records$duration
)

warning <- ""
if (ending_stock < 0) {
  warning <- "ending stock is negative; check constraints or sign conventions"
}

audit <- data.frame(
  initial_stock = initial_stock,
  cumulative_inflow = cumulative_inflow,
  cumulative_outflow = cumulative_outflow,
  net_accumulation = net_accumulation,
  ending_stock = ending_stock,
  cumulative_exposure = cumulative_exposure,
  population_weighted_exposure = population_weighted_exposure,
  method = "discrete time-step accumulation",
  unit_check = "flow multiplied by duration gives stock units; intensity multiplied by duration gives exposure units",
  warning = warning
)

dir.create("outputs/tables", recursive = TRUE, showWarnings = FALSE)
write.csv(audit, "outputs/tables/r_flow_to_stock_audit.csv", row.names = FALSE)

print(audit)

This workflow makes the measurement window, duration, and weighting structure visible in the accumulation output.

Back to top ↑

Haskell Workflow: Typed Flow and Stock Records

Haskell can represent flow, stock, duration, and exposure with separate types, reducing the risk of mixing rates and accumulated quantities.

module Main where

newtype Duration = Duration Double deriving (Show)
newtype Flow = Flow Double deriving (Show)
newtype Stock = Stock Double deriving (Show)
newtype ExposureIntensity = ExposureIntensity Double deriving (Show)
newtype Exposure = Exposure Double deriving (Show)
newtype PopulationWeight = PopulationWeight Double deriving (Show)

data FlowRecord = FlowRecord
  { duration :: Duration
  , inflow :: Flow
  , outflow :: Flow
  , exposureIntensity :: ExposureIntensity
  , populationWeight :: PopulationWeight
  } deriving (Show)

data StockExposureAudit = StockExposureAudit
  { initialStock :: Stock
  , cumulativeInflow :: Stock
  , cumulativeOutflow :: Stock
  , netAccumulation :: Stock
  , endingStock :: Stock
  , cumulativeExposure :: Exposure
  , populationWeightedExposure :: Double
  } deriving (Show)

records :: [FlowRecord]
records =
  [ FlowRecord (Duration 1.0) (Flow 12.0) (Flow 6.0) (ExposureIntensity 20.0) (PopulationWeight 1000.0)
  , FlowRecord (Duration 1.0) (Flow 10.0) (Flow 7.0) (ExposureIntensity 18.0) (PopulationWeight 1100.0)
  , FlowRecord (Duration 1.0) (Flow 9.0)  (Flow 8.0) (ExposureIntensity 15.0) (PopulationWeight 1050.0)
  , FlowRecord (Duration 1.0) (Flow 8.0)  (Flow 9.0) (ExposureIntensity 13.0) (PopulationWeight 980.0)
  , FlowRecord (Duration 1.0) (Flow 7.0)  (Flow 9.0) (ExposureIntensity 11.0) (PopulationWeight 960.0)
  ]

stockContribution :: FlowRecord -> (Double, Double)
stockContribution row =
  let Duration dt = duration row
      Flow i = inflow row
      Flow o = outflow row
  in (i * dt, o * dt)

exposureContribution :: FlowRecord -> (Double, Double)
exposureContribution row =
  let Duration dt = duration row
      ExposureIntensity c = exposureIntensity row
      PopulationWeight p = populationWeight row
  in (c * dt, c * p * dt)

audit :: Stock -> [FlowRecord] -> StockExposureAudit
audit (Stock initial) rows =
  let inflows = map (fst . stockContribution) rows
      outflows = map (snd . stockContribution) rows
      exposures = map (fst . exposureContribution) rows
      popExposures = map (snd . exposureContribution) rows
      cumulativeIn = sum inflows
      cumulativeOut = sum outflows
      net = cumulativeIn - cumulativeOut
  in StockExposureAudit
      { initialStock = Stock initial
      , cumulativeInflow = Stock cumulativeIn
      , cumulativeOutflow = Stock cumulativeOut
      , netAccumulation = Stock net
      , endingStock = Stock (initial + net)
      , cumulativeExposure = Exposure (sum exposures)
      , populationWeightedExposure = sum popExposures
      }

main :: IO ()
main = print (audit (Stock 50.0) records)

The typed structure keeps stock units, flow units, duration, and exposure interpretation distinct.

Back to top ↑

SQL Workflow: Accumulation Assumption Registry

SQL can document flow-to-stock assumptions, especially when cumulative claims support reporting, governance, public communication, or decision review.

CREATE TABLE accumulation_assumption_registry (
    assumption_key TEXT PRIMARY KEY,
    assumption_name TEXT NOT NULL,
    mathematical_role TEXT NOT NULL,
    systems_modeling_role TEXT NOT NULL,
    review_warning TEXT NOT NULL
);

INSERT INTO accumulation_assumption_registry VALUES
(
  'initial_stock',
  'Initial stock',
  'Defines the starting value before interval accumulation begins.',
  'Keeps inherited burden separate from new accumulated change.',
  'Omitting the initial stock can make ending-stock claims misleading.'
);

INSERT INTO accumulation_assumption_registry VALUES
(
  'net_flow',
  'Net flow',
  'Defines the rate of stock change as inflow minus outflow.',
  'Explains why the stock grows, declines, or stabilizes.',
  'Sign conventions must be explicit.'
);

INSERT INTO accumulation_assumption_registry VALUES
(
  'gross_flows',
  'Gross flows',
  'Records cumulative inflow and cumulative outflow separately.',
  'Prevents large offsetting activity from being hidden by net change.',
  'Net change alone can conceal high stress or turnover.'
);

INSERT INTO accumulation_assumption_registry VALUES
(
  'exposure_window',
  'Exposure window',
  'Defines the interval over which exposure is accumulated.',
  'Clarifies whether the result is event-based, annual, lifetime, or moving-window exposure.',
  'Changing the window changes the cumulative claim.'
);

INSERT INTO accumulation_assumption_registry VALUES
(
  'unit_consistency',
  'Unit consistency',
  'Checks that rate multiplied by time produces stock or exposure units.',
  'Prevents flow, stock, concentration, and exposure from being confused.',
  'Unit mismatch indicates an invalid cumulative interpretation.'
);

SELECT
    assumption_name,
    mathematical_role,
    systems_modeling_role,
    review_warning
FROM accumulation_assumption_registry
ORDER BY assumption_key;

This registry keeps cumulative interpretation tied to initial conditions, net flows, gross activity, exposure windows, and unit consistency.

Back to top ↑

GitHub Repository

The companion repository for this article is designed as a reproducible mathematical-modeling workspace. It supports flow-to-stock audits, inflow-outflow accounting, cumulative exposure diagnostics, population-weighted exposure, unit checks, measurement-window documentation, typed flow and stock records, SQL assumption registries, generated outputs, and advanced mathematical audit reports.

Back to top ↑

Interpretive Limits and Responsible Use

Accumulation is powerful because it connects local rates to cumulative consequence. It is also easy to misuse. A cumulative total can be made larger or smaller by changing the measurement window. A net stock change can hide large gross flows. A rate reduction can be mistaken for stock reduction. A cumulative exposure measure can be treated as harm without a justified response model. A stock can inherit burden from before the reporting interval.

Responsible use requires several checks. State the stock. Define inflows and outflows. Specify sign conventions. Record the initial condition. Report cumulative inflow, cumulative outflow, and net accumulation separately when possible. Define exposure intensity, duration, pathway, population, and weighting. State the measurement window and baseline. Check units. Document numerical method and data resolution. Distinguish accumulated exposure from modeled consequence.

The central modeling question is not only “What accumulated?” It is “What flowed, over what interval, into what stock or exposure measure, under what assumptions, with what units, and with what consequences for interpreting cumulative burden?”

Back to top ↑

Related Articles

• Calculus for Systems Modeling
• Improper Integrals and Unbounded Quantities
• Integration by Parts and Structured Decomposition
• Definite Integrals and Total Change
• The Fundamental Theorem of Calculus
• Antiderivatives and the Recovery of Accumulation
• Stocks, Flows, and Accumulation
• Carbon Accumulation and Emissions Pathways
• Environmental Systems Modeling
• Public Health Modeling

Back to top ↑

Further Reading

• Apostol, T.M. (1967) Calculus, Volume 1: One-Variable Calculus, with an Introduction to Linear Algebra. 2nd edn. New York: Wiley.
• Spivak, M. (2008) Calculus. 4th edn. Houston, TX: Publish or Perish.
• Courant, R. and John, F. (1999) Introduction to Calculus and Analysis, Volume I. Berlin: Springer.
• Abbott, S. (2015) Understanding Analysis. 2nd edn. New York: Springer.
• Sterman, J.D. (2000) Business Dynamics: Systems Thinking and Modeling for a Complex World. Boston, MA: Irwin/McGraw-Hill.
• Meadows, D.H. (2008) Thinking in Systems: A Primer. White River Junction, VT: Chelsea Green Publishing.
• Ford, A. (2010) Modeling the Environment. 2nd edn. Washington, DC: Island Press.
• Forrester, J.W. (1961) Industrial Dynamics. Cambridge, MA: MIT Press.
• Massachusetts Institute of Technology (MIT) OpenCourseWare (2010) Single Variable Calculus. Cambridge, MA: MIT OpenCourseWare.
• OpenStax (2016a) Calculus Volume 1. Houston, TX: OpenStax, Rice University.

Back to top ↑

References

• Abbott, S. (2015) Understanding Analysis. 2nd edn. New York: Springer.
• Apostol, T.M. (1967) Calculus, Volume 1: One-Variable Calculus, with an Introduction to Linear Algebra. 2nd edn. New York: Wiley.
• Courant, R. and John, F. (1999) Introduction to Calculus and Analysis, Volume I. Berlin: Springer.
• Ford, A. (2010) Modeling the Environment. 2nd edn. Washington, DC: Island Press.
• Forrester, J.W. (1961) Industrial Dynamics. Cambridge, MA: MIT Press.
• Massachusetts Institute of Technology (MIT) OpenCourseWare (2010) Single Variable Calculus. Cambridge, MA: MIT OpenCourseWare.
• Meadows, D.H. (2008) Thinking in Systems: A Primer. White River Junction, VT: Chelsea Green Publishing.
• OpenStax (2016a) Calculus Volume 1. Houston, TX: OpenStax, Rice University.
• Spivak, M. (2008) Calculus. 4th edn. Houston, TX: Publish or Perish.
• Sterman, J.D. (2000) Business Dynamics: Systems Thinking and Modeling for a Complex World. Boston, MA: Irwin/McGraw-Hill.

Back to top ↑