- DBPro Coding Challenges - - GameCreators Forum

Author

Message

qwe

22

Years of Service

User Offline

Joined: 3rd Sep 2003

Location: place

Posted: 9th Jul 2007 21:53

Link

the way my cars turn is choppy, theres only one kind of car, and theres an occasional glitch at the stoplights. unfortunately, we are about to leave for a family vacation, and i'm not sure if we'll have internet. if we do i'll update my code before midnight tonight

the reason it gets messy if you put in too many more cars is that once the car gets within a certain distance from the stoplight (25 units IIRC), it gets a place in line at the stoplight, and goes to that place. if there are too many cars, they come within 40 units, then realize their place in line is maybe 45 or 50 units away from the stoplight because there are so many cars there already, and back up to that place. so if you put in too many cars, you have to change that number (its on lines 905, 928, 951, and 974)

Quote: "the thing I thought was a bit funny in qwe's entry though, is that I swore I saw two cars got caught in the middle of a junction when the lights turned red, and they both reversed back to where they came from to clear the junction!"

hehe.. my cars are smart! any unfixed glitches happen, i think, because the cars are placed randomly along the road, so sometimes two cars happen to be placed inside of eachother or in the middle of an intersection at start, and then weird things happen

if i win and i'm not here (i'll be back thursday or friday), or if someone else wins and wants a suggestion for the next challenge,

let's brainstorm science ideas. physical simulations or Spelunx-like rooms that would be useful in the classroom and fun to play around with (learning in disguise)?

Back to top

Profile PM Email

Ric

22

Years of Service

User Offline

Joined: 11th Jul 2004

Location: object position x

Posted: 9th Jul 2007 22:30

Link

Quote: "it'd seem to me that people would still be able to post entries unti midnight today, and i'll judge tomorrow?
is this correct!?"

Yep - that's how it usually goes. If you leave it until tomorrow, people in different time zones may get something last minute done before they go to bed.

Back to top

Profile PM Email Website

Pricey

23

Years of Service

User Offline

Joined: 22nd Feb 2003

Location:

Posted: 9th Jul 2007 23:03

Link

okay sounds good

errr in your entry gadgetfreak, i couldn't work out how to place traffic lights?! i click the option but when i try to place them, nothing!

perhaps you would be able to clarify how to place them
or fix the bug before midnight?

because it is just one minor flaw in an excellent entry!

good luck everyone!

http://forum.thegamecreators.com/?m=forum_view&t=97191&b=32
God Of Thunder DEMO

Back to top

Profile PM Email

gadgetfreak 66

19

Years of Service

User Offline

Joined: 26th Sep 2006

Location: Kulin,WA, Australia

Posted: 10th Jul 2007 05:11

Link

They need to be on an 4-way intersection of road, traffic lights generally don't work otherwise.

Back to top

Profile PM Email

Spopovich

20

Years of Service

User Offline

Joined: 12th Jul 2006

Location:

Posted: 10th Jul 2007 06:51 Edited at: 10th Jul 2007 07:21

Link

First of all, all of the entries look great, so good luck to everyone. I fixed a few bugs and figured I should post my finished product. This has been a great challenge.

sync on:sync rate 0:sync

set display mode 1024,768,16

randomize timer()

make camera 1
color backdrop 1,0x00bfff
set camera range 1,1,500
position light 0,0,600,0
make light 1
position light 1,0,-600,0
make light 2
position light 2,0,0,-600
make light 3
position light 3,0,0,600
make light 4
position light 4,-600,0,0
make light 5
position light 5,600,0,0

move camera 1,-1000
pitch camera up 1,90
move camera 1,10
pitch camera down 1,90

global stopNum
global col_check
global car
global dim cars(1000,7)
global dim lights(100,2)

ink 0,0

make_first_car()

make_textures()

make_stoplight()

build_city2()

make_cars(300)

do

fastsync
   text 0,0,"FPS "+str$(screen fps())

if mouseclick()=1
      if clicked=0
         pitch camera up 1,90
         move camera 1,150
         pitch camera down 1,180
         clicked=1
      endif
   else
      if clicked=1
         clicked=0
         move camera 1,150
         pitch camera up 1,90
      endif
      control camera using arrowkeys 1,8,3
   endif

control_cars()

work_stoplight()

loop

function make_textures()

ink rgb(178,34,34),rgb(178,34,34)
   box 1,1,10,5
   ink rgb(100,100,100),rgb(100,100,100)
   line 1,1,10,1
   line 1,1,1,5
   line 1,5,10,5
   line 10,1,10,5
   get image 1,0,0,11,6,1
   cls

for j=0 to 150
      for i=0 to 51
         if j mod 2
            paste image 1,i*10,j*5
         else
            paste image 1,i*10-5,j*5
         endif
      next i
   next j

ink rgb(176,226,255),rgb(176,226,255)
   windows()

get image 1,0,0,500,750,1

cls

ink rgb(255,231,186),rgb(255,231,186)
   box 0,0,500,750
   ink rgb(54,100,139),rgb(54,100,139)
   windows()

get image 2,0,0,500,750,1

cls

ink rgb(20,20,20),rgb(20,20,20)
   box 0,0,500,750

ink rgb(176,226,255),rgb(176,226,255)
   windows()

get image 3,0,0,500,750,1
   cls
   sync

endfunction

function windows()

for j=5 to 145 step 10
      for i=1 to 50 step 5
         box i*10,j*5,i*10+5,j*5+10
      next i
   next j

endfunction

function control_cars()

for i=1 to car
      cars(i,5)=object position x(i+5000)
      cars(i,6)=object position z(i+5000)
   next i

col_check=col_check+1
   if col_check>5
      col_check=1
   endif

for i=1 to car

a=wrapvalue(object angle y(i+5000))

if cars(i,3)<>2
            if a>85 and z<95
               cars(i,1)=1
               yrotate object i+5000,90
            endif
            if a>265 and a<275
               cars(i,1)=1
               yrotate object i+5000,270
            endif
            if a>175 and a<185
               cars(i,1)=0
               yrotate object i+5000,180
            endif
            if (a>355 and a<365) or (a>-5 and a<5)
               cars(i,1)=0
               yrotate object i+5000,0
            endif
         endif

moved=0
         no=0
         x=cars(i,5)
         y=cars(i,6)

if y>650 or y<-650 or x>650 or x<-650
            move object i+5000,-1250
         endif

for j=1 to stopNum
            x2=object position x(j+1000)
            y2=object position z(j+1000)
            if x>x2-25 and x<x2+25 and y>y2-25 and y<y2+25
               no=1

if cars(i,3)=0
                  cars(i,3)=rnd(1)+1
                  if cars(i,3)=2
                     cars(i,4)=24
                  endif
               endif

if lights(j,2)=cars(i,1)
                  move object i+5000,1

if cars(i,3)=2
                     if cars(i,4)>0
                        turn object right i+5000,3.75`yrotate object i+5000,object angle y(i+5000)+10
                        cars(i,4)=cars(i,4)-1
                     else
                        if cars(i,1)=1
                           cars(i,1)=0
                        else
                           cars(i,1)=1
                        endif
                        cars(i,3)=1
                        cars(i,4)=100
                     endif
                  endif

moved=1
                  cars(i,2)=1
               else
                  if cars(i,2)=1
                     move object i+5000,1
                     moved=1
                     if cars(i,3)=2
                        if cars(i,4)>0
                           turn object right i+5000,3.75`yrotate object i+5000,object angle y(i+5000)+10
                           cars(i,4)=cars(i,4)-1
                        else
                           if cars(i,1)=1
                              cars(i,1)=0
                           else
                              cars(i,1)=1
                           endif
                           cars(i,3)=1
                           cars(i,4)=100
                        endif
                     endif
                  endif

endif
            endif
         next j

if no=0
            move object i+5000,1
            moved=1
            if cars(i,3)=2
               move object right i+5000,1.78
            endif
            cars(i,2)=0
            cars(i,3)=0
            cars(i,4)=0
         endif

if moved

for j=1 to car
                if i<>j
                   if car_col(i,j)

if cars(i,3)=2
                        cars(i,4)=cars(i,4)+1
                        turn object right i+5000,-3.75
                     endif

move object i+5000,-1

endif
                endif
            next j
         endif

next i

endfunction

function car_col(num1,num2)

dist=(cars(num2,5)-cars(num2,5))^2+(cars(num1,6)-cars(num2,6))^2

if dist<400
      yes=object collision(num1+5000,num2+5000)
   endif

endfunction yes

function make_cars(n)

for i=1 to n
      make_car()
   next i

endfunction

function make_first_car()

make object box 4999,5,5,10
   hide object 4999
   make object collision box 4999,-2.5,0,-7.5,2.5,5,7.5,1

endfunction

function make_car()

car=car+1

instance object car+5000,4999
   show object car+5000

repeat
      redo=0
      a=rnd(3)
      select a
         case 0
            x=rnd(9)-5
            x=x*100
            cars(car,1)=0

y=rnd(1100)-550
            yrotate object car+5000,0
         endcase
         case 1
            x=rnd(9)-5
            x=x*100
            cars(car,1)=0

y=rnd(1100)-550

yrotate object car+5000,180
         endcase
         case 2
            y=rnd(9)-5
            y=y*100
            cars(car,1)=1
            x=rnd(1100)-550

yrotate object car+5000,90
         endcase
         case 3
            y=rnd(9)-5
            y=y*100
            cars(car,1)=1

x=rnd(1100)-550

yrotate object car+5000,270
         endcase
      endselect

position object car+5000,x,2.6,y
      move object right car+5000,10

for j=1 to stopNum
         x2=object position x(j+1000)
         y2=object position z(j+1000)
         if x>x2-25 and x<x2+25 and y>y2-25 and y<y2+25
            redo=1
         endif
      next j

for j=1 to car-2
          if object collision(car+5000,j+5000)
            redo=1
          endif
      next j

until redo=0

cars(car,2)=1

endfunction

function build_city2()

make object plain 1,1200,1100
   pitch object down 1,90
   position object 1,0,0,-50
   color object 1,0

num=1

for x=0 to 10
      for y=0 to 10
         num=num+1
         make object box num,60,rnd(200-50)+50,60
         position object num,x*100-550,object size y(num)/2,y*100-550
         texture object num,rnd(2)+1
      next y
   next x

for x=-5 to 4
      for y=-5 to 4
         make_stoplights(x*100,y*100)
      next y
   next x

endfunction

function make_stoplights(x,y)

stopNum=stopNum+1
   objnum=stopNum+1000

lights(stopNum,1)=rnd(300)
   lights(stopNum,2)=rnd(1)
   clone object objnum,1000
   position object objnum,x,object size y(objnum)/5,y

endfunction

function make_stoplight()

objnum=1000

make object cylinder objnum,10
   scale limb objnum,0,5,100,5
   make object box objnum+1,2,3,2
   make mesh from object 1,objnum+1
   delete object objnum+1
   add limb objnum,1,1
   position object objnum,0,5,0
   offset limb objnum,1,0,5,0
   delete mesh 1
   make object sphere objnum+1,1
   make mesh from object 2,objnum+1
   delete object objnum+1

for i=2 to 13
      add limb objnum,i,2
   next i
   delete mesh 2

for i=1 to 3
      offset limb objnum,i+1,0,i+3,-1
   next i
   for i=1 to 3
      offset limb objnum,i+4,0,i+3,1
   next i
   for i=1 to 3
      offset limb objnum,i+7,1,i+3,0
   next i
   for i=1 to 3
      offset limb objnum,i+10,-1,i+3,0
   next i

