for k = 1:101
    x(k) = (k-1)*0.01*2*pi;
    u0(k) = sin(x(k));
    z(k) = 0;
end
t = 0.3;
for k = 1:101
    x(k) = (k-1)*0.01*2*pi;
    a = 1/t;
    b = x(k)/t;
    xl = 0;
    xr = 2*pi;
    xm = (xl+xr)/2;
    for ip = 1:10
        if ((b-a*xl-sin(xl))*(b-a*xm-sin(xm)) <= 0)
            xr = xm;
        else
            xl = xm;
        end
        xm = (xl+xr)/2;
    end
    x0 = xm;
    u1(k) = sin(x0);
end
t = 0.7;
for k = 1:101
    x(k) = (k-1)*0.01*2*pi;
    a = 1/t;
    b = x(k)/t;
    xl = 0;
    xr = 2*pi;
    xm = (xl+xr)/2;
    for ip = 1:10
        if ((b-a*xl-sin(xl))*(b-a*xm-sin(xm)) <= 0)
            xr = xm;
        else
            xl = xm;
        end
        xm = (xl+xr)/2;
    end
    x0 = xm;
    u2(k) = sin(x0);
end
plot(x,u0,'k-',x,u1,'b-',x,u2,'r-',x,z,'k--')
axis([0 2*pi -1.5 1.5])