using System;
using System.Collections.Generic;
using System.ComponentModel;
using System.Data;
using System.Drawing;
using System.Linq;
using System.Text;
using System.Threading.Tasks;
using System.Windows.Forms;
namespace WindowsFormsApplication1
{
public partial class Form1 : Form
{
public Form1()
{
InitializeComponent();
}
public double[,] mult(double[,] a, double[,] b)
{
double[,] c = new double[8, 4];
for (int i = 0; i < 8; ++i)
{
for (int j = 0; j < 4; ++j)
{
for (int k = 0; k < 4; ++k)
{
c[i, j] += a[i, k] * b[k, j];
}
}
}
return c;
}
double[,] cube = { {0.0, 0.0, 0.0, 1.0},
{100.0, 0.0, 0.0, 1.0},
{100.0, 0.0, 100.0, 1.0},
{0.0, 0.0, 100.0, 1.0},
{0.0, 100.0, 0.0, 1.0},
{100.0, 100.0, 0.0, 1.0},
{100.0, 100.0, 100.0, 1.0},
{0.0, 100.0, 100.0, 1.0} };
double[] pt = { 1.0, 1.0, 1.0 };
double n1, n2, n3;
double alpha = 0.0;
private void Form1_Paint(object sender, PaintEventArgs e)
{
double[,] rotate = new double[4, 4];
n1 = pt[0] / Math.Sqrt(pt[0] * pt[0] + pt[1] * pt[1] + pt[2] * pt[2]);
n2 = pt[1] / Math.Sqrt(pt[0] * pt[0] + pt[1] * pt[1] + pt[2] * pt[2]);
n3 = pt[2] / Math.Sqrt(pt[0] * pt[0] + pt[1] * pt[1] + pt[2] * pt[2]);
rotate[0, 0] = n1 * n1 + (1.0 - n1 * n1) * Math.Cos(alpha);
rotate[1, 0] = n1 * n2 * (1.0 - Math.Cos(alpha)) - n3 * Math.Sin(alpha);
rotate[2, 0] = n1 * n3 * (1.0 - Math.Cos(alpha)) + n2 * Math.Sin(alpha);
rotate[3, 0] = 0.0;
rotate[0, 1] = n1 * n2 * (1.0 - Math.Cos(alpha)) + n3 * Math.Sin(alpha);
rotate[1, 1] = n2 * n2 + (1.0 - n2 * n2) * Math.Cos(alpha);
rotate[2, 1] = n2 * n3 * (1.0 - Math.Cos(alpha)) - n1 * Math.Sin(alpha);
rotate[3, 1] = 0.0;
rotate[0, 2] = n1 * n3 * (1.0 - Math.Cos(alpha)) - n2 * Math.Sin(alpha);
rotate[1, 2] = n2 * n3 * (1.0 - Math.Cos(alpha)) + n1 * Math.Sin(alpha);
rotate[2, 2] = n3 * n3 + (1.0 - n3 * n3) * Math.Cos(alpha);
rotate[3, 2] = 0.0;
rotate[0, 3] = 0.0;
rotate[1, 3] = 0.0;
rotate[2, 3] = 0.0;
rotate[3, 3] = 1.0;
double[,] output = mult(cube, rotate);
Graphics g = e.Graphics;
g.TranslateTransform(this.ClientRectangle.Width / 2, this.ClientRectangle.Height / 2);
g.DrawLine(new Pen(Color.Black), new PointF((float)output[0, 0], (float)output[0, 1]), new PointF((float)output[1, 0], (float)output[1, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[1, 0], (float)output[1, 1]), new PointF((float)output[2, 0], (float)output[2, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[2, 0], (float)output[2, 1]), new PointF((float)output[3, 0], (float)output[3, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[3, 0], (float)output[3, 1]), new PointF((float)output[0, 0], (float)output[0, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[4, 0], (float)output[4, 1]), new PointF((float)output[5, 0], (float)output[5, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[5, 0], (float)output[5, 1]), new PointF((float)output[6, 0], (float)output[6, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[6, 0], (float)output[6, 1]), new PointF((float)output[7, 0], (float)output[7, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[7, 0], (float)output[7, 1]), new PointF((float)output[4, 0], (float)output[4, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[0, 0], (float)output[0, 1]), new PointF((float)output[4, 0], (float)output[4, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[1, 0], (float)output[1, 1]), new PointF((float)output[5, 0], (float)output[5, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[2, 0], (float)output[2, 1]), new PointF((float)output[6, 0], (float)output[6, 1]));
g.DrawLine(new Pen(Color.Black), new PointF((float)output[3, 0], (float)output[3, 1]), new PointF((float)output[7, 0], (float)output[7, 1]));
}
private void timer1_Tick(object sender, EventArgs e)
{
this.Invalidate();
this.DoubleBuffered = true;
alpha += 0.01;
}
private void Form1_Load(object sender, EventArgs e)
{
}
}
}