Would this solve your problem?

#include <concepts>
#include <cmath>
#include <utility>
#include <tuple>

double funct(std::floating_point auto... ts){
    const auto tuple = std::make_tuple(std::forward<decltype(ts)>(ts)...);
    const double z = (std::get<0>(tuple) - std::get<2>(tuple)) / std::get<3>(tuple);
    return std::get<1>(tuple) * std::exp(-0.5 * z * z);
}

int main()
{
    funct(0.0, 0.0, 0.0, 0.0);
    return 0;
}