s=1
k=5
e=9
sy=1
ky=5
ey=3

b=[sy ky ky ey 0 0 0 0]';

A=[s^3 s^2 s 1 0 0 0 0
   k^3 k^2 k 1 0 0 0 0
   0 0 0 0 k^3 k^2 k 1
   0 0 0 0 e^3 e^2 e 1
   3*k^2 2*k 1 0 -3*k^2 -2*k -1 0
   6*k 2 0 0 -6*k -2 0 0
   6*s 2 0 0 0 0 0 0
   0 0 0 0 6*e 2 0 0]
v=A\b
s1=v(1)*s^3+v(2)*s^2+v(3)*s+v(4)
s1=v(1)*k^3+v(2)*k^2+v(3)*k+v(4)

s2=v(5)*k^3+v(6)*k^2+v(7)*k+v(8)
s2=v(5)*e^3+v(6)*e^2+v(7)*e+v(8)
%%%%%%%
s=1
k=5
e=9
f=12
sy=1
ky=5
ey=3
fy=7

b=[sy ky ky ey ey fy 0 0 0 0 0 0]';

A=[s^3 s^2 s 1 0 0 0 0 0 0 0 0
   k^3 k^2 k 1 0 0 0 0 0 0 0 0
   0 0 0 0 k^3 k^2 k 1 0 0 0 0
   0 0 0 0 e^3 e^2 e 1 0 0 0 0
   0 0 0 0 0 0 0 0 e^3 e^2 e 1
   0 0 0 0 0 0 0 0 f^3 f^2 f 1
   3*k^2 2*k 1 0 -3*k^2 -2*k -1 0 0 0 0 0
   0 0 0 0 3*e^2 2*e 1 0 -3*e^2 -2*e -1 0
   6*k 2 0 0 -6*k -2 0 0 0 0 0 0
   0 0 0 0 6*e 2 0 0 -6*e -2 0 0
   6*s 2 0 0 0 0 0 0 0 0 0 0
   0 0 0 0 0 0 0 0 6*f 2 0 0]
v=A\b
s1=v(1)*s^3+v(2)*s^2+v(3)*s+v(4)
s1=v(1)*k^3+v(2)*k^2+v(3)*k+v(4)

s2=v(5)*k^3+v(6)*k^2+v(7)*k+v(8)
s2=v(5)*e^3+v(6)*e^2+v(7)*e+v(8)


s3=v(9)*e^3+v(10)*e^2+v(11)*e+v(12)
s3=v(9)*f^3+v(10)*f^2+v(11)*f+v(12)

x=[s k e f]
y=[sy ky ey fy]

plot(x,y,'*')
hold on
plot(0,0)
plot(15,8,'w.')