for i=2 to 13 step 3
      color limb objnum,i,rgb(0,100,0)
   next i

for i=3 to 13 step 3
      color limb objnum,i,rgb(100,100,0)
   next i

for i=4 to 13 step 3
      color limb objnum,i,rgb(100,0,0)
   next i

hide object objnum

endfunction

function work_stoplight()

for j=1 to stopNum

change=0

objnum=j+1000

lights(j,1)=lights(j,1)+1

if lights(j,1)=275
         if lights(j,2)=1
            lights(j,2)=4
         else
            lights(j,2)=3
         endif
      endif

if lights(j,1)>300
         lights(j,1)=0
         if lights(j,2)=4
            lights(j,2)=0
         else
            lights(j,2)=1
         endif
         change=1
      endif

if change=1 or lights(j,1)=250
         for i=2 to 13 step 3
            color limb objnum,i,rgb(0,100,0)
         next i

for i=3 to 13 step 3
            color limb objnum,i,rgb(100,100,0)
         next i

for i=4 to 13 step 3
            color limb objnum,i,rgb(100,0,0)
         next i

select lights(j,2)
            case 0

if lights(j,1)<250
                  color limb objnum,2,rgb(0,255,0)
                  color limb objnum,5,rgb(0,255,0)
               else
                  color limb objnum,3,rgb(255,255,0)
                  color limb objnum,6,rgb(255,255,0)
               endif

color limb objnum,10,rgb(255,0,0)
               color limb objnum,13,rgb(255,0,0)
            endcase
            case 1
               if lights(j,1)<250
                  color limb objnum,8,rgb(0,255,0)
                  color limb objnum,11,rgb(0,255,0)
               else
                  color limb objnum,9,rgb(255,255,0)
                  color limb objnum,12,rgb(255,255,0)
               endif
               color limb objnum,4,rgb(255,0,0)
               color limb objnum,7,rgb(255,0,0)
            endcase
         endselect
      endif

next j

endfunction

+ Code Snippet


sync on:sync rate 0:sync

set display mode 1024,768,16

randomize timer()



make camera 1
color backdrop 1,0x00bfff
set camera range 1,1,500
position light 0,0,600,0
make light 1
position light 1,0,-600,0
make light 2
position light 2,0,0,-600
make light 3
position light 3,0,0,600
make light 4
position light 4,-600,0,0
make light 5
position light 5,600,0,0



move camera 1,-1000
pitch camera up 1,90
move camera 1,10
pitch camera down 1,90




global stopNum
global col_check
global car
global dim cars(1000,7)
global dim lights(100,2)

ink 0,0


make_first_car()

make_textures()

make_stoplight()

build_city2()

make_cars(300)




do

   fastsync
   text 0,0,"FPS "+str$(screen fps())

   if mouseclick()=1
      if clicked=0
         pitch camera up 1,90
         move camera 1,150
         pitch camera down 1,180
         clicked=1
      endif
   else
      if clicked=1
         clicked=0
         move camera 1,150
         pitch camera up 1,90
      endif
      control camera using arrowkeys 1,8,3
   endif

   control_cars()

   work_stoplight()



loop



function make_textures()

   ink rgb(178,34,34),rgb(178,34,34)
   box 1,1,10,5
   ink rgb(100,100,100),rgb(100,100,100)
   line 1,1,10,1
   line 1,1,1,5
   line 1,5,10,5
   line 10,1,10,5
   get image 1,0,0,11,6,1
   cls



   for j=0 to 150
      for i=0 to 51
         if j mod 2
            paste image 1,i*10,j*5
         else
            paste image 1,i*10-5,j*5
         endif
      next i
   next j

   ink rgb(176,226,255),rgb(176,226,255)
   windows()

   get image 1,0,0,500,750,1

   cls

   ink rgb(255,231,186),rgb(255,231,186)
   box 0,0,500,750
   ink rgb(54,100,139),rgb(54,100,139)
   windows()

   get image 2,0,0,500,750,1

   cls


   ink rgb(20,20,20),rgb(20,20,20)
   box 0,0,500,750

   ink rgb(176,226,255),rgb(176,226,255)
   windows()

   get image 3,0,0,500,750,1
   cls
   sync

endfunction

function windows()


   for j=5 to 145 step 10
      for i=1 to 50 step 5
         box i*10,j*5,i*10+5,j*5+10
      next i
   next j

endfunction

function control_cars()

   for i=1 to car
      cars(i,5)=object position x(i+5000)
      cars(i,6)=object position z(i+5000)
   next i


   col_check=col_check+1
   if col_check>5
      col_check=1
   endif


   for i=1 to car


         a=wrapvalue(object angle y(i+5000))


         if cars(i,3)<>2
            if a>85 and z<95
               cars(i,1)=1
               yrotate object i+5000,90
            endif
            if a>265 and a<275
               cars(i,1)=1
               yrotate object i+5000,270
            endif
            if a>175 and a<185
               cars(i,1)=0
               yrotate object i+5000,180
            endif
            if (a>355 and a<365) or (a>-5 and a<5)
               cars(i,1)=0
               yrotate object i+5000,0
            endif
         endif



         moved=0
         no=0
         x=cars(i,5)
         y=cars(i,6)

         if y>650 or y<-650 or x>650 or x<-650
            move object i+5000,-1250
         endif


         for j=1 to stopNum
            x2=object position x(j+1000)
            y2=object position z(j+1000)
            if x>x2-25 and x<x2+25 and y>y2-25 and y<y2+25
               no=1

               if cars(i,3)=0
                  cars(i,3)=rnd(1)+1
                  if cars(i,3)=2
                     cars(i,4)=24
                  endif
               endif

               if lights(j,2)=cars(i,1)
                  move object i+5000,1

                  if cars(i,3)=2
                     if cars(i,4)>0
                        turn object right i+5000,3.75`yrotate object i+5000,object angle y(i+5000)+10
                        cars(i,4)=cars(i,4)-1
                     else
                        if cars(i,1)=1
                           cars(i,1)=0
                        else
                           cars(i,1)=1
                        endif
                        cars(i,3)=1
                        cars(i,4)=100
                     endif
                  endif


                  moved=1
                  cars(i,2)=1
               else
                  if cars(i,2)=1
                     move object i+5000,1
                     moved=1
                     if cars(i,3)=2
                        if cars(i,4)>0
                           turn object right i+5000,3.75`yrotate object i+5000,object angle y(i+5000)+10
                           cars(i,4)=cars(i,4)-1
                        else
                           if cars(i,1)=1
                              cars(i,1)=0
                           else
                              cars(i,1)=1
                           endif
                           cars(i,3)=1
                           cars(i,4)=100
                        endif
                     endif
                  endif

               endif
            endif
         next j

         if no=0
            move object i+5000,1
            moved=1
            if cars(i,3)=2
               move object right i+5000,1.78
            endif
            cars(i,2)=0
            cars(i,3)=0
            cars(i,4)=0
         endif



         if moved

            for j=1 to car
                if i<>j
                   if car_col(i,j)

                     if cars(i,3)=2
                        cars(i,4)=cars(i,4)+1
                        turn object right i+5000,-3.75
                     endif

                     move object i+5000,-1

                   endif
                endif
            next j
         endif


   next i



endfunction

function car_col(num1,num2)

   dist=(cars(num2,5)-cars(num2,5))^2+(cars(num1,6)-cars(num2,6))^2

   if dist<400
      yes=object collision(num1+5000,num2+5000)
   endif

endfunction yes


function make_cars(n)

   for i=1 to n
      make_car()
   next i

endfunction

function make_first_car()

   make object box 4999,5,5,10
   hide object 4999
   make object collision box 4999,-2.5,0,-7.5,2.5,5,7.5,1

endfunction

function make_car()

   car=car+1

   instance object car+5000,4999
   show object car+5000

   repeat
      redo=0
      a=rnd(3)
      select a
         case 0
            x=rnd(9)-5
            x=x*100
            cars(car,1)=0

            y=rnd(1100)-550
            yrotate object car+5000,0
         endcase
         case 1
            x=rnd(9)-5
            x=x*100
            cars(car,1)=0

            y=rnd(1100)-550

            yrotate object car+5000,180
         endcase
         case 2
            y=rnd(9)-5
            y=y*100
            cars(car,1)=1
            x=rnd(1100)-550

            yrotate object car+5000,90
         endcase
         case 3
            y=rnd(9)-5
            y=y*100
            cars(car,1)=1

            x=rnd(1100)-550

            yrotate object car+5000,270
         endcase
      endselect

      position object car+5000,x,2.6,y
      move object right car+5000,10

      for j=1 to stopNum
         x2=object position x(j+1000)
         y2=object position z(j+1000)
         if x>x2-25 and x<x2+25 and y>y2-25 and y<y2+25
            redo=1
         endif
      next j

      for j=1 to car-2
          if object collision(car+5000,j+5000)
            redo=1
          endif
      next j

   until redo=0

   cars(car,2)=1


endfunction


function build_city2()

   make object plain 1,1200,1100
   pitch object down 1,90
   position object 1,0,0,-50
   color object 1,0

   num=1

   for x=0 to 10
      for y=0 to 10
         num=num+1
         make object box num,60,rnd(200-50)+50,60
         position object num,x*100-550,object size y(num)/2,y*100-550
         texture object num,rnd(2)+1
      next y
   next x

   for x=-5 to 4
      for y=-5 to 4
         make_stoplights(x*100,y*100)
      next y
   next x

endfunction

function make_stoplights(x,y)

   stopNum=stopNum+1
   objnum=stopNum+1000


   lights(stopNum,1)=rnd(300)
   lights(stopNum,2)=rnd(1)
   clone object objnum,1000
   position object objnum,x,object size y(objnum)/5,y

endfunction

function make_stoplight()

   objnum=1000



   make object cylinder objnum,10
   scale limb objnum,0,5,100,5
   make object box objnum+1,2,3,2
   make mesh from object 1,objnum+1
   delete object objnum+1
   add limb objnum,1,1
   position object objnum,0,5,0
   offset limb objnum,1,0,5,0
   delete mesh 1
   make object sphere objnum+1,1
   make mesh from object 2,objnum+1
   delete object objnum+1



   for i=2 to 13
      add limb objnum,i,2
   next i
   delete mesh 2

   for i=1 to 3
      offset limb objnum,i+1,0,i+3,-1
   next i
   for i=1 to 3
      offset limb objnum,i+4,0,i+3,1
   next i
   for i=1 to 3
      offset limb objnum,i+7,1,i+3,0
   next i
   for i=1 to 3
      offset limb objnum,i+10,-1,i+3,0
   next i

   for i=2 to 13 step 3
      color limb objnum,i,rgb(0,100,0)
   next i

   for i=3 to 13 step 3
      color limb objnum,i,rgb(100,100,0)
   next i

   for i=4 to 13 step 3
      color limb objnum,i,rgb(100,0,0)
   next i

   hide object objnum



endfunction

