Evaluate one arithmetic/scientific expression safely. Use explicit * and variables. For dimensional work, pass original unit-bearing variable strings (for example 12.4 kN*m) and always set output_unit; never add manual conversion factors to the expression. Returns exact and high-precision numeric forms when unitless.

Convert a dimensionful value such as 32 degF, 120 km/hour, or 5 psi to a compatible target unit. Dimensional mismatches are errors rather than guessed conversions.

Solve up to eight algebraic equations for explicitly named variables. Each equation may use left = right or a zero-valued expression. Returns every solution in exact and numerical form; choose the real or complex domain explicitly.

Differentiate, integrate, take a limit, or expand a series symbolically. For definite integrals provide both bounds; for limits and series provide point. order selects derivative order or the number of series terms.

Perform exact matrix and linear-algebra operations on matrices up to 25x25. Supply other_matrix for multiplication or as the right-hand side of solve. String entries such as 1/3 preserve exact values. Condition-number results include numerical-stability warnings.

Compute descriptive statistics, Pearson correlation, linear regression, a Student-t confidence interval for a mean, or an independent Welch t-test. Provide data_y for paired correlation/regression or the second t-test sample.

Propagate independent standard uncertainties through an expression. Each measurement has value, uncertainty, and optional unit. Returns nominal value, standard uncertainty, relative uncertainty, and a coverage interval; use output_unit when the result is dimensionful.

Look up a physical constant from SciPy's bundled CODATA table by common alias or name. Returns its value, unit, standard uncertainty, source metadata, and an optional compatible unit conversion. Examples: c, G, hbar, Boltzmann constant, electron mass.