function [bc_params] = gen_bezier(P, S, theta1, theta2, d1, d2)

    bc_params = [P (P+d1*[cos(theta1);sin(theta1)]) (S-d2*[cos(theta2);sin(theta2)]) S];
end