function work_stoplight()

   for j=1 to stopNum

      change=0

      objnum=j+1000

      lights(j,1)=lights(j,1)+1

      if lights(j,1)=275
         if lights(j,2)=1
            lights(j,2)=4
         else
            lights(j,2)=3
         endif
      endif

      if lights(j,1)>300
         lights(j,1)=0
         if lights(j,2)=4
            lights(j,2)=0
         else
            lights(j,2)=1
         endif
         change=1
      endif

      if change=1 or lights(j,1)=250
         for i=2 to 13 step 3
            color limb objnum,i,rgb(0,100,0)
         next i

         for i=3 to 13 step 3
            color limb objnum,i,rgb(100,100,0)
         next i

         for i=4 to 13 step 3
            color limb objnum,i,rgb(100,0,0)
         next i

         select lights(j,2)
            case 0

               if lights(j,1)<250
                  color limb objnum,2,rgb(0,255,0)
                  color limb objnum,5,rgb(0,255,0)
               else
                  color limb objnum,3,rgb(255,255,0)
                  color limb objnum,6,rgb(255,255,0)
               endif

               color limb objnum,10,rgb(255,0,0)
               color limb objnum,13,rgb(255,0,0)
            endcase
            case 1
               if lights(j,1)<250
                  color limb objnum,8,rgb(0,255,0)
                  color limb objnum,11,rgb(0,255,0)
               else
                  color limb objnum,9,rgb(255,255,0)
                  color limb objnum,12,rgb(255,255,0)
               endif
               color limb objnum,4,rgb(255,0,0)
               color limb objnum,7,rgb(255,0,0)
            endcase
         endselect
      endif

   next j

endfunction

68% of rupublicans don't beleive in evolution but, only 5% of monkeys believe in republicans. ~Stephen Colbert 2007

Back to top

Profile PM

qwe

22

Years of Service

User Offline

Joined: 3rd Sep 2003

Location: place

Posted: 11th Jul 2007 02:47

Link

^nice simulation. however, they never turn left:p

Back to top

Profile PM Email

Spopovich

20

Years of Service

User Offline

Joined: 12th Jul 2006

Location:

Posted: 11th Jul 2007 03:26

Link

Yea left turns were gonna be a lot harder to implement because they would have to pass oncoming traffic, which wouldn't have ended well. Oh and when I posted it last night, I had left a bug in it so that it wouldn't load. I was in bed thinking about it and I was like oh shoot. I had to get up and fixed it.

68% of rupublicans don't beleive in evolution but, only 5% of monkeys believe in republicans. ~Stephen Colbert 2007

Back to top

Profile PM

qwe

22

Years of Service

User Offline

Joined: 3rd Sep 2003

Location: place

Posted: 11th Jul 2007 03:38

Link

i have to come all the way to the inn to access the internet. i finally get here and he hasn't picked the winner yet

maybe i can buy points

i

pricey

Back to top

Profile PM Email

Pricey

23

Years of Service

User Offline

Joined: 22nd Feb 2003

Location:

Posted: 11th Jul 2007 11:25 Edited at: 11th Jul 2007 20:05

Link

okay! sorry about the delay! having a few computer problems!
just need to get db pro installed so that i can play the entries through again to make my decision!

EDIT

okayyy people!

i have decided my winner shalst be

+ Code Snippet

gadget freak
this was definitely the coolest entry and really impressed me
i loved the way you could build the roads, the fact that it can work with virtually any map you make, and is not restricted to a grid system is just great! this could easily be turned into a great GTA or SimCity style game! well done!

my runner up is

+ Code Snippet

acelepage
woww! i was impressed by the technicality of this entry!! i doubt it would be too good for a game, but for observing traffic patterns and planning road systems it could become a great tool. i was impressed by the way it could build the roads from the data commands, but maybe this would make it hard to edit?

everyone else well done!
there is the making of some great projects there, but i think the main thing that held them back was the dependency on a simple grid system.

qwe: good entry apart from the thing about reversing back when the light changes mid-junction which I'm sure is not safe, isn't that what the amber light is for? the amber light should come on to allow any cars that are already in the junction to get across before it goes red.

Spopovich: good entry, but a few bugs with cars becoming stuck in things and crashing at junctions! again, maybe a detection for amber lights could be used to make sure this doesn't happen!

thanks everyone for all your great entries!

i hope i judged that fairly...

http://forum.thegamecreators.com/?m=forum_view&t=97191&b=32
God Of Thunder DEMO

Back to top

Profile PM Email

Spopovich

20

Years of Service

User Offline

Joined: 12th Jul 2006

Location:

Posted: 11th Jul 2007 20:21

Link

Congrats gadget freak your entry was amazing. A fun challenge, can't wait for the next one.

68% of rupublicans don't beleive in evolution but, only 5% of monkeys believe in republicans. ~Stephen Colbert 2007

Back to top

Profile PM

gadgetfreak 66

19

Years of Service

User Offline

Joined: 26th Sep 2006

Location: Kulin,WA, Australia

Posted: 12th Jul 2007 08:53 Edited at: 12th Jul 2007 08:58

Link

Alright! Thanks Pricey!

Now, ideas anyone? Cause as it has been previously stated, my imaginative ability is basically non-existant, so if anyone has had an interesting idea lately, speak up now!

Back to top

Profile PM Email

gadgetfreak 66

19

Years of Service

User Offline

Joined: 26th Sep 2006

Location: Kulin,WA, Australia

Posted: 12th Jul 2007 08:53 Edited at: 12th Jul 2007 08:56

Link

(edited)
Sigh. Dial-up. The bane of my existence. Sorry everyone, double post.

Back to top

Profile PM Email

Cave Man

19

Years of Service

User Offline

Joined: 22nd Aug 2006

Location: North Carolina, US

Posted: 12th Jul 2007 16:41 Edited at: 12th Jul 2007 17:23

Link

Do a ecosystem challenge that involves AI and stuff. That would be fun.

Back to top

Profile PM

LongFist

22

Years of Service

User Offline

Joined: 25th May 2004

Location: Winston-Salem, NC

Posted: 12th Jul 2007 18:10 Edited at: 12th Jul 2007 18:10

Link

Try as I might, I can never seem to get an entry in on time! But I graduate from College next week, so I'll have more time to allow for development - - - I think!

These were some really incredible entries - I'm always amazed at what you guys accomplish without media - I'm still "crutching" it with media editors (but it makes my classroom presentations much more entertaining!), but just trying to meet the challenges has helped me learn LOADs!

An ecosystem challenge - you mean like a "food-chain" sort of thing? Or maybe a "SimEarth (tm)" type of thing? What's on your mind - this type of challenge can get pretty detailed amazingly quickly. And my creative waterfall just kicked in...

Can't wait to see what the next challenge brings. (sniff!) I've got a folder marked "DBPro Challenges" that holds four unfinished projects - complete with the competitor's projects! Gad, but I have soooo much to learn!!!

LongFist
SuperProgrammer For Hire

Back to top

Profile PM Email Website

da power pwnerer

19

Years of Service

User Offline

Joined: 28th Jul 2006

Location: Pittsburgh, PA

Posted: 12th Jul 2007 18:31 Edited at: 12th Jul 2007 18:34

Link

how about an 3d rpg game comp

with a battle system
money system magic system and all that good stuff

Back to top

Profile PM Email

Spopovich

20

Years of Service

User Offline

Joined: 12th Jul 2006

Location:

Posted: 12th Jul 2007 19:04

Link

How about Tetris, or just a puzzle game in general.

68% of rupublicans don't beleive in evolution but, only 5% of monkeys believe in republicans. ~Stephen Colbert 2007

Back to top

Profile PM

Cave Man

19

Years of Service

User Offline

Joined: 22nd Aug 2006

Location: North Carolina, US

Posted: 12th Jul 2007 19:33

Link

Quote: "you mean like a "food-chain" sort of thing?"

Yea, thats what i ment. It would be a good AI challenge with some of the animals hunting, and the other animals running away. I hope gadget freak picks this.

Back to top

Profile PM

qwe

22

Years of Service

User Offline

Joined: 3rd Sep 2003

Location: place

Posted: 12th Jul 2007 20:16

Link

>>let's brainstorm science ideas. physical simulations or Spelunx-like rooms that would be useful in the classroom and fun to play around with (learning in disguise)?>>

ecosystem challenge works^^

Back to top

Profile PM Email

Pricey

23

Years of Service

User Offline

Joined: 22nd Feb 2003

Location:

Posted: 12th Jul 2007 20:22

Link

i like the sound of that!

http://forum.thegamecreators.com/?m=forum_view&t=97191&b=32
God Of Thunder DEMO

Back to top

Profile PM Email

da power pwnerer

19

Years of Service

User Offline

Joined: 28th Jul 2006

Location: Pittsburgh, PA

Posted: 12th Jul 2007 20:59

Link

or my rpg idea

Back to top

Profile PM Email

acelepage

22

Years of Service

User Offline

Joined: 2nd Jun 2004

Location: Plattsville

Posted: 12th Jul 2007 23:25

Link

So, what's the verdict? Are you setting the challenge?
Ecosystem sounds pretty open ended. When I think of ecosystem, I think of the water cycle (rain, stream, river, ocean, evaporation, clouds).
What form of program? Game, simulation, 2D, 3D, whatever?

( 2b || !2b ), that is the question. The answer: true

Back to top

Profile PM Email

McLaine

20

Years of Service

User Offline

Joined: 20th Feb 2006

Location:

Posted: 13th Jul 2007 01:41

Link

Hey guys,

I'm back.

It's not my fault!

Back to top

Profile PM Email

da power pwnerer

19

Years of Service

User Offline

Joined: 28th Jul 2006

Location: Pittsburgh, PA

Posted: 13th Jul 2007 04:26 Edited at: 13th Jul 2007 04:26

Link

wahts up

(You probably don't recongnize me. I'm one of the FPSCreator users, but i'm considering buying DBPro, im currently using the demo. But don't worry, i'm not a total n00b. I'm a fast learner and i've used the trial on my other PC. It's amazing what you can get done in 30 days

Well then, back on track

)welcome back

Back to top

Profile PM Email

Ankillito

19

Years of Service

User Offline

Joined: 10th Dec 2006

Location: Litigious California

Posted: 13th Jul 2007 06:30

Link

what about just a sim type game in general. To vague?

I like the puzzle idea - i just might find time to enter something for that...

http://forum.thegamecreators.com/xt/xt_apollo_pic.php?i=1234759

Back to top

Profile PM

gadgetfreak 66

19

Years of Service

User Offline

Joined: 26th Sep 2006

Location: Kulin,WA, Australia

Posted: 13th Jul 2007 07:50

Link

RPG, or ecosystem sim...both sound interesting. An RPG would be fun, but making an RPG would be hell without media, so I think i might go with the ecosystem thing.

Challenge: Ecosystem Sim

This can be anything that shows a relationship between different things in nature. There are a few different things you could do, but I think the most fun ones to do would be mass element sims, like the growth of a forest, herd AI, seasonal migration, or some land where you can watch some herbivore block creatures eat grass, some carnivore triangle creatures eat the herbivores, and the carcases and manure help grow the grass, or something like that. Choose something thats fun or interesting to watch, anyway.

End of Challenge: Friday 27th of July

Back to top

Profile PM Email

qwe

22

Years of Service

User Offline

Joined: 3rd Sep 2003

Location: place

Posted: 13th Jul 2007 07:59 Edited at: 13th Jul 2007 08:00

Link

<edit>

Back to top

Profile PM Email

da power pwnerer

19

Years of Service

User Offline

Joined: 28th Jul 2006

Location: Pittsburgh, PA

Posted: 13th Jul 2007 08:50

Link

awwwww man. Oh well, there's always next time

I might actually try this challenge

http://Freewebs.com/noobisoft
Come to Noobisoft's website today!

Back to top

Profile PM Email

Ric

22

Years of Service

User Offline

