function out1 = sdl_carem_parameter_to_coefficient_volumetric(eta_v_0,x_eta_peak,eta_v_peak,x_max,eta_v_xmax)
t2 = x_eta_peak.^2;
t3 = x_eta_peak.^3;
t4 = x_max.^2;
t5 = x_max.^3;
t6 = 1.0./t2;
t7 = t2.*x_max;
t8 = t4.*x_eta_peak.*2.0;
t9 = -t8;
t10 = t5+t7+t9;
t11 = 1.0./t10;
out1 = [eta_v_0,...
    -(t11.*(eta_v_0.*t3+eta_v_0.*t5.*2.0-eta_v_peak.*t5.*2.0-eta_v_xmax.*t3...
        -eta_v_0.*t4.*x_eta_peak.*3.0+eta_v_peak.*t4.*x_eta_peak.*3.0))./x_eta_peak,...
    t6.*t11.*(eta_v_0.*t3.*2.0+eta_v_0.*t5-eta_v_0.*t7.*3.0-eta_v_peak.*t5...
        +eta_v_peak.*t7.*3.0-eta_v_xmax.*t3.*2.0),...
    -t6.*t11.*(eta_v_0.*t2+eta_v_0.*t4-eta_v_peak.*t4-eta_v_xmax.*t2...
        -eta_v_0.*x_eta_peak.*x_max.*2.0+eta_v_peak.*x_eta_peak.*x_max.*2.0)];
