1 July 29,2015 2 P4 3 syms p1 p2 p3 4 phi(:,1) = 1/6*p1.*(4*p1-1).*(4*p1-2).*(4*p1-3); 5 phi(:,2) = 1/6*p2.*(4*p2-1).*(4*p2-2).*(4*p2-3); 6 phi(:,3) = 1/6*p3.*(4*p3-1).*(4*p3-2).*(4*p3-3); 7 phi(:,4) = 32*p1.*p2.*p3.*(4*p1-1); 8 phi(:,5) = 32*p1.*p2.*p3.*(4*p2-1); 9 phi(:,6) = 32*p1.*p2.*p3.*(4*p3-1); 10 phi(:,7) = 8/3*p1.*p2.*(4*p1-1).*(4*p1-2); 11 phi(:,8) = 4*p1.*p2.*(4*p1-1).*(4*p2-1); 12 phi(:,9) = 8/3*p1.*p2.*(4*p2-1).*(4*p2-2); 13 phi(:,10) = 8/3*p1.*p3*(4*p1-1).*(4*p1-2); 14 phi(:,11) = 4*p1.*p3.*(4*p1-1).*(4*p3-1); 15 phi(:,12) = 8/3*p1.*p3.*(4*p3-1).*(4*p3-1); 16 phi(:,13) = 8/3*p2.*p3.*(4*p2-1).*(4*p2-2); 17 phi(:,14) = 4*p2.*p3.*(4*p2-1).*(4*p3-1); 18 phi(:,15) = 8/3*p2.*p3.*(4*p3-1).*(4*p3-2);
1 July 29,2015 2 P4 3 syms p1 p2 p3 4 phi(:,1) = 1/6*p1.*(4*p1-1).*(4*p1-2).*(4*p1-3); 5 phi(:,2) = 1/6*p2.*(4*p2-1).*(4*p2-2).*(4*p2-3); 6 phi(:,3) = 1/6*p3.*(4*p3-1).*(4*p3-2).*(4*p3-3); 7 phi(:,4) = 8/3*p2.*p3.*(4*p2-1).*(4*p2-2); 8 phi(:,5) = 8/3*p1.*p3.*(4*p3-1).*(4*p3-1); 9 phi(:,6) = 8/3*p1.*p2.*(4*p1-1).*(4*p1-2); 10 phi(:,7) = 4*p2.*p3.*(4*p2-1).*(4*p3-1); 11 phi(:,8) = 4*p1.*p3.*(4*p1-1).*(4*p3-1); 12 phi(:,9) = 4*p1.*p2.*(4*p1-1).*(4*p2-1); 13 phi(:,10) = 8/3*p2.*p3.*(4*p3-1).*(4*p3-2); 14 phi(:,11) = 8/3*p1.*p3.*(4*p1-1).*(4*p1-2); 15 phi(:,12) = 8/3*p1.*p2.*(4*p2-1).*(4*p2-2); 16 phi(:,13) = 32*p1.*p2.*p3.*(4*p1-1); 17 phi(:,14) = 32*p1.*p2.*p3.*(4*p2-1); 18 phi(:,15) = 32*p1.*p2.*p3.*(4*p3-1);