function [g dg d2g] = masterp1(id,x) n=size(x,2); for i=1:n if (id == 1) g(1,i) = 0.5*(1-x(1,i)); g(2,i) = 1-g(1,i); dg(1,i) = -0.5; dg(2,i) = 0.5; d2g(1,i) = 0; d2g(2,i) = 0; elseif (id == 2) g(1,i) = 1-x(1,i)-x(2,i); g(2,i) = x(1,i); g(3,i) = x(2,i); dg(1,1,i) = -1; dg(1,2,i) = -1; dg(2,1,i) = 1; dg(2,2,i) = 0; dg(3,1,i) = 0; dg(3,2,i) = 1; d2g(1:3,1:2,i)=zeros(3,2); end end end