Joined: 11th Jul 2004

Location: object position x

Posted: 13th Jul 2007 10:47

Link

Congrats Gadgetfreak on the last challenge - it was a great entry. This new challenge sounds very interesting. Could be quite educational too!

Back to top

Profile PM Email Website

Ankillito

19

Years of Service

User Offline

Joined: 10th Dec 2006

Location: Litigious California

Posted: 14th Jul 2007 00:03

Link

grrr...... I can't come up with any cool ideas!

http://forum.thegamecreators.com/xt/xt_apollo_pic.php?i=1234759

Back to top

Profile PM

Jimmy

22

Years of Service

User Offline

Joined: 20th Aug 2003

Location: Back in the USA

Posted: 15th Jul 2007 01:01 Edited at: 15th Jul 2007 01:03

Link

Quote: "interesting to watch, anyway."

Ok, maybe a little description is necessary. This little gem is called Eco Destructo.

They dance. Sweetly.

+ Code Snippet

sync on : sync rate 60
backdrop on : color backdrop rgb(132, 135, 140)
fog color rgb(28, 73, 140)
set camera range 10, 48000
fog distance 48000
set text font "arial" : set text size 12 : set text transparent
global Looptime : Global lastTime : UpdateDelta()
autocam off

`start coding
randomize timer()

for i = 1 to 50
	make object sphere i, 1+rnd(5)
	color object i, rgb(0,255,0)
	position object i, -100+rnd(200), 0, -100+rnd(200)
next i

position camera 0, 150, 0
point camera 0, 0, 0

do
	UpdateDelta()
	
	EcoDance()
	
	text 10, 10, str$(screen fps())
	sync
loop

end

Function EcoDance()
	for i = 1 to 50
		yrotate object i, rnd(360)
		move object i, Delta(2+rnd(5))
	next i
endfunction

`Delta
Function UpdateDelta()
TheTimer=Timer()
LoopTime=TheTimer-lastTime
lastTime=TheTimer
Endfunction

