clear allclcTrebuie sa adaugi 0.05*i la l0, l1, l2, l3, BS2, CS3, unde i=nr din condical0=0.9;l1=0.35;l2=0.7;l3=0.6;BS2=0.4;CS3=0.2;teta=pi/9;phi1=pi/3;omega1=5;epsilon1=2;A2=2*l2*(l1*cos(phi1)-l0)B2=2*l1*l2*sin(phi1)C2=l0*l0+l1*l1+l2*l2-l3*l3-2*l0*l1*cos(phi1)A3=2*l3*(l1*cos(phi1)-l0)B3=2*l1*l3*sin(phi1)C3=l0*l0+l1*l1+l3*l3-l2*l2-2*l0*l1*cos(phi1)t1=(B2-sqrt(A2^2+B2^2-C2^2))/(A2-C2)t2=(B2+sqrt(A2^2+B2^2-C2^2))/(A2-C2)phi2=2*atan(t1)phi2p=2*atan(t2)t3=(B3+sqrt(A3^2+B3^2-C3^2))/(A3-C3)t4=(B3-sqrt(A3^2+B3^2-C3^2))/(A3-C3)phi3=2*atan(t3)phi3p=2*atan(t4)omega2=omega1*(l1/l2)*(l2*sin(phi1-phi2)-l0*sin(phi1))/(l1*sin(phi1-phi2)+l0*sin(phi2))omega3=omega1*(l1/l3)*(l3*sin(phi1-phi3)-l0*sin(phi1))/(l1*sin(phi1-phi3)+l0*sin(phi3))i21=omega2/omega1;i31=omega3/omega1;epsilon2=(epsilon1*(l1/l2)*(l2*sin(phi1-phi2)-l0*sin(phi1))/(l1*sin(phi1-phi2)+l0*sin(phi2)))+((l1*l2*(omega1-omega2)^2)*cos(phi1-phi2)-(omega2^2)*l0*l2*cos(phi2)-(omega1^2)*l0*l1*cos(phi1))/(l1*l2*sin(phi1-phi2)+l0*l2*sin(phi2))epsilon3=(epsilon1*(l1/l3)*(l2*sin(phi1-phi3)-l0*sin(phi1))/(l1*sin(phi1-phi3)+l0*sin(phi3)))+((l1*l3*(omega1-omega3)^2)*cos(phi1-phi3)-(omega3^2)*l0*l3*cos(phi3)-(omega1^2)*l0*l1*cos(phi1))/(l1*l3*sin(phi1-phi3)+l0*l3*sin(phi3))xB=l1*cos(phi1)yB=l1*sin(phi1)xS2=l1*cos(phi1)+BS2*cos(phi2+teta)yS2=l1*sin(phi1)+BS2*sin(phi2+teta)xC=l1*cos(phi1)+l2*cos(phi2)yC=l1*sin(phi1)+l2*sin(phi2)xS3=xC+CS3*cos(phi3)yS3=yC+CS3*sin(phi3)vxB=-omega1*l1*sin(phi1);vyB=omega1*l1*cos(phi1);vB=sqrt(vxB^2+vyB^2)vxS2=-omega1*l1*sin(phi1)-omega2*BS2*sin(phi2+teta);vyS2=omega1*l1*cos(phi1)+omega2*BS2*cos(phi2+teta);vS2=sqrt(vxS2^2+vyS2^2)vxC=-omega1*l1*sin(phi1)-omega2*l2*sin(phi2);vyC=omega1*l1*cos(phi1)+omega2*l2*cos(phi2);vC=sqrt(vxC^2+vyC^2)vxS3=vxC-omega3*CS3*sin(phi3);vyS3=vyC+omega3*CS3*cos(phi3);vS3=sqrt(vxS3^2+vyS3^2)axB=-epsilon1*l1*sin(phi1)-(omega1^2)*l1*cos(phi1);ayB=epsilon1*l1*cos(phi1)-(omega1^2)*l1*sin(phi1);aB=sqrt(axB^2+ayB^2)axS2=-epsilon1*l1*sin(phi1)-(omega1^2)*l1*cos(phi1)-epsilon2*BS2*sin(phi2+teta)-(omega2^2)*BS2*cos(phi2+teta);ayS2=epsilon1*l1*cos(phi1)-(omega1^2)*l1*sin(phi1)+epsilon2*BS2*cos(phi2+teta)-(omega2^2)*BS2*sin(phi2+teta);aS2=sqrt(axS2^2+ayS2^2)axC=-epsilon1*l1*sin(phi1)-(omega1^2)*l1*cos(phi1)-epsilon2*l2*sin(phi2)-(omega2^2)*l2*cos(phi2);ayC=epsilon1*l1*cos(phi1)-(omega1^2)*l1*sin(phi1)+epsilon2*l2*cos(phi2)-(omega2^2)*l2*sin(phi2);aC=sqrt(axC^2+ayC^2)axS3=axC-epsilon3*CS3*sin(phi3)-(omega3^2)*CS3*cos(phi3);ayS3=ayC+epsilon3*CS3*cos(phi3)-(omega3^2)*CS3*sin(phi3);aS3=sqrt(axS3^2+ayS3^2)