%?%單純3D心
%?clear
%?clc
%?[xyz]=meshgrid（linspace（-1.31.3））;
%?val=（x.^2?+?（9/4）*y.^2?+?z.^2?-?1）.^3?-?x.^2.*z.^3?-?（9/80）*y.^2.*z.^3;
%?isosurface（xyzval0）
%?axis?equal
%?view（-1024）
%?colormap?flag
%?camlight
%?lighting?phong
%?title（‘\fontsize{16}I?Love?You!‘）;
%?axis?off



%?%動畫黑白3D心
%?%構造體積方程和坐標軸，畫出圖形;
%?[XYZ]?=?meshgrid（linspace（-33101））;
%?%3D心型圖方程如下;
%?F?=?-X.^2.*Z.^3-（9/80）.*Y.^2.*Z.^3+（X.^2+（9/4）.*Y.^2+Z.^2-1）.^3;
%?hFigure?=?figure;
%?sz?=?get（hFigure?‘Position‘）;
%?set（hFigure?‘Position‘?[sz（1）-0.15*sz（3）?sz（2）?1.3*sz（3）?sz（4）]）;
%?set（hFigure‘color‘‘w‘?‘menu‘‘none‘）
%?
%?hAxes?=?axes（‘Parent‘hFigure‘NextPlot‘‘add‘‘DataAspectRatio‘[1?1?1]‘XLim‘[30?120]‘YLim‘[35?65]‘ZLim‘[30?75]）;
%?view（[-39?30]）;
%?axis?off
%?%?制作出動態(tài)的隱形效果;
%?hidden?on
%?%?畫出網格，制作網格動態(tài)效果;
%?
%?%?快渲染心得背面:
%?p?=?patch（isosurface（F-0.001））;
%?set（p‘FaceColor‘‘w‘‘EdgeColor‘‘w‘）;
%?
%?%?構造Y-Z平面描完函數在該平面的點:
%?for?iX?=?[35?38?41?45?48?51?54?57?61?64?67]
%?????plane?=?reshape（F（:iX:）101101）;
%?????cData?=?contourc（plane[0?0]）;
%?????xData?=?iX.*ones（1cData（21））;
%?????plot3（hAxesxDatacData（22:end）cData（12:end）‘k‘）;
%?????pause（.1）?drawnow
%?end
%?%?構造X-Z平面描完函數在該平面的點:
%?for?iY?=?[41?44?47?51?55?58?61]
%?????plane?=?reshape（F（iY::）101101）;
%?????cData?=?contourc（plane[0?0]）;
%?????yData?=?iY.*ones（1cData（21））;
%?????plot3（hAxescData（22:end）yDatacData（12:end）‘k‘）;
%?????pause（.1）?drawnow
%?end
%?%?構造X-Y平面描完函數在該平面的點:
%?for?iZ?=?[36?38?40?42?44?46?48?50?52?54?56?58?60?62?64?66?69?71]
%?????plane?=?F（::iZ）;
%?????cData?=?contourc（plane[0?0]）;
%?????startIndex?=?1;
%?????if?size（cData2）?>?（cData（21）+1）
%?????????startIndex?=?cData（21）+2;
%?????????zData?=?iZ.*ones（1cData（21））;
%?????????plot3（hAxescData（12:（startIndex-1））cData（22:（startIndex-1））zData‘k‘）;
%?????end
%?????zData?=?iZ.*ones（1cData（2startIndex））;
%?????plot3（hAxescData（1（startIndex+1）:end）cData（2（startIndex+1）:end）zData‘k‘）;
%?????pause（.1）?drawnow
%?end
%?
%?%函數已經畫完接下來為文字部分