Analisis de posicion

%Análisis de posición

        E2=input('E2: ');

        E3=input('E3: ');

        E4=input('E4: ');

        theta2=input('theta2:');

        vel.ang2=input('vel.ang2:');

        alfa2=input('alfa2');

         %Condicion de grashof

        if ((E2>E3)&(E2>E4))

            s=E2;

        elseif ((E2

            l=E2;

        else

            p=E2;

        end

        if ((E3>E2)&(E3>E4))

            s=E3;

        elseif ((E3

            l=E3;

        else

            p=E3;

        end

        if ((E4>E2)&(E4>E3))

            s=E4;

        elseif ((E4

            l=E4;

        else

            P=E4;

        end

        if   s+l<=p+p

            disp('se cumple la condicion de grashof');

        theta3_1=asind((E2*sind(theta2)-E4)/(E3));

        E1=E2*cosd(theta2-E4*cosd(theta3_1));

        fprintf('Las posiciones en abierto son \n');

        fprintf('theta3_1: %f \n', theta3_1);

        fprintf('E1: %f \n\n', E1);

        theta3_2=asind((-asind(theta2)-E4)/(E3))+pi;

        E1=E2*cosd(theta2-E4*cosd(theta3_2));

        fprintf('Las posiciones en cerrado son ');

        fprintf('theta3_2: %f \n', theta3_2);

        fprintf('E1: %f \n\n' , E1);

        %Análisis de velocidad angular

        W3_1=((E2*cosd(theta2))/(E3*cosd(theta3_1)))*vel.ang2;

        E1=((-E2*vel.ang2*sind(theta2))+(E3*W3_1*sind(theta3_1)));

        fprintf('Las velocidades angulares en abierto son \n ');

        fprintf('W3_1: %f \n', W3_1);

        fprintf('E1: %f \n\n', E1);

...