Heya folks, its been almost 5 years since my last submission so I decided to get out and do some fun stuff.
A lot of you will be doing math and calculus once you enter university. One important thing youll learn is coordinate systems!
One such system is the Spherical Polar Coordinates
https://en.wikipedia.org/wiki/Spherical_coordinate_system
Using just two angles we are able to address any point on a R= 1 sphere.. meaning if we were to go through every possible angle--to a degree-- we could figure out the converted points on our cartesian plane!
Woah! Using a bit of logic i decided.. well let me create a sphere first.. but it didnt look well.. it needs an aspect of volume. So for each point i calculate i determine a predetermined width and create a polygon based on our resolution of the sphere.. i even made it.. so the higher the resolution the smaller the blocks looK!
| Turing: |
%%submitted by MihaiG%% var tmax := 5%define the theta subdivission var pmax := 5 %define the phi subdivision var vs := 360 %define how much of the circle around horizontal axis you want to draw var hs := 180 %define how much of the circle around vertical axis you want to draw %ie if you do 180/90 you get a quadrant var pts : array 0 .. (tmax + 1) * (pmax + 1), 1 .. 4, 1 .. 3 of real %defintion of the point storage, first set is the address, second is for defining the four corners, each additional point var a : int := 1 %counter for indexing the array, we could use something like a := theta*tmax + phy var sx := (360 / pmax) / 2 var sy := (180 / tmax) / 2 View.Set ("offscreenonly") setscreen ("graphics:max;max,nobuttonbar") for theta : 0 .. tmax for phi : 0 .. pmax % pts (a, 1, 1) := sind (hs * (theta) / tmax - sy) * cosd (vs * (phi) / pmax - sx) pts (a, 1, 2) := sind (hs * (theta) / tmax - sy) * sind (vs * (phi) / pmax - sx) pts (a, 1, 3) := cosd (hs * (theta) / tmax - sy) % pts (a, 3, 1) := sind (hs * (theta) / tmax + sy) * cosd (vs * (phi) / pmax + sx) pts (a, 3, 2) := sind (hs * (theta) / tmax + sy) * sind (vs * (phi) / pmax + sx) pts (a, 3, 3) := cosd (hs * (theta) / tmax + sy) % pts (a, 2, 1) := sind (hs * (theta) / tmax + sy) * cosd (vs * (phi) / pmax - sx) pts (a, 2, 2) := sind (hs * (theta) / tmax + sy) * sind (vs * (phi) / pmax - sx) pts (a, 2, 3) := cosd (hs * (theta) / tmax + sy) % pts (a, 4, 1) := sind (hs * (theta) / tmax - sy) * cosd (vs * (phi) / pmax + sx) pts (a, 4, 2) := sind (hs * (theta) / tmax - sy) * sind (vs * (phi) / pmax + sx) pts (a, 4, 3) := cosd (hs * (theta) / tmax - sy) % a += 1 end for end for var x, y : array 1 .. 4 of int loop cls put tmax put pmax for i : 1 .. a - 1 for p : 1 .. 4 pts (i, p, 1) := pts (i, p, 1) * 0.999847695 + pts (i, p, 2) * 0.0174524064 pts (i, p, 2) := pts (i, p, 2) * 0.999847695 - ((pts (i, p, 1) - pts (i, p, 2) * 0.0174524064) / 0.999847695) * 0.0174524064 pts (i, p, 3) := pts (i, p, 3) * 0.999847695 + pts (i, p, 2) * 0.0174524064 pts (i, p, 2) := pts (i, p, 2) * 0.999847695 - ((pts (i, p, 3) - pts (i, p, 2) * 0.0174524064) / 0.999847695) * 0.0174524064 x (p) := round (5 * pts (i, p, 1) * 10 / (pts (i, p, 3) * 0.00001 + 0.2) + maxx / 2) y (p) := round (5 * pts (i, p, 2) * 10 / (pts (i, p, 3) * 0.00001 + .2) + maxy / 2) Draw.Dot (x (p), y (p), 7) end for Draw.Polygon (x, y, 4, 7) end for View.Update end loop |
The rotation stuff ill leave as an exercise to the reader.. If you have any questions or improvements let me know! This was about 2-4 hours of work through about 3 different revisions!
[/img]