闲来无事，用matlab来演绎一下数学的浪漫 \(≧▽≦)/ 提供matlab代码及图形，最后还有一个程序，大家根据图形来猜一下这像什么 ↖(^ω^)↗ （看m文件的命名，或许你就恍然大悟呦~）
先上图
u=ezplot('17*x.^2-16*abs(x)*y+17*y.^2=225');
axis([-5 20 -13 14])
set(u,'color','r')
未有严谨证据证明心脏线是由笛卡尔发明。
心脏线亦为蚶线的一种。在曼德博集合正中间的图形便是一个心脏线。心脏线的英文名称“Cardioid”是 de Castillon 在1741年的《Philosophical Transactions of the Royal Society》发表的；意为“像心脏的”。
数学表达
极坐标方程
水平方向： r=a(1-cosθ)或 r=a(1+cosθ) (a>0)
垂直方向： r=a(1-sinθ)或 r=a(1+sinθ) (a>0)
直角坐标方程
心形线的平面直角坐标系方程表达式分别为 x^2+y^2+a*x=a*sqrt(x^2+y^2) 和x^2+y^2-a*x=a*sqrt(x^2+y^2）
参数方程
x=a*(2*cos(t)-cos(2*t))
y=a*(2*sin(t)-sin(2*t))
所围面积为3/2*PI*a^2，形成的弧长为8a
故事
《数学的故事》里面说到了数学家笛卡尔的爱情故事。笛卡尔于1596年出生在法国，欧洲大陆爆发黑死病时他流浪到瑞典，
1649年，斯德哥尔摩的街头，52岁的笛卡尔邂逅了18岁的瑞典公主克里斯汀。几天后，他意外的接到通知，国王聘请他做小公主的数学老师。跟随前来通知的侍卫一起来到皇宫，他见到了在街头偶遇的女孩子。从此，他当上了小公主的数学老师。
小公主的数学在笛卡尔的悉心指导下突飞猛进，笛卡尔向她介绍了自己 研究的新领域--直角坐标系。每天形影不离的相处使他们彼此产生爱慕之心，公主的父亲国王知道了后勃然大怒，下令将笛卡尔处死，小公主克里斯汀苦苦哀求后，国王将其流放回法国，克里斯汀公主也被父亲软禁起来。
笛卡尔回法国后不久便染上重病，他日日给公主写信，因被国王拦截，克里斯汀一直没收到笛卡尔的信。笛卡尔在给克里斯汀寄出第十三封信后就气绝身亡了，这第十三封信内容只有短短的一个公式：r=a(1-sinθ）。国王看不懂，觉得他们俩之间并不是总是说情话的，将全城的数学家召集到皇宫，但没有一个人能解开，他不忍心看着心爱的女儿整日闷闷不乐，就把这封信交给一直闷闷不乐的克里斯汀。
公主看到后，立即明了恋人的意图，她马上着手把方程的图形画出来，看到图形，她开心极了，她知道恋人仍然爱着她，原来方程的图形是一颗心的形状。这也就是著名的“心形线”。
国王死后，克里斯汀登基，立即派人在欧洲四处寻找心上人，无奈斯人已故，先她一步走了，徒留她孤零零在人间...
据说这封享誉世界的另类情书还保存在欧洲笛卡尔的纪念馆里。
再来一个三维立体的
[x,y,z]=meshgrid(linspace(-5,5,120));
f=(x.^2+(9*y.^2)./4+z.^2-1).^3-((9*y.^2).*(z.^3))./80-(x.^2).*(z.^3);
p=patch(isosurface(x,y,z,f,0));
set(p,'FaceColor','r')
grid on
daspect([1 1 1])
view(-30,24)
哈哈，大家开动脑洞，来看看这是什么
%Rosa_rugosa.m
a=1;
x=(-pi:pi/100:pi);
y=a*(1-sin(x));
polar(y,'r')
