clear;
clf;

n = 1;

LEPx = [0 0 0];
LEPy = [0 0 0];

LGPx = [0.25 0.5 3];
LGPy = [0.25 0.5 3];

RGPx = [0.75  2  2];
RGPy = [0.25 -1 -1];

REPx = [1 1 1];
REPy = [0 0 0];


clrs = 'rgbmck';

for i = 1:n
    a0 = LEPx(i);
    a1 = 3 * (LGPx(i) - LEPx(i));
    a2 = 3 * (LEPx(i) + RGPx(i) - 2*LGPx(i));
    a3 = REPx(i) - LEPx(i) + 3 * (LGPx(i) - RGPx(i));

    b0 = LEPy(i);
    b1 = 3 * (LGPy(i) - LEPy(i));
    b2 = 3 * (LEPy(i) + RGPy(i) - 2*LGPy(i));
    b3 = REPy(i) - LEPy(i) + 3 * (LGPy(i) - RGPy(i));

    tt = 0 : 0.01 : 1;
    for j = 1 : length(tt)
        t = tt(j);
        x(j) = a0 + a1*t + a2*t^2 + a3*t^3;
        y(j) = b0 + b1*t + b2*t^2 + b3*t^3;
    end

    hold on;
    plot(x, y, clrs(i));
    plot(LEPx(i), LEPy(i), strcat(clrs(i),'o'));
    plot(REPx(i), REPy(i), strcat(clrs(i),'o'));
    plot(LGPx(i), LGPy(i), strcat(clrs(i),'*'));    
    plot(RGPx(i), RGPy(i), strcat(clrs(i),'*'));    
end