A股上市公司传智教育(股票代码 003032)旗下技术交流社区北京昌平校区

 找回密码
 加入黑马

QQ登录

只需一步,快速开始

© xu不是许 中级黑马   /  2015-10-27 23:29  /  1500 人查看  /  7 人回复  /   0 人收藏 转载请遵从CC协议 禁止商业使用本文




  1. #include <stdio.h>
  2. #include <math.h>

  3. float f(float x, float y, float z) {
  4.     float a = x * x + 9.0f / 4.0f * y * y + z * z - 1;
  5.     return a * a * a - x * x * z * z * z - 9.0f / 80.0f * y * y * z * z * z;
  6. }

  7. float h(float x, float z) {
  8.     for (float y = 1.0f; y >= 0.0f; y -= 0.001f)
  9.         if (f(x, y, z) <= 0.0f)
  10.             return y;
  11.     return 0.0f;
  12. }

  13. int main() {
  14.     for (float z = 1.5f; z > -1.5f; z -= 0.05f) {
  15.         for (float x = -1.5f; x < 1.5f; x += 0.025f) {
  16.             float v = f(x, 0.0f, z);
  17.             if (v <= 0.0f) {
  18.                 float y0 = h(x, z);
  19.                 float ny = 0.01f;
  20.                 float nx = h(x + ny, z) - y0;
  21.                 float nz = h(x, z + ny) - y0;
  22.                 float nd = 1.0f / sqrtf(nx * nx + ny * ny + nz * nz);
  23.                 float d = (nx + ny - nz) * nd * 0.5f + 0.5f;
  24.                 putchar(".:-=+*#%@"[(int)(d * 5.0f)]);
  25.             }
  26.             else
  27.                 putchar(' ');
  28.         }
  29.         putchar('\n');
  30.     }
  31. }
复制代码

7 个回复

倒序浏览
不赖,好高大上
回复 使用道具 举报
好牛叉呀!厉害!学了!
回复 使用道具 举报
给点注解 会更完美
回复 使用道具 举报
真当没人上知乎?
回复 使用道具 举报
厉害!!!
回复 使用道具 举报
艺术感很强啊
回复 使用道具 举报 1 0
6666666666666666666
回复 使用道具 举报
您需要登录后才可以回帖 登录 | 加入黑马