M = [2 -3 0 1;-2 3 0 0;1 -2 1 0;1 -1 0 0];


p00=[0 0 0];
p01=[0 1 0];
p10=[1 0 0];
p11=[1 1 0];
r000=[1 1 1];
r001=[1 1 1];
r010=[1 -1 1];
r011=[-1 1 -1];
r100=[1 -1 -1];
r101=[-1 1 1];
r110=[1 1 -1];
r111=[1 1 -1];

TW00=[0 0 0];
TW01=[0 0 0];
TW10=[0 0 0];
TW11=[0 0 0];

Gx=[[p00(1) p01(1) r001(1) r011(1)];
    [p10(1) p11(1) r101(1) r111(1)];
    [r000(1) r010(1) TW00(1) TW01(1)];
    [r100(1) r110(1) TW10(1) TW11(1)]];

Gy=[[p00(2) p01(2) r001(2) r011(2)];
    [p10(2) p11(2) r101(2) r111(2)];
    [r000(2) r010(2) TW00(2) TW01(2)];
    [r100(2) r110(2) TW10(2) TW11(2)]];

Gz=[[p00(3) p01(3) r001(3) r011(3)];
    [p10(3) p11(3) r101(3) r111(3)];
    [r000(3) r010(3) TW00(3) TW01(3)];
    [r100(3) r110(3) TW10(3) TW11(3)]];
t= 0:0.01:1;
s= 0:0.01:1;

T = [t.^3;t.^2;t;ones(size(t))] ;

S = [s.^3;s.^2;s;ones(size(s))] ;

Qx = (S.')*(M.')*Gx*M*T;
Qy = (S.')*(M.')*Gy*M*T;
Qz = (S.')*(M.')*Gz*M*T;

mesh(Qx,Qy,Qz)



%zad 1
hermit([0,0,0], [0,1,0], [1,0,0], [1,1,0], [1,1,1], [1,1,1], [1,-1,1],[-1,1,-1],[1,-1,-1],[-1,1,1],[1,1,-1],[1,1,-1],[0,0,0],[0,0,0],[0,0,0],[0,0,0]);

P00=[0,0,0];P01=[0,1,0];P10=[1,0,0];P11=[1,1,0],R000=[1,0,1];R001=[0,1,1],R010=[1,0,1];R011=[0,1,-1];R100=[1,0,-1];R101=[0,1,1];R110=[1,0,-1];R111=[0,1,-1]
TW = [0,0,0];
hermit(P00,P01,P10,P11,R000,R001,R010,R011,R100,R101,R110,R111,TW,TW,TW,TW)
TW = [1,2,3];
hermit(P00,P01,P10,P11,R000,R001,R010,R011,R100,R101,R110,R111,TW,TW,TW,TW)