四杆机构运动分析

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（