FriesischScott

using Distributions

function P(x::Float64, d::Int)
    if d == 0
        return 1.0
    end
    if d == 1
        return x
    end

FriesischScott

using Distributions
using FastGaussQuadrature

function P(x::Float64, d::Int)
    if d == 0
        return 1.0
    end

    if d == 1
        return x

FriesischScott

using Distributions
using SpecialFunctions
using Plots

"""
Generator
"""
function φ(x::Real, ϑ::Real)
    return (1+x)^(−1/ϑ)
end

FriesischScott

function carbonation(P, T) {
    return (P + 1.013) * Math.exp(-10.73797 + (2617.25 / (T + 273.15))) * 10
}