create a loop to solve for all theta values in a 4 bar mechanism in matlab

main.m
% Define range of theta values to iterate through
theta_range = 0:pi/180:2*pi;

% Preallocate arrays to store end effector positions for each set of joint angles
x_end = zeros(length(theta_range)^4,1);
y_end = zeros(length(theta_range)^4,1);

% Loop through all combinations of joint angles
count = 1;
for i = 1:length(theta_range)
    for j = 1:length(theta_range)
        for k = 1:length(theta_range)
            for l = 1:length(theta_range)
                
                % Calculate end effector position and store in arrays
                [x_end(count), y_end(count)] = four_bar_fwd(theta_range(i),theta_range(j),theta_range(k),theta_range(l));
                
                count = count + 1;
            end
        end
    end
end

% Plot end effector positions
figure;
plot(x_end,y_end,'.');
xlabel('X position');
ylabel('Y position');
title('End Effector Positions for all Joint Angle Combinations');

% Define function for forward kinematics of 4 bar mechanism
function [x_end, y_end] = four_bar_fwd(theta1,theta2,theta3,theta4)
% Define mechanism parameters
a = 1;
b = 2;
c = 3;
d = 2.5;

% Calculate end effector position using forward kinematics equations
x_end = a*cos(theta1) + b*cos(theta2) - c*cos(theta3) - d*cos(theta4);
y_end = a*sin(theta1) + b*sin(theta2) - c*sin(theta3) - d*sin(theta4);
1320 chars
42 lines