Function Delta(Value1#)
 	ResultValue#=Value1#*0.060*Looptime
Endfunction ResultValue#

Bloggy

Back to top

Profile PM Website

gadgetfreak 66

19

Years of Service

User Offline

Joined: 26th Sep 2006

Location: Kulin,WA, Australia

Posted: 15th Jul 2007 13:59

Link

err.. is this an entry? If so, what's the relationship it's demonstrating?

Back to top

Profile PM Email

Kohaku

22

Years of Service

User Offline

Joined: 3rd May 2004

Location: The not very United Kingdom

Posted: 15th Jul 2007 15:06 Edited at: 15th Jul 2007 15:10

Link

That is an excellent display of ecosystemology.

You are not alone.

Back to top

Profile PM Website

NA170425

User Deleted

Posted: 15th Jul 2007 19:54

Link

Can you guess what it simulates?

REM Project: Planets
REM Created: 7/15/2007 9:43:47 AM
REM
REM ***** Main Source File *****
REM

#constant Planets = 20
#constant pi = 3.14159265

sync on
sync rate 60
autocam off
hide mouse
color backdrop 0
randomize timer()

type Planet

object as dword
  radius as float
  mass as float

vx as float
  vy as float
  vz as float

endtype

dim Planet() as Planet

local count as dword
local id as dword

make matrix 1, 300, 300, 30, 30
position matrix 1, -100, 0, -100
update matrix 1

for count = 0 to Planets
  id = makePlanet()
next count

position camera -200, 100, -200
point camera 0, 0, 0

local dx as float
local dy as float
local dz as float
local total as float
local g as double float
local gravityconstant as double float

gravityconstant = 0.0000001

do

text 0, 0, "Screen Fps: " + str$(screen fps())
  text 0, 20, "Universal Gravity Constant: " + str$(gravityconstant*10000000, 1) + " * 10^-7"
  text 0, 30, "(arrowkeys to change)"
  if upkey() then inc gravityconstant, 0.0000000015
  if downkey() then dec gravityconstant, 0.0000000015
  text 0, 50, "Press enter to create more planets"
  if returnkey() then makePlanet()
  text 0, 60, "Press delete to erase all the planets"
  if keystate(211) then deletePlanets()
  text 0, 70, "Total Planets: " + str$(array count(Planet())+1)

for Planet1 = 0 to array count(Planet())
    for Planet2 = 0 to array count(Planet())
      if Planet1 <> Planet2

dx = object position x(Planet(Planet2).object) - object position x(Planet(Planet1).object)
        dy = object position y(Planet(Planet2).object) - object position y(Planet(Planet1).object)
        dz = object position x(Planet(Planet2).object) - object position z(Planet(Planet1).object)
        total = sqrt(dx^2 + dy^2 + dz^2)
        dx = dx / total
        dy = dy / total
        dz = dz / total

g = (Planet(Planet1).mass * Planet(Planet2).mass * gravityconstant) / total

inc Planet(Planet1).vx, dx*g
        inc Planet(Planet1).vy, dy*g
        inc Planet(Planet1).vz, dz*g

endif
    next Planet2
    
    position object Planet(Planet1).object, object position x(Planet(Planet1).object) + Planet(Planet1).vx, object position y(Planet(Planet1).object) + Planet(Planet1).vy, object position z(Planet(Planet1).object) + Planet(Planet1).vz
    
  next Planet1
  
  sync
  
loop

function makePlanet()

local id as dword

add to queue Planet()
  id = array count(Planet())
  Planet().radius = rnd(9) + 1
  Planet().mass = (4/3) * pi * Planet().radius ^ 3
  Planet().object = make_object_sphere(Planet().radius)
  color object Planet().object, rnd(0x00FFFFFF)
  position object Planet(id).object, (rnd(200) - 100), (rnd(200) - 100), (rnd(200) - 100)
  
endfunction id

function deletePlanets()

local totalplanets as integer
  local id as integer

totalplanets = array count(Planet())
  
  while id <= totalplanets
    delete object Planet(id).object
    inc id
  endwhile
  
  empty array Planet()
  
endfunction

function make_object_sphere(size as float)

local id as dword

id = free_object()
  make object sphere id, size

endfunction id

function free_object()

local id as dword = 1

while object exist(id)
    inc id
  endwhile

endfunction id

+ Code Snippet

REM Project: Planets
REM Created: 7/15/2007 9:43:47 AM
REM
REM ***** Main Source File *****
REM

#constant Planets = 20
#constant pi = 3.14159265

sync on
sync rate 60
autocam off
hide mouse
color backdrop 0
randomize timer()

type Planet

  object as dword
  radius as float
  mass as float

  vx as float
  vy as float
  vz as float

endtype

dim Planet() as Planet

local count as dword
local id as dword

make matrix 1, 300, 300, 30, 30
position matrix 1, -100, 0, -100
update matrix 1

for count = 0 to Planets
  id = makePlanet()
next count

position camera -200, 100, -200
point camera 0, 0, 0

local dx as float
local dy as float
local dz as float
local total as float
local g as double float
local gravityconstant as double float

gravityconstant = 0.0000001

do

  text 0, 0, "Screen Fps: " + str$(screen fps())
  text 0, 20, "Universal Gravity Constant: " + str$(gravityconstant*10000000, 1) + " * 10^-7"
  text 0, 30, "(arrowkeys to change)"
  if upkey() then inc gravityconstant, 0.0000000015
  if downkey() then dec gravityconstant, 0.0000000015
  text 0, 50, "Press enter to create more planets"
  if returnkey() then makePlanet()
  text 0, 60, "Press delete to erase all the planets"
  if keystate(211) then deletePlanets()
  text 0, 70, "Total Planets: " + str$(array count(Planet())+1)

  for Planet1 = 0 to array count(Planet())
    for Planet2 = 0 to array count(Planet())
      if Planet1 <> Planet2

        dx = object position x(Planet(Planet2).object) - object position x(Planet(Planet1).object)
        dy = object position y(Planet(Planet2).object) - object position y(Planet(Planet1).object)
        dz = object position x(Planet(Planet2).object) - object position z(Planet(Planet1).object)
        total = sqrt(dx^2 + dy^2 + dz^2)
        dx = dx / total
        dy = dy / total
        dz = dz / total

        g = (Planet(Planet1).mass * Planet(Planet2).mass * gravityconstant) / total

        inc Planet(Planet1).vx, dx*g
        inc Planet(Planet1).vy, dy*g
        inc Planet(Planet1).vz, dz*g

      endif
    next Planet2
    
    position object Planet(Planet1).object, object position x(Planet(Planet1).object) + Planet(Planet1).vx, object position y(Planet(Planet1).object) + Planet(Planet1).vy, object position z(Planet(Planet1).object) + Planet(Planet1).vz
    
  next Planet1
  
  sync
  
loop

function makePlanet()

  local id as dword

  add to queue Planet()
  id = array count(Planet())
  Planet().radius = rnd(9) + 1
  Planet().mass = (4/3) * pi * Planet().radius ^ 3
  Planet().object = make_object_sphere(Planet().radius)
  color object Planet().object, rnd(0x00FFFFFF)
  position object Planet(id).object, (rnd(200) - 100), (rnd(200) - 100), (rnd(200) - 100)
  
endfunction id

function deletePlanets()

  local totalplanets as integer
  local id as integer

  totalplanets = array count(Planet())
  
  while id <= totalplanets
    delete object Planet(id).object
    inc id
  endwhile
  
  empty array Planet()
  
endfunction

function make_object_sphere(size as float)

  local id as dword

  id = free_object()
  make object sphere id, size

endfunction id

function free_object()

  local id as dword = 1

  while object exist(id)
    inc id
  endwhile

endfunction id

This sig has been viewed

times.

Back to top

Profile PM Email Website

Pricey

23

Years of Service

User Offline

Joined: 22nd Feb 2003

Location:

Posted: 15th Jul 2007 20:17

Link

i don't know what it simulates
but it stimulates me greatly

i think i'll be entering this one!
watch this space yo

http://forum.thegamecreators.com/?m=forum_view&t=97191&b=32
God Of Thunder DEMO

Back to top

Profile PM Email

gadgetfreak 66

19

Years of Service

User Offline

Joined: 26th Sep 2006

Location: Kulin,WA, Australia

Posted: 16th Jul 2007 07:50

Link

@ CodeDragon:

Thats pretty cool, I spent ages watching it, but the group always tends to move off the screen, or just get really far away, so an auto-camera or camera controls would be nice.

Back to top

Profile PM Email

NA170425

User Deleted

Posted: 16th Jul 2007 14:28

Link

I have no idea why it does that. I added a big planet that never moves, hopefully the gravitation force of it will stop them from wandering. I also added some basic camera controls, autocam on the entire group would be hard I think...

You can change the distance exponent now too. In real life it's 2.00 but for some reason that value causes the program to sliently crash, so I kept it as 1.

REM Project: Planets
REM Created: 7/15/2007 9:43:47 AM
REM
REM ***** Main Source File *****
REM

#constant Planets = 20
#constant pi = 3.14159265

sync on
sync rate 60
autocam off
hide mouse
color backdrop 0
randomize timer()

type Planet

object as dword
  radius as double float
  mass as double float

vx as double float
  vy as double float
  vz as double float

endtype

dim Planet() as Planet

local count as dword
local id as dword

make matrix 1, 300, 300, 30, 30
position matrix 1, -150, 0, -150
update matrix 1

for count = 0 to Planets
  id = makePlanet()
next count

local dx as double float
local dy as double float
local dz as double float
local total as double float
local g as double float
local gravityconstant as double float
local camang as float
local camheight as float
local disexp as double float

gravityconstant = 0.0000001
disexp = 1

do

text 0, 0, "Screen Fps: " + str$(screen fps())
  text 0, 20, "Universal Gravity Constant: " + str$(gravityconstant*10000000, 1) + " * 10^-7"
  text 0, 30, "(arrowkeys to change)"
  if upkey() then inc gravityconstant, 0.0000000015
  if downkey() then dec gravityconstant, 0.0000000015
  text 0, 40, "Distance Exponent: " + str$(disexp, 2)
  text 0, 50, "(left and right arrowkeys to change)"
  if rightkey() then inc disexp, 0.01
  if leftkey() then dec disexp, 0.01
  text 0, 70, "Press enter to create more planets"
  if returnkey() then makePlanet()
  text 0, 80, "Press delete to erase all the planets"
  if keystate(211) then deletePlanets()
  text 0, 90, "Total Planets: " + str$(array count(Planet())+1)

for Planet1 = 1 to array count(Planet())
    for Planet2 = 0 to array count(Planet())
      if Planet1 <> Planet2

dx = object position x(Planet(Planet2).object) - object position x(Planet(Planet1).object)
        dy = object position y(Planet(Planet2).object) - object position y(Planet(Planet1).object)
        dz = object position x(Planet(Planet2).object) - object position z(Planet(Planet1).object)
        total = sqrt(dx^2 + dy^2 + dz^2)

dx = dx / total
        dy = dy / total
        dz = dz / total

g = (Planet(Planet1).mass * Planet(Planet2).mass * gravityconstant) / total^disexp

inc Planet(Planet1).vx, dx*g
        inc Planet(Planet1).vy, dy*g
        inc Planet(Planet1).vz, dz*g

endif
    next Planet2
    
    position object Planet(Planet1).object, object position x(Planet(Planet1).object) + Planet(Planet1).vx, object position y(Planet(Planet1).object) + Planet(Planet1).vy, object position z(Planet(Planet1).object) + Planet(Planet1).vz
    
  next Planet1
  
  inc camang, -mousemovex()/2.0
  inc camheight, -mousemovey()/2.0
  set camera to follow 0, 0, 0, camang, 250, camheight, 1, 0
  point camera 0, 0, 0
  
  sync
  
loop

function makePlanet()

local id as dword

add to queue Planet()
  id = array count(Planet())
  Planet().radius = rnd(9) + 1
  if id = 0 then Planet().radius = 20
  Planet().mass = (4/3) * pi * Planet().radius ^ 3
  Planet().object = make_object_sphere(Planet().radius)
  color object Planet().object, rnd(0x00FFFFFF)
  if id > 0 then position object Planet(id).object, (rnd(200) - 100), (rnd(200) - 100), (rnd(200) - 100)
  
endfunction id

function deletePlanets()

local totalplanets as integer
  local id as integer

totalplanets = array count(Planet())
  
  while id <= totalplanets
    delete object Planet(id).object
    inc id
  endwhile
  
  empty array Planet()
  
endfunction

function make_object_sphere(size as float)

local id as dword

id = free_object()
  make object sphere id, size

endfunction id

function free_object()

local id as dword = 1

while object exist(id)
    inc id
  endwhile

endfunction id

+ Code Snippet

REM Project: Planets
REM Created: 7/15/2007 9:43:47 AM
REM
REM ***** Main Source File *****
REM

#constant Planets = 20
#constant pi = 3.14159265

sync on
sync rate 60
autocam off
hide mouse
color backdrop 0
randomize timer()

type Planet

  object as dword
  radius as double float
  mass as double float

  vx as double float
  vy as double float
  vz as double float

endtype

dim Planet() as Planet

local count as dword
local id as dword

make matrix 1, 300, 300, 30, 30
position matrix 1, -150, 0, -150
update matrix 1

for count = 0 to Planets
  id = makePlanet()
next count

local dx as double float
local dy as double float
local dz as double float
local total as double float
local g as double float
local gravityconstant as double float
local camang as float
local camheight as float
local disexp as double float

gravityconstant = 0.0000001
disexp = 1

do

  text 0, 0, "Screen Fps: " + str$(screen fps())
  text 0, 20, "Universal Gravity Constant: " + str$(gravityconstant*10000000, 1) + " * 10^-7"
  text 0, 30, "(arrowkeys to change)"
  if upkey() then inc gravityconstant, 0.0000000015
  if downkey() then dec gravityconstant, 0.0000000015
  text 0, 40, "Distance Exponent: " + str$(disexp, 2)
  text 0, 50, "(left and right arrowkeys to change)"
  if rightkey() then inc disexp, 0.01
  if leftkey() then dec disexp, 0.01
  text 0, 70, "Press enter to create more planets"
  if returnkey() then makePlanet()
  text 0, 80, "Press delete to erase all the planets"
  if keystate(211) then deletePlanets()
  text 0, 90, "Total Planets: " + str$(array count(Planet())+1)

  for Planet1 = 1 to array count(Planet())
    for Planet2 = 0 to array count(Planet())
      if Planet1 <> Planet2

        dx = object position x(Planet(Planet2).object) - object position x(Planet(Planet1).object)
        dy = object position y(Planet(Planet2).object) - object position y(Planet(Planet1).object)
        dz = object position x(Planet(Planet2).object) - object position z(Planet(Planet1).object)
        total = sqrt(dx^2 + dy^2 + dz^2)

        dx = dx / total
        dy = dy / total
        dz = dz / total

        g = (Planet(Planet1).mass * Planet(Planet2).mass * gravityconstant) / total^disexp

        inc Planet(Planet1).vx, dx*g
        inc Planet(Planet1).vy, dy*g
        inc Planet(Planet1).vz, dz*g

      endif
    next Planet2
    
    position object Planet(Planet1).object, object position x(Planet(Planet1).object) + Planet(Planet1).vx, object position y(Planet(Planet1).object) + Planet(Planet1).vy, object position z(Planet(Planet1).object) + Planet(Planet1).vz
    
  next Planet1
  
  inc camang, -mousemovex()/2.0
  inc camheight, -mousemovey()/2.0
  set camera to follow 0, 0, 0, camang, 250, camheight, 1, 0
  point camera 0, 0, 0
  
  sync
  
loop

function makePlanet()

  local id as dword

  add to queue Planet()
  id = array count(Planet())
  Planet().radius = rnd(9) + 1
  if id = 0 then Planet().radius = 20
  Planet().mass = (4/3) * pi * Planet().radius ^ 3
  Planet().object = make_object_sphere(Planet().radius)
  color object Planet().object, rnd(0x00FFFFFF)
  if id > 0 then position object Planet(id).object, (rnd(200) - 100), (rnd(200) - 100), (rnd(200) - 100)
  
endfunction id

function deletePlanets()

  local totalplanets as integer
  local id as integer

  totalplanets = array count(Planet())
  
  while id <= totalplanets
    delete object Planet(id).object
    inc id
  endwhile
  
  empty array Planet()
  
endfunction

function make_object_sphere(size as float)

  local id as dword

  id = free_object()
  make object sphere id, size

endfunction id

function free_object()

  local id as dword = 1

  while object exist(id)
    inc id
  endwhile

endfunction id

This sig has been viewed

times.

Back to top

Profile PM Email Website

Ric

22

Years of Service

User Offline

Joined: 11th Jul 2004

Location: object position x

Posted: 16th Jul 2007 19:11

Link

Not quite sure how it relates to ecology, but it covers the 'interesting' requirement! Crashes after about 30s here - most likely you have values somewhere which are exploding out of bounds.

Back to top

Profile PM Email Website

LongFist

22

Years of Service

User Offline

Joined: 25th May 2004

Location: Winston-Salem, NC

Posted: 17th Jul 2007 23:45 Edited at: 17th Jul 2007 23:49

Link

Hey GadgetFreak: What're we gonna' do? (Or did I miss that?)

LongFist
SuperProgrammer For Hire

Back to top

Profile PM Email Website

qwe

22

Years of Service

User Offline

Joined: 3rd Sep 2003

Location: place

Posted: 18th Jul 2007 01:20 Edited at: 18th Jul 2007 01:30

Link

cool code dragon. i programmed something very similar a while ago, this also simulates gravity

REM Particles lab
REM By qwe

REM Constants
global maxp = 20
global border as double float				:	border = 1000.0
global Gravity as double float			:	Gravity = 1.0
global Electrostatic as double float	:	Electrostatic = 1.0

REM initialize
GOSUB INI

REM Test
`remstart
for a = 0 to 16
	p(a).t = 1
	p(a).cr = rnd(255)
	p(a).cg = rnd(255)
	p(a).cb = rnd(255)
next a

p(0).ir = 10.0
p(0).im = 1.0

REM Inner

p(1).ir = 1.0
p(1).im = 0.000001
p(1).ix = 10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(1).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * 0.7071

p(2).ir = 1.0
p(2).im = 0.000001
p(2).ix = -10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(2).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * -0.7071

p(3).ir = 1.0
p(3).im = 0.000001
p(3).iz = 10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(3).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * 0.7071

p(4).ir = 1.0
p(4).im = 0.000001
p(4).iz = -10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(4).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * -0.7071

REM Outer

p(5).ir = 1.0
p(5).im = 0.000001
p(5).ix = 20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(5).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * 0.7071

p(6).ir = 1.0
p(6).im = 0.000001
p(6).ix = -20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(6).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * -0.7071

p(7).ir = 1.0
p(7).im = 0.000001
p(7).iz = 20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(7).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * 0.7071

p(8).ir = 1.0
p(8).im = 0.000001
p(8).iz = -20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(8).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * -0.7071

REM OuterOuter

p(9).ir = 1.0
p(9).im = 0.000001
p(9).ix = 30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(9).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * 0.7071

p(10).ir = 1.0
p(10).im = 0.000001
p(10).ix = -30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(10).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * -0.7071

p(11).ir = 1.0
p(11).im = 0.000001
p(11).iz = 30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(11).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * 0.7071

p(12).ir = 1.0
p(12).im = 0.000001
p(12).iz = -30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(12).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * -0.7071

REM OuterOuterOuter

p(13).ir = 1.0
p(13).im = 0.000001
p(13).ix = 40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(13).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * 0.7071

p(14).ir = 1.0
p(14).im = 0.000001
p(14).ix = -40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(14).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * -0.7071

p(15).ir = 1.0
p(15).im = 0.000001
p(15).iz = 40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(15).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * 0.7071

p(16).ir = 1.0
p(16).im = 0.000001
p(16).iz = -40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(16).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * -0.7071

`remend

remstart
for a = 0 to maxp
	p(a).t = 1
	p(a).cr = rnd(255)
	p(a).cg = rnd(255)
	p(a).cb = rnd(255)
	p(a).ir = 10.0
	p(a).im = 0.01
	p(a).ix = rnd(int(border)*2)-int(border)
	p(a).iy = rnd(int(border)*2)-int(border)
	p(a).iz = rnd(int(border)*2)-int(border)
	
	p(a).ixv = (rnd(2)-1) * ((rnd(100) + 0.0)/1000.0)
	p(a).iyv = (rnd(2)-1) * ((rnd(100) + 0.0)/1000.0)
	p(a).izv = (rnd(2)-1) * ((rnd(100) + 0.0)/1000.0)
next a
remend
remstart
for a = 0 to 1
	p(a).t = 1
	p(a).cr = rnd(255)
	p(a).cg = rnd(255)
	p(a).cb = rnd(255)
	p(a).ir = 0.3333333
	p(a).im = 0.01
next a
p(0).ix = 5.0
remend

REM Main
display_readouts = 0
reini()
do
	REM Display data
	set cursor 0,0
	ink rgb(255,255,255),0
	print "FPS: "; screen fps()
	print "Time elapsed: "; left$(str$(time), 5); " seconds"
	print "Total energy: "; total_energy
	print
	if display_readouts = 1 then show_readouts()
	
	REM User control
	do_move()
	do_look()
	do_console()
	
	REM run particles
	run_particles()
	
   sync
loop

function show_readouts()
	for a = 0 to maxp
		if p(a).t <> 0 and a < 3
			ink rgb(p(a).cr, p(a).cg, p(a).cb), 0
			print "X Position of particle ";a;": "; p(a).x
			print "Y Position of particle ";a;": "; p(a).y
			print "Z Position of particle ";a;": "; p(a).z
			print
			print "X Velocity of particle ";a;": "; p(a).xv
			print "Y Velocity of particle ";a;": "; p(a).yv
			print "Z Velocity of particle ";a;": "; p(a).zv
			print
			print "Mass of particle ";a;": "; p(a).m
			print "Charge of particle ";a;": "; p(a).c
			print "Radius of particle ";a;": "; p(a).r
			print
		endif
	next a
endfunction

function apply_test_keys(a)
	REM Test
	if keystate(2)
		p(a).xv = p(a).xv * 1.01
		p(a).yv = p(a).yv * 1.01
		p(a).zv = p(a).zv * 1.01
	endif
	if keystate(3)
		p(a).xv = p(a).xv * 0.99
		p(a).yv = p(a).yv * 0.99
		p(a).zv = p(a).zv * 0.99
	endif          
	if keystate(4)
		dec p(a).yv, 0.01
	endif
	if keystate(5)
		inc p(a).m, 0.0001
	endif
	if keystate(6)
		p(a).xv = p(a).xv * 0.0
		p(a).yv = p(a).yv * 0.0
		p(a).zv = p(a).zv * 0.0
	endif
endfunction

function run_particles()
	inc time, 0.01
	
	for a = 0 to maxp
		if p(a).t <> 0
			REM apply forces and collision with border
			
			if p(a).m <> 0.0 then apply_gravity(a)
			if p(a).c <> 0.0 then apply_electromagnetism(a)
			
			apply_collision_walls(a)
			apply_test_keys(a)
			
		endif
	next a

for a = 0 to maxp
		if p(a).t <> 0
			REM apply velocity and update 3d representation
			
			update_particle_positions(a)
			update_particle_objects(a)
			
		endif
	next a
	
	total_energy = 0.0
	for a = 0 to maxp
		if p(a).t <> 0
			REM calculate total energy of system as a means of error checking
			calculate_total_energy(a)
			inc total_energy, p(a).te
		endif
	next a
endfunction

function calculate_total_energy(a)
	` KE PE TE
	p(a).pe = 0.0
	p(a).ke = 0.0
	p(a).te = 0.0
	
	REM calculate KE (Kinetic Energy)
	REM KE = 1/2 * M * V^2
	p(a).ke = 0.5 * p(a).m * (sqrt(p(a).xv^2 + p(a).yv^2 + p(a).zv^2) ^ 2.0)
	
	REM Calculate PE (Potential energy)
	for b = 0 to maxp
		if p(b).t <> 0 and  b <> a `for each active particle besides a
			
			REM calculate distance
			dx = p(b).x - p(a).x
			dy = p(b).y - p(a).y
			dz = p(b).z - p(a).z
			Distance = sqrt(dx*dx + dy*dy + dz*dz)
			REM calculate potential for gravity
			REM V = G * M1 * M2 / D
			inc p(a).pe, (Gravity * p(a).m * p(b).m) /  Distance
			
		endif
	next b
	
	REM Calculate TE (Total energy)
	p(a).te = p(a).ke + p(a).pe
	
endfunction

function apply_gravity(a)
	local Force as double float
	local Distance as double float
	local dx as double float
	local dy as double float
	local dz as double float
	local Acceleration as double float
			
	for b = 0 to maxp
		if p(b).t <> 0 and p(b).m <> 0.0 and b <> a `for each active particle besides particle a with a mass
			
			REM calculate distance
			dx = p(b).x - p(a).x
			dy = p(b).y - p(a).y
			dz = p(b).z - p(a).z
			Distance = sqrt(dx*dx + dy*dy + dz*dz)
			REM calculate force
			`F = (G*m1*m2)/(DD)
			Force = (Gravity * p(a).m * p(b).m) / (Distance * Distance)
			REM calculate acceleration
			`A = F / M
			Acceleration = Force / p(a).m
			REM create unit vector from origin
			REM with magnitude equal to appropriate acceleratoin
			REM and direction pointing to particle B's relative position
			`we divide each value by D so that it is a unit vector
			`and multiply by acceleration so its magnitude equals acceleration
			inc p(a).xv, dx * (Acceleration/Distance)
			inc p(a).yv, dy * (Acceleration/Distance)
			inc p(a).zv, dz * (Acceleration/Distance)
			
		endif
	next b
			
endfunction

function apply_electromagnetism(a)
	local Force as double float
	local Distance as double float
	local dx as double float
	local dy as double float
	local dz as double float
	local Acceleration as double float
			
	for b = 0 to maxp
		if p(b).t <> 0 and p(b).c <> 0.0 and  b <> a `for each active particle besides particle a with a charge
			
			REM calculate distance
			dx = p(b).x - p(a).x
			dy = p(b).y - p(a).y
			dz = p(b).z - p(a).z
			Distance = sqrt(dx*dx + dy*dy + dz*dz)
			REM calculate force (electromagnetism)
			`F = (Kc*C1*C2)/(DD)
			Force = (Electrostatic * abs(p(a).c) * abs(p(b).c)) / (Distance * Distance)
			REM calculate acceleration
			`A = F / M
			Acceleration = Force` / p(a).m
			REM whether A accelerates away or towards B depends on likeness of charge
			if p(a).c > 0.0	`charge on A and B is same (positive)
				if p(b).c > 0.0 then Acceleration = -1.0 * Acceleration
			endif
			if p(a).c < 0.0	`charge on A and B is same (negative)
				if p(b).c < 0.0 then Acceleration = -1.0 * Acceleration
			endif
			REM create unit vector from origin
			REM with magnitude equal to appropriate acceleratoin
			REM and direction pointing to particle B's relative position
			`we divide each value by D so that it is a unit vector
			`and multiply by acceleration so its magnitude equals acceleration
			inc p(a).xv, dx * (Acceleration/Distance)
			inc p(a).yv, dy * (Acceleration/Distance)
			inc p(a).zv, dz * (Acceleration/Distance)
			
		endif
	next b

endfunction

function do_console()
	REM toggle sync rate
	if keystate(74)
		sync rate 5
	endif
	if keystate(78)
		sync rate 0
	endif
	if keystate(156)
		sync rate 50
	endif
	REM sc and oldsc and console are global
	oldsc = sc
	sc = scancode()
	REM Toggle console with tab key
	if sc = 15
		if oldsc = 0
			if console = 1
				console = 0
			else
				console = 1
				clear entry buffer
				console$(15) = ""
			endif
		endif
	endif
	REM Show console and type in console
	if console
		ink rgb(32,32,32),0
		box 0,0,sw,300
		set cursor 0,0
		for a=0 to 15
			ink rgb(64,255,64),0
			if left$(console$(a), 3) = ">> " then ink rgb(255,64,64),0
			print console$(a)
		next a
		
		olda$  = a$
		a$ = entry$()
		clear entry buffer
		if 9 = asc(a$) then a$="q"
		if 13 = asc(a$) then a$=""
		if 28 = oldsc and sc = 0 then do_console_command()
		if 8 = asc(a$)
		   a# = len(console$(15))
		   b$ = left$(console$(15), a# - 1)
		   console$(15) = b$
		else
		   console$(15) = console$(15)+a$
		endif
	endif
endfunction

function shift_console_items_up()
	for a = 0 to 14
		console$(a) = console$(a+1)
	next a
	console$(15) = ""
endfunction

function do_console_command()
	shift_console_items_up()
	c$ = console$(14)
	REM DO COMMANDS
	if c$ = "ini"
		reini()
	endif
	if c$ = "spd 0"
		sync rate 0
	endif
	if c$ = "spd 5"
		sync rate 5
	endif
	if c$ = "spd 50"
		sync rate 50
	endif
	if c$ = "read 0"
		display_readouts = 0
	endif
	if c$ = "read 1"
		display_readouts = 1
	endif
endfunction

INI:

REM Types
type particles_kinds
	name as string
	m as double float
	c as double float
	r as double float
	cr
	cg
	cb
	t
endtype
type particles
	x as double float
	y as double float
	z as double float
	xv as double float
	yv as double float
	zv as double float
	ix as double float
	iy as double float
	iz as double float
	ixv as double float
	iyv as double float
	izv as double float
	im as double float
	ic as double float
	ir as double float
	m as double float
	c as double float
	r as double float
	ke as double float
	pe as double float
	te as double float
	cr
	cg
	cb
	t
endtype

REM Display
sync on
sync rate 0
global sw	`screen width
global sh	`screen height
global display_readouts = 1
sw = screen width()
sh = screen height()

REM initiate console
global sc		`Scan code (for console)
global oldsc	`Previous scan code (for console)
global olda$	`previous entry (for console)
global console	`console toggle

DIM console$(16)
for a=0 to 15
	console$(a)=" "
next a

REM initiate globals
global time as double float
global total_energy as double float

REM initiate particles
DIM pk(32) as particles_kinds
DIM p(maxp+1) as particles

REM initiate camera
make camera 1
backdrop on 1
color backdrop 1, rgb(0,0,0)
Set Camera Range 1, 2000
position camera 1,0,50,50
point camera 1,0,0,0

REM initiate level
`load object "level/level.x",1
make object cube 2, border * 2.0
set object wireframe 2, 1
set object cull 2, 0
color object 2, rgb(0,255,0)
make light 1
set point light 1, 0, 0, 0
color light 1, 255, 255, 255
set light range 1, 100000

return

function update_particle_objects(a)
	position object a+100, p(a).x, p(a).y, p(a).z
	point object a+100, camera position x(1), camera position y(1), camera position z(1)
endfunction

function update_particle_positions(a)
	inc p(a).x, p(a).xv
	inc p(a).y, p(a).yv
	inc p(a).z, p(a).zv
endfunction

function reini()
	time = 0
	
	for a = 0 to maxp
		if object exist(a+100) then delete object a+100
		
		if p(a).t <> 0
			REM set initial properties
			p(a).x = p(a).ix
			p(a).y = p(a).iy
			p(a).z = p(a).iz
			p(a).xv = p(a).ixv
			p(a).yv = p(a).iyv
			p(a).zv = p(a).izv
			p(a).m = p(a).im
			p(a).c = p(a).ic
			p(a).r = p(a).ir
			REM make object
			make object sphere a+100, p(a).r, 4, 4
			position object a+100, p(a).ix, p(a).iy, p(a).iz
			color object a+100, rgb(p(a).cr, p(a).cg, p(a).cb)
		endif
	next a
endfunction

function apply_collision_walls(a)	
			
	if p(a).x > border or p(a).x < -1*border
		p(a).xv = p(a).xv * -1.0
	endif
	if p(a).y > border or p(a).y < -1*border
		p(a).yv = p(a).yv * -1.0
	endif
	if p(a).z > border or p(a).z < -1*border
		p(a).zv = p(a).zv * -1.0
	endif
			
endfunction

function do_move()
	local spd as float
	
	if console = 0
		REM var
		spd=0.5
		
		REM set speed and if hold LShift then go twice as fast
		if keystate(42)
			spd=spd*10
		endif
		
		REM move around
		if keystate(17)
			move camera 1, spd
		endif
		if keystate(31)
			move camera 1, -spd
		endif
		if keystate(30)
			position camera 1, newxvalue(camera position x(1), camera angle y(1)-90.0, spd), camera position y(1), newzvalue(camera position z(1), camera angle y(1)-90.0, spd)
		endif
		if keystate(32)
			position camera 1, newxvalue(camera position x(1), camera angle y(1)+90.0, spd), camera position y(1), newzvalue(camera position z(1), camera angle y(1)+90.0, spd)
		endif
		
		REM push upward or downward
		if keystate(200) or keystate(57)
			position camera 1, camera position x(1), camera position y(1) + spd, camera position z(1)
		endif
		if keystate(208) or keystate(56)
			position camera 1, camera position x(1), camera position y(1) - spd, camera position z(1)
		endif
	endif
endfunction

function do_look()
	local mx as float
	local my as float
	local mousedamp as float
	local cam as integer
	
	REM var
	mousedamp = 8.0
	
	REM look around, unless right clicking
	if mouseclick() > 1
		show mouse
	else
		hide mouse
		`still cursor
		position mouse 500,500
		`dampen amount to move
		mx=mousemovex()/mousedamp
		my=mousemovey()/mousedamp
		`turn camera object
		xrotate camera 1, wrapvalue(camera angle x(1)+my)
		yrotate camera 1, wrapvalue(camera angle y(1)+mx)
		`bound camera object angle
		if camera angle x(1)>89 and camera angle x(1)<180 then xrotate camera 1,89
		if camera angle x(1)<271 and camera angle x(1)>180 then xrotate camera 1,271
	endif
endfunction

+ Code Snippet

REM Particles lab
REM By qwe

REM Constants
global maxp = 20
global border as double float				:	border = 1000.0
global Gravity as double float			:	Gravity = 1.0
global Electrostatic as double float	:	Electrostatic = 1.0

REM initialize
GOSUB INI

REM Test
`remstart
for a = 0 to 16
	p(a).t = 1
	p(a).cr = rnd(255)
	p(a).cg = rnd(255)
	p(a).cb = rnd(255)
next a

p(0).ir = 10.0
p(0).im = 1.0

REM Inner

p(1).ir = 1.0
p(1).im = 0.000001
p(1).ix = 10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(1).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * 0.7071

p(2).ir = 1.0
p(2).im = 0.000001
p(2).ix = -10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(2).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * -0.7071

p(3).ir = 1.0
p(3).im = 0.000001
p(3).iz = 10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(3).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * 0.7071

p(4).ir = 1.0
p(4).im = 0.000001
p(4).iz = -10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(4).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * -0.7071

REM Outer

p(5).ir = 1.0
p(5).im = 0.000001
p(5).ix = 20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(5).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * 0.7071

p(6).ir = 1.0
p(6).im = 0.000001
p(6).ix = -20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(6).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * -0.7071

p(7).ir = 1.0
p(7).im = 0.000001
p(7).iz = 20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(7).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * 0.7071

p(8).ir = 1.0
p(8).im = 0.000001
p(8).iz = -20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(8).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * -0.7071

REM OuterOuter

p(9).ir = 1.0
p(9).im = 0.000001
p(9).ix = 30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(9).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * 0.7071

p(10).ir = 1.0
p(10).im = 0.000001
p(10).ix = -30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(10).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * -0.7071

p(11).ir = 1.0
p(11).im = 0.000001
p(11).iz = 30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(11).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * 0.7071

p(12).ir = 1.0
p(12).im = 0.000001
p(12).iz = -30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(12).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * -0.7071

REM OuterOuterOuter

p(13).ir = 1.0
p(13).im = 0.000001
p(13).ix = 40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(13).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * 0.7071

p(14).ir = 1.0
p(14).im = 0.000001
p(14).ix = -40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(14).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * -0.7071

p(15).ir = 1.0
p(15).im = 0.000001
p(15).iz = 40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(15).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * 0.7071

p(16).ir = 1.0
p(16).im = 0.000001
p(16).iz = -40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(16).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * -0.7071

`remend

remstart
for a = 0 to maxp
	p(a).t = 1
	p(a).cr = rnd(255)
	p(a).cg = rnd(255)
	p(a).cb = rnd(255)
	p(a).ir = 10.0
	p(a).im = 0.01
	p(a).ix = rnd(int(border)*2)-int(border)
	p(a).iy = rnd(int(border)*2)-int(border)
	p(a).iz = rnd(int(border)*2)-int(border)
	
	p(a).ixv = (rnd(2)-1) * ((rnd(100) + 0.0)/1000.0)
	p(a).iyv = (rnd(2)-1) * ((rnd(100) + 0.0)/1000.0)
	p(a).izv = (rnd(2)-1) * ((rnd(100) + 0.0)/1000.0)
next a
remend
remstart
for a = 0 to 1
	p(a).t = 1
	p(a).cr = rnd(255)
	p(a).cg = rnd(255)
	p(a).cb = rnd(255)
	p(a).ir = 0.3333333
	p(a).im = 0.01
next a
p(0).ix = 5.0
remend

REM Main
display_readouts = 0
reini()
do
	REM Display data
	set cursor 0,0
	ink rgb(255,255,255),0
	print "FPS: "; screen fps()
	print "Time elapsed: "; left$(str$(time), 5); " seconds"
	print "Total energy: "; total_energy
	print
	if display_readouts = 1 then show_readouts()
	
	REM User control
	do_move()
	do_look()
	do_console()
	
	REM run particles
	run_particles()
	
   sync
loop

function show_readouts()
	for a = 0 to maxp
		if p(a).t <> 0 and a < 3
			ink rgb(p(a).cr, p(a).cg, p(a).cb), 0
			print "X Position of particle ";a;": "; p(a).x
			print "Y Position of particle ";a;": "; p(a).y
			print "Z Position of particle ";a;": "; p(a).z
			print
			print "X Velocity of particle ";a;": "; p(a).xv
			print "Y Velocity of particle ";a;": "; p(a).yv
			print "Z Velocity of particle ";a;": "; p(a).zv
			print
			print "Mass of particle ";a;": "; p(a).m
			print "Charge of particle ";a;": "; p(a).c
			print "Radius of particle ";a;": "; p(a).r
			print
		endif
	next a
endfunction

function apply_test_keys(a)
	REM Test
	if keystate(2)
		p(a).xv = p(a).xv * 1.01
		p(a).yv = p(a).yv * 1.01
		p(a).zv = p(a).zv * 1.01
	endif
	if keystate(3)
		p(a).xv = p(a).xv * 0.99
		p(a).yv = p(a).yv * 0.99
		p(a).zv = p(a).zv * 0.99
	endif          
	if keystate(4)
		dec p(a).yv, 0.01
	endif
	if keystate(5)
		inc p(a).m, 0.0001
	endif
	if keystate(6)
		p(a).xv = p(a).xv * 0.0
		p(a).yv = p(a).yv * 0.0
		p(a).zv = p(a).zv * 0.0
	endif
endfunction

function run_particles()
	inc time, 0.01
	
	for a = 0 to maxp
		if p(a).t <> 0
			REM apply forces and collision with border
			
			if p(a).m <> 0.0 then apply_gravity(a)
			if p(a).c <> 0.0 then apply_electromagnetism(a)
			
			apply_collision_walls(a)
			apply_test_keys(a)
			
		endif
	next a

	for a = 0 to maxp
		if p(a).t <> 0
			REM apply velocity and update 3d representation
			
			update_particle_positions(a)
			update_particle_objects(a)
			
		endif
	next a
	
	total_energy = 0.0
	for a = 0 to maxp
		if p(a).t <> 0
			REM calculate total energy of system as a means of error checking
			calculate_total_energy(a)
			inc total_energy, p(a).te
		endif
	next a
endfunction

function calculate_total_energy(a)
	` KE PE TE
	p(a).pe = 0.0
	p(a).ke = 0.0
	p(a).te = 0.0
	
	REM calculate KE (Kinetic Energy)
	REM KE = 1/2 * M * V^2
	p(a).ke = 0.5 * p(a).m * (sqrt(p(a).xv^2 + p(a).yv^2 + p(a).zv^2) ^ 2.0)
	
	REM Calculate PE (Potential energy)
	for b = 0 to maxp
		if p(b).t <> 0 and  b <> a `for each active particle besides a
			
			REM calculate distance
			dx = p(b).x - p(a).x
			dy = p(b).y - p(a).y
			dz = p(b).z - p(a).z
			Distance = sqrt(dx*dx + dy*dy + dz*dz)
			REM calculate potential for gravity
			REM V = G * M1 * M2 / D
			inc p(a).pe, (Gravity * p(a).m * p(b).m) /  Distance
			
		endif
	next b
	
	REM Calculate TE (Total energy)
	p(a).te = p(a).ke + p(a).pe
	
endfunction

function apply_gravity(a)
	local Force as double float
	local Distance as double float
	local dx as double float
	local dy as double float
	local dz as double float
	local Acceleration as double float
			
	for b = 0 to maxp
		if p(b).t <> 0 and p(b).m <> 0.0 and b <> a `for each active particle besides particle a with a mass
			
			REM calculate distance
			dx = p(b).x - p(a).x
			dy = p(b).y - p(a).y
			dz = p(b).z - p(a).z
			Distance = sqrt(dx*dx + dy*dy + dz*dz)
			REM calculate force
			`F = (G*m1*m2)/(DD)
			Force = (Gravity * p(a).m * p(b).m) / (Distance * Distance)
			REM calculate acceleration
			`A = F / M
			Acceleration = Force / p(a).m
			REM create unit vector from origin
			REM with magnitude equal to appropriate acceleratoin
			REM and direction pointing to particle B's relative position
			`we divide each value by D so that it is a unit vector
			`and multiply by acceleration so its magnitude equals acceleration
			inc p(a).xv, dx * (Acceleration/Distance)
			inc p(a).yv, dy * (Acceleration/Distance)
			inc p(a).zv, dz * (Acceleration/Distance)
			
		endif
	next b
			
endfunction

function apply_electromagnetism(a)
	local Force as double float
	local Distance as double float
	local dx as double float
	local dy as double float
	local dz as double float
	local Acceleration as double float
			
	for b = 0 to maxp
		if p(b).t <> 0 and p(b).c <> 0.0 and  b <> a `for each active particle besides particle a with a charge
			
			REM calculate distance
			dx = p(b).x - p(a).x
			dy = p(b).y - p(a).y
			dz = p(b).z - p(a).z
			Distance = sqrt(dx*dx + dy*dy + dz*dz)
			REM calculate force (electromagnetism)
			`F = (Kc*C1*C2)/(DD)
			Force = (Electrostatic * abs(p(a).c) * abs(p(b).c)) / (Distance * Distance)
			REM calculate acceleration
			`A = F / M
			Acceleration = Force` / p(a).m
			REM whether A accelerates away or towards B depends on likeness of charge
			if p(a).c > 0.0	`charge on A and B is same (positive)
				if p(b).c > 0.0 then Acceleration = -1.0 * Acceleration
			endif
			if p(a).c < 0.0	`charge on A and B is same (negative)
				if p(b).c < 0.0 then Acceleration = -1.0 * Acceleration
			endif
			REM create unit vector from origin
			REM with magnitude equal to appropriate acceleratoin
			REM and direction pointing to particle B's relative position
			`we divide each value by D so that it is a unit vector
			`and multiply by acceleration so its magnitude equals acceleration
			inc p(a).xv, dx * (Acceleration/Distance)
			inc p(a).yv, dy * (Acceleration/Distance)
			inc p(a).zv, dz * (Acceleration/Distance)
			
		endif
	next b

endfunction

function do_console()
	REM toggle sync rate
	if keystate(74)
		sync rate 5
	endif
	if keystate(78)
		sync rate 0
	endif
	if keystate(156)
		sync rate 50
	endif
	REM sc and oldsc and console are global
	oldsc = sc
	sc = scancode()
	REM Toggle console with tab key
	if sc = 15
		if oldsc = 0
			if console = 1
				console = 0
			else
				console = 1
				clear entry buffer
				console$(15) = ""
			endif
		endif
	endif
	REM Show console and type in console
	if console
		ink rgb(32,32,32),0
		box 0,0,sw,300
		set cursor 0,0
		for a=0 to 15
			ink rgb(64,255,64),0
			if left$(console$(a), 3) = ">> " then ink rgb(255,64,64),0
			print console$(a)
		next a
		
		olda$  = a$
		a$ = entry$()
		clear entry buffer
		if 9 = asc(a$) then a$="q"
		if 13 = asc(a$) then a$=""
		if 28 = oldsc and sc = 0 then do_console_command()
		if 8 = asc(a$)
		   a# = len(console$(15))
		   b$ = left$(console$(15), a# - 1)
		   console$(15) = b$
		else
		   console$(15) = console$(15)+a$
		endif
	endif
endfunction

function shift_console_items_up()
	for a = 0 to 14
		console$(a) = console$(a+1)
	next a
	console$(15) = ""
endfunction

function do_console_command()
	shift_console_items_up()
	c$ = console$(14)
	REM DO COMMANDS
	if c$ = "ini"
		reini()
	endif
	if c$ = "spd 0"
		sync rate 0
	endif
	if c$ = "spd 5"
		sync rate 5
	endif
	if c$ = "spd 50"
		sync rate 50
	endif
	if c$ = "read 0"
		display_readouts = 0
	endif
	if c$ = "read 1"
		display_readouts = 1
	endif
endfunction

INI:

REM Types
type particles_kinds
	name as string
	m as double float
	c as double float
	r as double float
	cr
	cg
	cb
	t
endtype
type particles
	x as double float
	y as double float
	z as double float
	xv as double float
	yv as double float
	zv as double float
	ix as double float
	iy as double float
	iz as double float
	ixv as double float
	iyv as double float
	izv as double float
	im as double float
	ic as double float
	ir as double float
	m as double float
	c as double float
	r as double float
	ke as double float
	pe as double float
	te as double float
	cr
	cg
	cb
	t
endtype

REM Display
sync on
sync rate 0
global sw	`screen width
global sh	`screen height
global display_readouts = 1
sw = screen width()
sh = screen height()

REM initiate console
global sc		`Scan code (for console)
global oldsc	`Previous scan code (for console)
global olda$	`previous entry (for console)
global console	`console toggle

DIM console$(16)
for a=0 to 15
	console$(a)=" "
next a

REM initiate globals
global time as double float
global total_energy as double float

REM initiate particles
DIM pk(32) as particles_kinds
DIM p(maxp+1) as particles

REM initiate camera
make camera 1
backdrop on 1
color backdrop 1, rgb(0,0,0)
Set Camera Range 1, 2000
position camera 1,0,50,50
point camera 1,0,0,0

REM initiate level
`load object "level/level.x",1
make object cube 2, border * 2.0
set object wireframe 2, 1
set object cull 2, 0
color object 2, rgb(0,255,0)
make light 1
set point light 1, 0, 0, 0
color light 1, 255, 255, 255
set light range 1, 100000

return

function update_particle_objects(a)
	position object a+100, p(a).x, p(a).y, p(a).z
	point object a+100, camera position x(1), camera position y(1), camera position z(1)
endfunction

function update_particle_positions(a)
	inc p(a).x, p(a).xv
	inc p(a).y, p(a).yv
	inc p(a).z, p(a).zv
endfunction

function reini()
	time = 0
	
	for a = 0 to maxp
		if object exist(a+100) then delete object a+100
		
		if p(a).t <> 0
			REM set initial properties
			p(a).x = p(a).ix
			p(a).y = p(a).iy
			p(a).z = p(a).iz
			p(a).xv = p(a).ixv
			p(a).yv = p(a).iyv
			p(a).zv = p(a).izv
			p(a).m = p(a).im
			p(a).c = p(a).ic
			p(a).r = p(a).ir
			REM make object
			make object sphere a+100, p(a).r, 4, 4
			position object a+100, p(a).ix, p(a).iy, p(a).iz
			color object a+100, rgb(p(a).cr, p(a).cg, p(a).cb)
		endif
	next a
endfunction

function apply_collision_walls(a)	
			
	if p(a).x > border or p(a).x < -1*border
		p(a).xv = p(a).xv * -1.0
	endif
	if p(a).y > border or p(a).y < -1*border
		p(a).yv = p(a).yv * -1.0
	endif
	if p(a).z > border or p(a).z < -1*border
		p(a).zv = p(a).zv * -1.0
	endif
			
endfunction

function do_move()
	local spd as float
	
	if console = 0
		REM var
		spd=0.5
		
		REM set speed and if hold LShift then go twice as fast
		if keystate(42)
			spd=spd*10
		endif
		
		REM move around
		if keystate(17)
			move camera 1, spd
		endif
		if keystate(31)
			move camera 1, -spd
		endif
		if keystate(30)
			position camera 1, newxvalue(camera position x(1), camera angle y(1)-90.0, spd), camera position y(1), newzvalue(camera position z(1), camera angle y(1)-90.0, spd)
		endif
		if keystate(32)
			position camera 1, newxvalue(camera position x(1), camera angle y(1)+90.0, spd), camera position y(1), newzvalue(camera position z(1), camera angle y(1)+90.0, spd)
		endif
		
		REM push upward or downward
		if keystate(200) or keystate(57)
			position camera 1, camera position x(1), camera position y(1) + spd, camera position z(1)
		endif
		if keystate(208) or keystate(56)
			position camera 1, camera position x(1), camera position y(1) - spd, camera position z(1)
		endif
	endif
endfunction

function do_look()
	local mx as float
	local my as float
	local mousedamp as float
	local cam as integer
	
	REM var
	mousedamp = 8.0
	
	REM look around, unless right clicking
	if mouseclick() > 1
		show mouse
	else
		hide mouse
		`still cursor
		position mouse 500,500
		`dampen amount to move
		mx=mousemovex()/mousedamp
		my=mousemovey()/mousedamp
		`turn camera object
		xrotate camera 1, wrapvalue(camera angle x(1)+my)
		yrotate camera 1, wrapvalue(camera angle y(1)+mx)
		`bound camera object angle
		if camera angle x(1)>89 and camera angle x(1)<180 then xrotate camera 1,89
		if camera angle x(1)<271 and camera angle x(1)>180 then xrotate camera 1,271
	endif
endfunction

press, on the numpad, "enter" for 60 FPS, or "-" for slow/10 FPS, or "+" for 0 FPS (fast as possible)

the beginning of the code sets all the particles and their settings (color, gravity)

it sets up a central planet, with a bunch of moons orbiting

the following code sets up numerous randomly placed planets (set the amount by setting the 'maxp' variable)

REM Particles lab
REM By qwe
REM -10.46MeV * ( e^(-r/b) / (r/b) ) where b = 1.43 fm

REM Constants
global maxp = 20
global border as double float				:	border = 1000.0
global Gravity as double float			:	Gravity = 1.0
global Electrostatic as double float	:	Electrostatic = 1.0

REM initialize
GOSUB INI

REM Test
remstart
for a = 0 to 16
	p(a).t = 1
	p(a).cr = rnd(255)
	p(a).cg = rnd(255)
	p(a).cb = rnd(255)
next a

p(0).ir = 10.0
p(0).im = 1.0

REM Inner

p(1).ir = 1.0
p(1).im = 0.000001
p(1).ix = 10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(1).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * 0.7071

p(2).ir = 1.0
p(2).im = 0.000001
p(2).ix = -10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(2).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * -0.7071

p(3).ir = 1.0
p(3).im = 0.000001
p(3).iz = 10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(3).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * 0.7071

p(4).ir = 1.0
p(4).im = 0.000001
p(4).iz = -10.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(4).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 10.0 )) * -0.7071

