l1=input（‘enter‘‘l1‘‘:‘）;
l2=input（‘enter‘‘l2‘‘:‘）;
l3=input（‘enter‘‘l3‘‘:‘）;
l4=input（‘enter‘‘l4‘‘:‘）;
w1=input（‘enter‘‘w1‘‘:‘）;
d=pi/180;
for?j=1:360
????m=j*d;
????Y（j:）=fsolve（@（b）[l2*cos（b（1）*d）-l3*cos（b（2）*d）+l1*cos（m）-l4;
????l2*sin（b（1）*d）-l3*sin（b（2）*d）+l1*sin（m）][20;30]‘options‘）;
????y（j:）=Y（j:）*d;
????w（:j）=（[-l2*sin（y（j1））?l3*sin（y（j2））;
????l2*cos（y（j1））?-l3*cos（y（j2））]）\[l1*sin（m）;l1*cos（m）]*w1;
????C=[-l2*sin（y（j1））?l3*sin（y（j2））;
????????l2*cos（y（j1））?-l3*cos（y（j2））];
????D=w1*[w1*l1*cos（m）;
????????w1*l1*sin（m）]-[-w（