function out1 = sdl_carem_parameter_to_coefficient_throttle(S_leak,D_plate,theta_rest)
t2 = cos(theta_rest);
t3 = D_plate.^2;
t4 = pi.^2;
t5 = theta_rest.^2;
t6 = theta_rest.^3;
t8 = theta_rest.*1.8e+2;
t10 = theta_rest-9.0e+1;
t11 = (theta_rest.*pi)./1.8e+2;
t7 = 1.0./t2;
t9 = -t8;
t12 = 1.0./t10;
t13 = cos(t11);
t14 = sin(t11);
t15 = t5+t9+8.1e+3;
t17 = t3.*t5.*t14.*pi;
t18 = t3.*t13.*theta_rest.*5.4e+2;
t16 = 1.0./t15;
out1 = [t7.*t12.*t16.*(S_leak.*t2.*2.916e+6+t2.*t3.*pi.*7.29e+5-t3.*t13.*pi.*7.29e+5...
        +S_leak.*t2.*t5.*1.08e+3-S_leak.*t2.*t6.*4.0-S_leak.*t2.*theta_rest.*9.72e+4...
        +t2.*t3.*t5.*pi.*2.7e+2-t2.*t3.*t6.*pi-t2.*t3.*theta_rest.*pi.*2.43e+4...
        +t3.*t13.*theta_rest.*pi.*2.43e+4+t3.*t4.*t5.*t14.*4.5e+1...
        -t3.*t4.*t14.*theta_rest.*4.05e+3).*(-1.0./4.0),...
    t7.*t12.*t16.*(t17+t18-t3.*t14.*pi.*4.05e+3-t3.*t14.*theta_rest.*pi.*4.5e+1).*4.5e+1,...
    (t7.*t12.*(t17+t18+t3.*t13.*4.86e+4-t3.*t14.*pi.*1.62e+4+...
        t3.*t14.*theta_rest.*pi.*9.0e+1).*-4.5e+1)./(pi.*8.1e+3+t5.*pi-theta_rest.*pi.*1.8e+2),...
    (t7.*t12.*(t3.*t13.*3.6e+2-t3.*t14.*pi.*9.0e+1+...
        t3.*t14.*theta_rest.*pi).*8.1e+3)./(t4.*8.1e+3+t4.*t5-t4.*theta_rest.*1.8e+2)];