REM Outer

p(5).ir = 1.0
p(5).im = 0.000001
p(5).ix = 20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(5).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * 0.7071

p(6).ir = 1.0
p(6).im = 0.000001
p(6).ix = -20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(6).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * -0.7071

p(7).ir = 1.0
p(7).im = 0.000001
p(7).iz = 20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(7).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * 0.7071

p(8).ir = 1.0
p(8).im = 0.000001
p(8).iz = -20.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(8).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 20.0 )) * -0.7071

REM OuterOuter

p(9).ir = 1.0
p(9).im = 0.000001
p(9).ix = 30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(9).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * 0.7071

p(10).ir = 1.0
p(10).im = 0.000001
p(10).ix = -30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(10).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * -0.7071

p(11).ir = 1.0
p(11).im = 0.000001
p(11).iz = 30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(11).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * 0.7071

p(12).ir = 1.0
p(12).im = 0.000001
p(12).iz = -30.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(12).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 30.0 )) * -0.7071

REM OuterOuterOuter

p(13).ir = 1.0
p(13).im = 0.000001
p(13).ix = 40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(13).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * 0.7071

p(14).ir = 1.0
p(14).im = 0.000001
p(14).ix = -40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(14).iyv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * -0.7071

p(15).ir = 1.0
p(15).im = 0.000001
p(15).iz = 40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(15).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * 0.7071

