HW4_Script.m

%HW_4 Script

setdefaults;

%Part B: using the goldmin function
theta = 0;
[x1,fx1,ea1,iter1] = goldmin(@(x) collar_potential_energy(x,theta),-10,10,0.0001);

%Part C: Creating For-loop solving for potential energy

theta = .9*pi/180;  
thetaVector = linspace(0,pi, 400);
fx2 = 0;

[x3,fx3,ea3,iter3] = goldmin(@(x) collar_potential_energy(x,theta),-10,10,0.0001);

for n = 1:length(thetaVector)
    theta = thetaVector(n);
[x2,fx2(n),ea2,iter2] = goldmin(@(x) collar_potential_energy(x,theta),-10,10,0.0001);
end

%Part D: Plotting the Data

plot(thetaVector, fx2);
    title('Plot of Minimum PE for x_C with given Angle');
    xlabel('Theta (radians)');
    ylabel('x_C (meters)');