function hermit2(Pvec,Rvec,TWvec)
%HERMIT Summary of this function goes here
%   Detailed explanation goes here
    Gx = [Pvec(1,1) Pvec(2,1) Rvec(2,1) Rvec(4,1);
        Pvec(3,1) Pvec(4,1) Rvec(6,1) Rvec(8,1);
        Rvec(1,1) Rvec(3,1) TWvec(1,1) TWvec(2,1);
        Rvec(5,1) Rvec(7,1) TWvec(3,1) TWvec(4,1)];
    Gy = [Pvec(1,2) Pvec(2,2) Rvec(2,2) Rvec(4,2);
        Pvec(3,2) Pvec(4,2) Rvec(6,2) Rvec(8,2);
        Rvec(1,2) Rvec(3,2) TWvec(1,2) TWvec(2,2);
        Rvec(5,2) Rvec(7,2) TWvec(3,2) TWvec(4,2)];
    Gz = [Pvec(1,3) Pvec(2,3) Rvec(2,3) Rvec(4,3);
        Pvec(3,3) Pvec(4,3) Rvec(6,3) Rvec(8,3);
        Rvec(1,3) Rvec(3,3) TWvec(1,3) TWvec(2,3);
        Rvec(5,3) Rvec(7,3) TWvec(3,3) TWvec(4,3)];
    s = 0:0.1:1;
    t = 0:0.1:1;
    S = [s.*s.*s;s.*s;s;ones(size(s))];
    T = [t.*t.*t;t.*t;t;ones(size(t))];
    M = [2 -3 0 1;-2 3 0 0;1 -2 1 0;1 -1 0 0];
    X = S'*M'*Gx*M*T;
    Y = S'*M'*Gy*M*T;
    Z = S'*M'*Gz*M*T;
    mesh(X,Y,Z)
end