p(16).ir = 1.0
p(16).im = 0.000001
p(16).iz = -40.0
`crit_velocity = sqrt( 2*G*m1 / dist )
p(16).ixv = (sqrt( (2.0 * Gravity * p(0).im) / 40.0 )) * -0.7071

remend

`remstart
for a = 0 to maxp
	p(a).t = 1
	p(a).cr = rnd(255)
	p(a).cg = rnd(255)
	p(a).cb = rnd(255)
	p(a).ir = 10.0
	p(a).im = 0.01
	p(a).ix = rnd(int(border)*2)-int(border)
	p(a).iy = rnd(int(border)*2)-int(border)
	p(a).iz = rnd(int(border)*2)-int(border)
	
	p(a).ixv = (rnd(2)-1) * ((rnd(100) + 0.0)/1000.0)
	p(a).iyv = (rnd(2)-1) * ((rnd(100) + 0.0)/1000.0)
	p(a).izv = (rnd(2)-1) * ((rnd(100) + 0.0)/1000.0)
next a
`remend
remstart
for a = 0 to 1
	p(a).t = 1
	p(a).cr = rnd(255)
	p(a).cg = rnd(255)
	p(a).cb = rnd(255)
	p(a).ir = 0.3333333
	p(a).im = 0.01
next a
p(0).ix = 5.0
remend