

C_calc=S_forw_normal*material_distribution(:,1);
Cn=(1./C_calc-1./C_low)./(1./C_high-1./C_low);
pic=S_inv_normal*Cn;
trisurf(simp,vert(:,1),vert(:,2),vert(:,3),pic);

C_high=S_forw_normal*(3*ones(size(matS)));
Cn=(1./C_calc-1./C_low)./(1./C_high-1./C_low);
pic=S_inv_normal*Cn;
trisurf(simp,vert(:,1),vert(:,2),vert(:,3),pic);
C_high=S_forw_normal*(4*ones(size(matS)));
Cn=(1./C_calc-1./C_low)./(1./C_high-1./C_low);
pic=S_inv_normal*Cn;
trisurf(simp,vert(:,1),vert(:,2),vert(:,3),pic);
trisurf(simp,vert(:,1),vert(:,2),vert(:,3),matS_2);
trisurf(simp,vert(:,1),vert(:,2),vert(:,3),matS_2./2);