function [ model_handle ] = PlotImplicitQuadratcSurface( C,Domain,step )
color='r';
xm= linspace(Domain(1),Domain(2),step);
ym = linspace(Domain(3),Domain(4),step);
zm = linspace(Domain(5),Domain(6),step);
[X,Y,Z] = meshgrid(xm,ym,zm);
fvalues=C(1)+C(2).*X+C(3).*Y+C(4).*Z+C(5).*X.^2+C(6).*Y.^2+C(7).*Z.^2+C(8).*X.*Y+C(9).*Y.*Z+C(10).*Z.*X;
model_handle = patch(isosurface(X,Y,Z,fvalues,0));
isonormals(X,Y,Z,fvalues,model_handle)
set(model_handle,'FaceColor',color,'EdgeColor','none');
alpha(0.7)
grid on;
view([1,1,1]);
axis equal;
camlight;
lighting gouraud
end
0 Comments
Sign in to comment.