Sorry your browser is not supported!

`Raceing Game Contest Entry
`Very suckey too :P

SYNC ON : SYNC RATE 100 : SET DISPLAY MODE 1024,768,32 : HIDE MOUSE

SET TEXT SIZE 28 : SET TEXT FONT "Impact"
SYNC : TEXT 450,320,"Loading..." : SYNC

MAKE OBJECT BOX 1,50,50,75
COLOR OBJECT 1,RGB(255,255,0)
POSITION OBJECT 1,70,0,-500

MAKE OBJECT BOX 2,50,50,75
COLOR OBJECT 2,RGB(255,0,0)
POSITION OBJECT 2,-70,0,-500

MAKE OBJECT BOX 3,280,20,1150
COLOR OBJECT 3,RGB(100,20,0)
POSITION OBJECT 3,0,-25,350

DO
   SET TEXT SIZE 18 : TEXT 5,5,"Just hit ^ untill you reach the end of the cube :P"

   IF UPKEY()=1 THEN MOVE OBJECT 1,5
   POSITION CAMERA OBJECT POSITION X(1)-15,OBJECT POSITION Y(1)+120,OBJECT POSITION Z(1)-230

   MOVE OBJECT 2,3

   SYNC

   IF OBJECT POSITION Z(1)=>900 THEN GOSUB WINNER
   IF OBJECT POSITION Z(2)=>900 AND OBJECT POSITION Z(1) < 980 THEN GOSUB LOOSER

LOOP

WINNER:
   CLS : SET TEXT SIZE 18 : SET TEXT FONT "Verdana"
   CENTER TEXT 450,320,"You win and get lots of pie!! :D" : SYNC
   SLEEP 2000 : GOSUB ENDING

LOOSER:
   CLS : SET TEXT SIZE 18 : SET TEXT FONT "Verdana"
   CENTER TEXT 450,320,"You loose and get lots of pie!! :D" : SYNC
   SLEEP 2000 : GOSUB ENDING

ENDING:
   CLS : SET TEXT SIZE 28 : SET TEXT FONT "Impact"
   CENTER TEXT 450,320,"Deleting game objects..." : SYNC
   DELETE OBJECT 1 : DELETE OBJECT 2 : DELETE OBJECT 3
   END

`Standard setup
Sync on : Sync rate 60
Set display mode 1024,768,16 : Hide mouse
Autocam off : Set camera range 1,0x7fffffff

`Load Physics Functions
`#include "ZPhysics.dba"

`Start World Physics
Start_World_Physics()
Gravity(0)=9.8

`Start Object Physics
wheels=4
Start_Object_Physics(1)
For w=1 to wheels+1
Start_Object_Physics(w)
   MoveConstant(w)=200000
   TractionConstant(w)=0.992
   AngularTractionConstant(w)=0.98
   DragConstant(w)=0.47
   Mass(w)=10.0
   Weight(w)=Mass(w)*Gravity(0)
Next w

`Create Player
for w=1 to wheels
   Make object box w,10,50,50
   color object w,rgb(100,100,100)
next w
make object cube 5,50 : color object 5,rgb(255,0,0)

`Create world
make matrix 1,10000,10000,50,50
`randomize matrix 1,150
Update matrix 1

`**Main loop**
Do

`Controls
if upkey()=1 then set_linear_force(5,dirX(5)*10,dirY(5),DirZ(5)*10)
if leftkey()=1
set_angular_force(1,0,-(MoveSpeed(1)/25),0)
set_angular_force(2,0,-(MoveSpeed(2)/25),0)
endif
if rightkey()=1
set_angular_force(1,0,(MoveSpeed(1)/25),0)
set_angular_force(2,0,(MoveSpeed(2)/25),0)
endif
if upkey()=0 or downkey()=1 then set_linear_force(5,0,0,0)
if leftkey()=0 and rightkey()=0 then set_angular_force(1,0,0,0) : set_angular_force(2,0,0,0)
AngleY(5)=Curveangle(AngleY(1),AngleY(5),10.0)
For w=3 to 4
   AngleY(w)=AngleY(5)
Next w
for w=1 to wheels
   Xrotate object w,wrapvalue(AngleX(w)+(MoveSpeed(w)*(3.14*Object Size Z(w))))
next w

`Update Physics
for w=1 to wheels
   VelX(w)=VelX(5) : VelY(w)=VelY(5) : VelZ(w)=VelZ(5)
   PositionX(w)=PositionX(5) : PositionY(w)=PositionY(5) : PositionZ(w)=PositionZ(5)
   PositionY(w)=get ground height(1,object position x(w),object position z(w))+Object Size Y(w)/2
next w
PositionY(5)=get ground height(1,positionx(5),positionz(5))+(Object Size Y(5)/2)+(Object Size Y(5)/5)
Update_Physics(1,315,60)
Update_Physics(2,45,60)
Update_Physics(3,135,60)
Update_Physics(4,215,60)
Update_Physics(5,180,10)

`Camera
lx#=positionX(5)
lz#=positionZ(5)
ca#=curveangle(angley(1),ca#,20.0)
cx#=lx#-sin(ca#)*600
cz#=lz#-cos(ca#)*600
cy#=positiony(5)+300
position camera cx#,cy#,cz#
point camera lx#,positiony(5),lz#

`Info
text 1020-text width(str$(screen fps())),740,str$(screen fps())

`**End loop**
Sync
Loop

`**Functions**
Function Start_World_Physics()
Dim Gravity(0) as float
Dim Windx(0) as float
Dim Windy(0) as float
Dim Windz(0) as float
Endfunction

Function Start_Object_Physics(o)
Dim Positionx(o) as Float : Dim Positiony(o) as Float : Dim Positionz(o) as Float
Dim OldPosX(o) as Float : Dim OldPosY(o) as Float : Dim OldPosZ(o) as Float
Dim Anglex(o) as Float : Dim Angley(o) as Float : Dim Anglez(o) as Float
Dim LongForceX(o) as Float : Dim LongForceY(o) as Float : Dim LongForceZ(o) as Float
Dim TracForceX(o) as Float : Dim TracForceY(o) as Float : Dim TracForceZ(o) as Float
Dim DragForceX(o) as Float : Dim DragForceY(o) as Float : Dim DragForceZ(o) as Float
Dim AngForceX(o) as Float : Dim AngForceY(o) as Float : Dim AngForceZ(o) as Float
Dim AngVelX(o) as Float : Dim AngVelY(o) as Float : Dim AngVelZ(o) as Float
Dim AngAccX(o) as Float : Dim AngAccY(o) as Float : Dim AngAccZ(o) as Float
Dim VelX(o) as Float : Dim VelY(o) as Float : Dim VelZ(o) as Float
Dim AccX(o) as Float : Dim AccY(o) as Float : Dim AccZ(o) as Float
Dim DirX(o) as Float : Dim DirY(o) as Float : Dim DirZ(o) as Float
Dim MoveTime(o) as Float : Dim MoveDist(o) as Float : Dim Power(o) as Float
Dim Movespeed(o) as Float : Dim MoveConstant(o) as Float : Dim Acceleration(o) as Float
Dim TractionConstant(o) as Float : Dim AngularTractionConstant(o) as Float
Dim DragConstant(o) as Float
Dim Mass(o) as Float : Dim Weight(o) as Float
Dim KineticEnergy(o) as Float: Dim WorkDone(o) as Float
Dim LatForce(o) as Float : Dim PushX(o) as Float : Dim PushY(o) as Float : Dim PushZ(o) as Float
Dim Force(o) as Float
Endfunction

Function Set_Linear_Force(o,x,y,z)
LongForceX(o)=x : LongForceY(o)=y : LongForceZ(o)=z
Endfunction

Function Set_Linear_Velocity(o,x,y,z)
VelX(o)=x : VelY(o)=y : VelZ(o)=z
Endfunction

Function Set_Angular_Force(o,x,y,z)
AngForceX(o)=x : AngForceY(o)=y : AngForceZ(o)=z
Endfunction

Function Set_Angular_Velocity(o,x,y,z)
AngVelX(o)=x : AngVelY(o)=y : AngVelZ(o)=z
Endfunction

Function Move_Object(o,spd)
VelX(o)=dirX(o)*spd
VelY(o)=dirY(o)*spd
VelZ(o)=dirZ(o)*spd
Endfunction

Function Push_Object(o,spd)
LongForceX(o)=dirX(o)*spd
LongForceY(o)=dirY(o)*spd
LongForceZ(o)=dirZ(o)*spd
Endfunction

Function Control_Object_Using_Forces(o,push,turn)
if upkey()=1 then Push_Object(o,push)
if downkey()=1 then Push_Object(o,-push)
if leftkey()=1 then set_angular_Force(o,0,-(MoveSpeed(o)/turn),0)
if rightkey()=1 then set_angular_Force(o,0,(MoveSpeed(o)/turn),0)
if upkey()=0 and downkey()=0 then Push_Object(o,0)
if leftkey()=0 and rightkey()=0 then set_angular_force(o,0,0,0)
Endfunction

Function Control_Object_Using_Velocities(o,spd,turn)
if upkey()=1 then Move_Object(o,spd)
if downkey()=1 then Move_Object(o,-spd)
if leftkey()=1 then set_angular_velocity(o,0,-turn,0)
if rightkey()=1 then set_angular_velocity(o,0,turn,0)
if upkey()=0 and downkey()=0 then Move_Object(o,0)
if leftkey()=0 and rightkey()=0 then set_angular_velocity(o,0,0,0)
Endfunction

Function Update_Physics(o,say#,out#)
OldPosX(o)=PositionX(o)-VelX(o) : OldPosY(o)=PositionY(o)-VelY(o) : OldPosZ(o)=PositionZ(o)-VelZ(o)
DirX(o)=sin(AngleY(o))*cos(AngleX(o)) : DirY(o)=-sin(AngleX(o)) : DirZ(o)=cos(AngleY(o))*cos(AngleX(o))
TracForceX(o)=DirX(o)*MoveConstant(o) : TracForceY(o)=DirY(o)*MoveConstant(o) : TracForceZ(o)=DirZ(o)*MoveConstant(o)
AccX(o)=LongForceX(o)/Mass(o) : AccY(o)=LongForceY(o)/Mass(o) : AccZ(o)=LongForceZ(o)/Mass(o)
VelX(o)=VelX(o)+AccX(o) : VelX(o)=VelX(o)*TractionConstant(o) : VelY(o)=VelY(o)+AccY(o) : VelY(o)=VelY(o)*TractionConstant(o) : VelZ(o)=VelZ(o)+AccZ(o) : VelZ(o)=VelZ(o)*TractionConstant(o)
MoveSpeed(o)=sqrt(VelX(o)^2 + VelY(o)^2 + VelZ(o)^2)
Acceleration(o)=Sqrt(AccX(o)^2+AccY(o)^2+AccZ(o)^2)
Force(o)=Sqrt(LongForceX(o)^2+LongForceY(o)^2+LongForceZ(o)^2)
DragForceX(o)=-DragConstant(o)*VelX(o)*MoveSpeed(o) : DragForceY(o)=-DragConstant(o)*VelY(o)*MoveSpeed(o) : DragForceZ(o)=-DragConstant(o)*VelZ(o)*MoveSpeed(o)
AngAccX(o)=AngForceX(o)/Mass(o) : AngAccY(o)=AngForceY(o)/Mass(o) : AngAccZ(o)=AngForceZ(o)/Mass(o)
AngVelX(o)=AngVelX(o)+AngAccX(o) : AngVelY(o)=AngVelY(o)+AngAccY(o) : AngVelZ(o)=AngVelZ(o)+AngAccZ(o)
AngVelX(o)=AngVelX(o)*AngularTractionConstant(o) : AngVelY(o)=AngVelY(o)*AngularTractionConstant(o) : AngVelZ(o)=AngVelZ(o)*AngularTractionConstant(o)
PositionX(o)=PositionX(o)+(DirX(o)*Movespeed(o)) : PositionY(o)=PositionY(o)+(DirY(o)*Movespeed(o)) : PositionZ(o)=PositionZ(o)+(DirZ(o)*Movespeed(o))
AngleX(o)=Wrapvalue(AngleX(o)+AngVelX(o)) : AngleY(o)=Wrapvalue(AngleY(o)+AngVelY(o)) : AngleZ(o)=Wrapvalue(AngleZ(o)+AngVelZ(o))
KineticEnergy(o)=(Mass(o)/2)*(Movespeed(o)^2)
if MoveSpeed(o)<0.9 then MoveSpeed(o)=MoveSpeed(o)*0.0
If MoveSpeed(o)<>0 then MoveTime(o)=Movetime(o)+1.0 else MoveTime(o)=0.0
If MoveTime(o)>0.0 then Power(o)=KineticEnergy(o)/MoveTime(o) else Power(o)=0.0
MoveDist(o)=MoveSpeed(o)*MoveTime(o) : WorkDone(o)=Force(o)*MoveDist(o)
LatForce(o)=abs(AngVelY(o))*(MoveSpeed(o)/Mass(o))
If AngVelY(o)>0.0
   PushX(o)=(sin(Wrapvalue(AngleY(o)-90))*LatForce(o))
   PushZ(o)=(cos(Wrapvalue(AngleY(o)-90))*LatForce(o))
endif
If AngVelY(o)<0.0
   PushX(o)=(sin(Wrapvalue(AngleY(o)+90))*LatForce(o))
   PushZ(o)=(cos(Wrapvalue(AngleY(o)+90))*LatForce(o))
endif
Position object o,Positionx(o)+PushX(o)+(Sin(wrapvalue(AngleY(o)+say#))*out#),Positiony(o)+PushY(o),Positionz(o)+PushZ(o)+(Cos(wrapvalue(AngleY(o)+say#))*out#)
Rotate object o,Anglex(o),Angley(o),Anglez(o)
Endfunction

SYNC ON : SYNC RATE 60 : HIDE MOUSE
SET DISPLAY MODE 1024,768,32

SET TEXT SIZE 28 : SET TEXT FONT "Impact"
SYNC : CENTER TEXT 512,350,"Loading..." : SYNC

INK RGB(200,200,200),0
BOX 1,1,128,128
FOR I=0 TO 2000
   COLR=RND(255)
   INK RGB(COLR,COLR,COLR),0
   DOT RND(128),RND(128)
NEXT I
GET IMAGE 1,1,1,128,128

CREATE BITMAP 1,128,128
SET CURRENT BITMAP 1
FOR I=1 TO 16
   FOR J=1 TO 16
      IF RND(1)=0
         INK RGB(255,0,0),0
      ELSE
         INK RGB(255,255,0),0
      ENDIF
      BOX ((J*8)-8),((I*8)-8),(J*8),(I*8)
   NEXT J
NEXT I
BLUR BITMAP 1,6
GET IMAGE 2,1,1,128,128
SET CURRENT BITMAP 0
DELETE BITMAP 1

MAKE MATRIX 1,5000,5000,5,5
PREPARE MATRIX TEXTURE 1,1,1,1
SET MATRIX TEXTURE 1,1,1
FOR I=1 TO 4
   SET MATRIX HEIGHT 1,I,1,1000
   SET MATRIX HEIGHT 1,I,4,1000
NEXT I
FOR I=2 TO 3
   SET MATRIX HEIGHT 1,1,I,1000
   SET MATRIX HEIGHT 1,2,I,1000
   SET MATRIX HEIGHT 1,3,I,1000
   SET MATRIX HEIGHT 1,4,I,1000
NEXT I
SET MATRIX HEIGHT 1,2,2,0
SET MATRIX HEIGHT 1,3,3,0
SET MATRIX HEIGHT 1,1,4,0
SET MATRIX HEIGHT 1,4,1,0
UPDATE MATRIX 1

LAVA=1
MAKE OBJECT PLAIN 1,7500,7500
TEXTURE OBJECT 1,2
XROTATE OBJECT 1,-90
SET OBJECT COLLISION ON 1
SET OBJECT LAVA,1,0,0,1,0,0,0
POSITION OBJECT LAVA,2500,750,2500

PLR=2
MAKE OBJECT SPHERE PLR,50,12,12
SET OBJECT COLLISION ON PLR

FOG ON : FOG COLOR RGB(255,60,0) : FOG DISTANCE 1000
BACKDROP ON : COLOR BACKDROP RGB(255,60,0)

START:
PLRA=90 : PLRY#=GET GROUND HEIGHT(1,1000,1000)+25
COLOR OBJECT PLR,RGBB(200)
POSITION OBJECT PLR,1000,PLRY#,1000
YROTATE OBJECT PLR,PLRA

DO
   MOVE OBJECT PLR,20
   IF LEFTKEY()=1 THEN PLRA=PLRA-2
   IF RIGHTKEY()=1 THEN PLRA=PLRA+2
   YROTATE OBJECT PLR,PLRA
   PLRY#=GET GROUND HEIGHT(1,OBJECT POSITION X(PLR),OBJECT POSITION Z(PLR))+25
   POSITION OBJECT PLR,OBJECT POSITION X(PLR),PLRY#,OBJECT POSITION Z(PLR)
   SET CAMERA TO OBJECT ORIENTATION PLR
   POSITION CAMERA OBJECT POSITION X(PLR),PLRY#+100,OBJECT POSITION Z(PLR)
   MOVE CAMERA -175 : POINT CAMERA OBJECT POSITION X(PLR),PLRY#+15,OBJECT POSITION Z(PLR)

   IF OBJECT COLLISION (PLR,LAVA) THEN GOSUB DEATH

   SYNC
LOOP

DEATH:
   INK RGB(255,255,255),0
   COLOR OBJECT PLR,0 : SYNC
   CENTER TEXT 512,350,"You Died" : SYNC
   SLEEP 1000 : CLS
   GOSUB START

`Standard setup
Sync on : Sync rate 60
Set display mode 1024,768,16 : Hide mouse
Autocam off : Set camera range 1,0x7fffffff

`Load Physics Functions
`#include "ZPhysics.dba"

`Start World Physics
Start_World_Physics()
Set_Standard_Gravity()
Set_Standard_Drag()

`Start Object Physics
wheels=4
Start_Object_Physics(1)
For w=1 to wheels+1
Start_Object_Physics(w)
   MoveConstant(w)=150000
   TractionConstant(w)=0.99
   AngularTractionConstant(w)=0.98
   Mass(w)=10.0
   Weigh_Object(w)
Next w
turn#=25.0

`Create Player
create bitmap 1,50,50 : set current bitmap 1
cls rgb(0,0,0)
ink rgb(0,0,0),0
box 0,0,50,50
ink rgb(100,100,100),0
for x=1 to 50
   for y=1 to 50
      if sqrt((x-50)^2+(y-50)^2)<50
         box x,y,x+1,y+1
      endif
   next x
next y
get image 1,0,0,50,50
set current bitmap 0 : delete bitmap 1
for w=1 to wheels
   Make object plain w,50,50 : yrotate object w,90 : fix object pivot w
   texture object w,1 : set object transparency w,1
next w
make object box 5,50,30,80 : color object 5,rgb(255,0,0)

`Create world --

`Textures-

`Grass
create bitmap 1,128,128 : set current bitmap 1
cls rgb(0,128,0)
for d=1 To 10000
ink rgb(0,rnd(255),rnd(50)),0
dot rnd(128),rnd(128)
next d
blur bitmap 1,4
get image 10,0,0,128,128
set current bitmap 0 : delete bitmap 1
create bitmap 1,128,128 : set current bitmap 1
cls rgb(128,128,128)
for d=1 To 10000
ink rgb(rnd(200),rnd(200),rnd(200)),0
dot rnd(128),rnd(128)
next d
blur bitmap 1,4
get image 11,0,0,128,128
set current bitmap 0 : delete bitmap 1
create bitmap 1,512,128 : set current bitmap 1
paste image 10,0,0
paste image 11,128,0
get image 15,0,0,512,128
set current bitmap 0 : delete bitmap 1

`Create and Mold
make matrix 1,20000,20000,50,50
wallh=1000
for x=1 to 49
   for z=1 to 49
      set matrix height 1,x,1,wallh : set matrix height 1,x,50,wallh
      set matrix height 1,1,z,wallh : set matrix height 1,50,z,wallh
   next x
next z
for x=2 to 49
   for z=2 to 49
      set matrix height 1,x,2,wallh : set matrix height 1,x,49,wallh
      set matrix height 1,2,z,wallh : set matrix height 1,49,z,wallh
   next x
next z
for x=10 to 40
   for z=10 to 40
      set matrix height 1,x,z,wallh
   next x
next z
for x=20 to 30
   for z=30 to 40
      set matrix height 1,x,z,0
   next x
next z
for z=40 to 49
   set matrix height 1,25,z,wallh
next z
for x=2 to 9
   set matrix height 1,x,21,30
   set matrix height 1,x,22,60
   set matrix height 1,x,23,150
   set matrix height 1,x,24,300
   set matrix height 1,x,25,700
next x

`Texture
prepare matrix texture 1,15,4,1
for x=1 to 49
   for z=1 to 49
      set matrix tile 1,x,z,2
   next x
next z
for x=3 to 8
   for z=3 to 47
      set matrix tile 1,x,z,1
   next x
next z
for x=9 to 23
   for z=41 to 47
      set matrix tile 1,x,z,1
   next x
next z
for x=20 to 29
   for z=31 to 38
      set matrix tile 1,x,z,1
   next x
next z
for x=26 to 47
   for z=41 to 47
      set matrix tile 1,x,z,1
   next x
next z
for x=41 to 47
   for z=3 to 47
      set matrix tile 1,x,z,1
   next x
next z
for x=3 to 47
   for z=3 to 8
      set matrix tile 1,x,z,1
   next x
next z
for x=20 to 23
   for z=31 to 47
      set matrix tile 1,x,z,1
   next x
next z
for x=26 to 29
   for z=31 to 47
      set matrix tile 1,x,z,1
   next x
next z

`Matrix Lightng
for z=1 to 49
for x=1 to 49
h8#=get matrix height(1,x,z-1)
h4#=get matrix height(1,x-1,z)
h#=get matrix height(1,x,z)
h2#=get matrix height(1,x,z)
x1#=(x-1)*25.0 : y1#=h#
x2#=(x+0)*25.0 : y2#=h4#
dx#=x2#-x1#
dy#=y2#-y1#
ax#=atanfull(dx#,dy#)
ax#=wrapvalue(90-ax#)
z1#=(z-1)*25.0 : y1#=h2#
z2#=(z+0)*25.0 : y2#=h8#
dz#=z2#-z1#
dy#=y2#-y1#
az#=atanfull(dz#,dy#)
az#=wrapvalue(90-az#)
nx#=sin(wrapvalue(ax#+45))
ny#=cos(ax#)
nz#=sin(wrapvalue(az#+45))
set matrix normal 1,x,z,nx#,ny#,nz#
next x
next z
update matrix 1

`Lighting
fogdist=30000
amb=30
lightrange=8000
lightheight=200
fogdist=30000
amb=30
lightrange=8000
lightheight=200
fog on
backdrop on
fog distance fogdist
fog color backcl
color backdrop backcl
set ambient light amb
color ambient light backcl
for l=1 to 5
   make light l
   color light l,rgb(255,255,255)
   set light range l,lightrange
next l
   position light 1,2400,lightheight,19000
   position light 2,8720,lightheight,12300
   position light 3,11950,lightheight,12050
   position light 4,19030,lightheight,19150
   position light 5,19050,lightheight,1600

`Create Light posts
for l=1 to 5
   make object box 100+l,10,lightheight*2,10
   make object cone 105+l,100
next l
   position object 101,2400,lightheight,19000
   position object 102,8720,lightheight,12300
   position object 103,11950,lightheight,12050
   position object 104,19030,lightheight,19150
   position object 105,19050,lightheight,1600
      position object 106,2400,lightheight*2,19000
      position object 107,8720,lightheight*2,12300
      position object 108,11950,lightheight*2,12050
      position object 109,19030,lightheight*2,19150
      position object 110,19050,lightheight*2,1600




`**Main loop**
Do

`Controls
if upkey()=1 or keystate(17)=1 then set_linear_force(5,dirX(5)*10,dirY(5),DirZ(5)*10)
if leftkey()=1
   set_angular_force(1,0,-(Movespeed(1)/turn#),0)
   set_angular_force(2,0,-(Movespeed(1)/turn#),0)
endif
if rightkey()=1
   set_angular_force(1,0,(Movespeed(1)/turn#),0)
   set_angular_force(2,0,(Movespeed(1)/turn#),0)
endif

if upkey()=0 and keystate(17)=0 then set_linear_force(5,0,0,0)
if leftkey()=0 and rightkey()=0 then set_angular_force(1,0,0,0) : set_angular_force(2,0,0,0)
Set_Next_Angle_Y(5,1,10.0)
For w=3 to 4
   Set_Next_Angle_Y(w,5,3.0)
Next w

`Update Physics
for w=1 to wheels
   Set_Linear_Velocity(w,VelX(5),VelY(5),VelZ(5))
   Position_At_Object(w,5)
   PositionY(w)=get ground height(1,object position x(w),object position z(w))+Object Size Y(w)/2
next w
PositionY(5)=get ground height(1,positionx(5),positionz(5))+(Object Size Y(5)/2)+(Object Size Y(5)/5)
Calculate_Offsets(1,0,315,60) : Calculate_Offsets(2,0,45,60)
Calculate_Offsets(3,0,135,60) : Calculate_Offsets(4,0,215,60)
for w=1 to 5
   Standard_Update(w)
next w

`Camera
lx#=positionX(5)
lz#=positionZ(5)
ca#=curveangle(angley(1),ca#,20.0)
cx#=lx#-sin(ca#)*600
cz#=lz#-cos(ca#)*600
cy#=positiony(5)+300
position camera cx#,cy#,cz#
point camera lx#,positiony(5),lz#

`Info
text 1020-text width(str$(screen fps())),740,str$(screen fps())

`**End loop**
Sync
Loop

`**Functions**
Function Start_World_Physics()
Dim Gravity(0) as float
Dim Drag(0) as Float
Dim Windx(0) as float
Dim Windy(0) as float
Dim Windz(0) as float
Endfunction

Function Start_Object_Physics(o)
Dim Positionx(o) as Float : Dim Positiony(o) as Float : Dim Positionz(o) as Float
Dim OldPosX(o) as Float : Dim OldPosY(o) as Float : Dim OldPosZ(o) as Float
Dim Anglex(o) as Float : Dim Angley(o) as Float : Dim Anglez(o) as Float
Dim LongForceX(o) as Float : Dim LongForceY(o) as Float : Dim LongForceZ(o) as Float
Dim TracForceX(o) as Float : Dim TracForceY(o) as Float : Dim TracForceZ(o) as Float
Dim DragForceX(o) as Float : Dim DragForceY(o) as Float : Dim DragForceZ(o) as Float
Dim AngForceX(o) as Float : Dim AngForceY(o) as Float : Dim AngForceZ(o) as Float
Dim AngVelX(o) as Float : Dim AngVelY(o) as Float : Dim AngVelZ(o) as Float
Dim AngAccX(o) as Float : Dim AngAccY(o) as Float : Dim AngAccZ(o) as Float
Dim VelX(o) as Float : Dim VelY(o) as Float : Dim VelZ(o) as Float
Dim AccX(o) as Float : Dim AccY(o) as Float : Dim AccZ(o) as Float
Dim DirX(o) as Float : Dim DirY(o) as Float : Dim DirZ(o) as Float
Dim OffsetX(o) as Float : Dim OffsetY(o) as Float : Dim OffsetZ(o) as Float
Dim MoveTime(o) as Float : Dim MoveDist(o) as Float : Dim Power(o) as Float
Dim Movespeed(o) as Float : Dim MoveConstant(o) as Float : Dim Acceleration(o) as Float
Dim TractionConstant(o) as Float : Dim AngularTractionConstant(o) as Float
Dim Mass(o) as Float : Dim Weight(o) as Float
Dim KineticEnergy(o) as Float: Dim WorkDone(o) as Float
Dim LatForce(o) as Float : Dim PushX(o) as Float : Dim PushY(o) as Float : Dim PushZ(o) as Float
Dim Force(o) as Float
Endfunction

Function Set_Standard_Gravity()
Gravity(0)=9.8
Endfunction

Function Set_Standard_Drag()
Drag(0)=0.47
Endfunction

Function Weigh_Object(o)
Weight(o)=Mass(o)*Gravity(0)
Endfunction

Function Set_Linear_Force(o,x,y,z)
LongForceX(o)=x : LongForceY(o)=y : LongForceZ(o)=z
Endfunction

Function Set_Linear_Velocity(o,x,y,z)
VelX(o)=x : VelY(o)=y : VelZ(o)=z
Endfunction

Function Set_Angular_Force(o,x,y,z)
AngForceX(o)=x : AngForceY(o)=y : AngForceZ(o)=z
Endfunction

Function Set_Angular_Velocity(o,x,y,z)
AngVelX(o)=x : AngVelY(o)=y : AngVelZ(o)=z
Endfunction

Function Move_Object(o,spd)
VelX(o)=dirX(o)*spd
VelY(o)=dirY(o)*spd
VelZ(o)=dirZ(o)*spd
Endfunction

Function Push_Object(o,spd)
LongForceX(o)=dirX(o)*spd
LongForceY(o)=dirY(o)*spd
LongForceZ(o)=dirZ(o)*spd
Endfunction

Function Control_Object_Using_Forces(o,push,turn)
if upkey()=1 then Push_Object(o,push)
if downkey()=1 then Push_Object(o,-push)
if leftkey()=1 then set_angular_Force(o,0,-(MoveSpeed(o)/turn),0)
if rightkey()=1 then set_angular_Force(o,0,(MoveSpeed(o)/turn),0)
if upkey()=0 and downkey()=0 then Push_Object(o,0)
if leftkey()=0 and rightkey()=0 then set_angular_force(o,0,0,0)
Standard_Update(o)
Endfunction

Function Control_Object_Using_Velocities(o,spd,turn)
if upkey()=1 then Move_Object(o,spd)
if downkey()=1 then Move_Object(o,-spd)
if leftkey()=1 then set_angular_velocity(o,0,-turn,0)
if rightkey()=1 then set_angular_velocity(o,0,turn,0)
if upkey()=0 and downkey()=0 then Move_Object(o,0)
if leftkey()=0 and rightkey()=0 then set_angular_velocity(o,0,0,0)
Standard_Update(o)
Endfunction

Function Set_Next_Angle_X(a,b,smth#)
AngleX(a)=Curveangle(AngleX(b),AngleX(a),smth#)
Endfunction

Function Set_Next_Angle_Y(a,b,smth#)
AngleY(a)=Curveangle(AngleY(b),AngleY(a),smth#)
Endfunction

Function Set_Next_Angle_Z(a,b,smth#)
AngleZ(a)=Curveangle(AngleZ(b),AngleZ(a),smth#)
Endfunction

Function Position_At_Object(a,b)
PositionX(a)=PositionX(b)
PositionY(a)=PositionY(b)
PositionZ(a)=PositionZ(b)
Endfunction

Function Get_Old_Position(o)
OldPosX(o)=PositionX(o)-VelX(o)
OldPosY(o)=PositionY(o)-VelY(o)
OldPosZ(o)=PositionZ(o)-VelZ(o)
Endfunction

Function Get_Object_Direction(o)
DirX(o)=sin(AngleY(o))*cos(AngleX(o))
DirY(o)=-sin(AngleX(o))
DirZ(o)=cos(AngleY(o))*cos(AngleX(o))
Endfunction

Function Get_Traction_Force(o)
TracForceX(o)=DirX(o)*MoveConstant(o)
TracForceY(o)=DirY(o)*MoveConstant(o)
TracForceZ(o)=DirZ(o)*MoveConstant(o)
Endfunction

Function Get_Drag_Force(o)
DragForceX(o)=-Drag(0)*VelX(o)*MoveSpeed(o)
DragForceY(o)=-Drag(0)*VelY(o)*MoveSpeed(o)
DragForceZ(o)=-Drag(0)*VelZ(o)*MoveSpeed(o)
Endfunction

Function Get_Force(o)
Force(o)=Sqrt(LongForceX(o)^2+LongForceY(o)^2+LongForceZ(o)^2)
Endfunction

Function Get_Speed(o)
MoveSpeed(o)=sqrt(VelX(o)^2 + VelY(o)^2 + VelZ(o)^2)
Endfunction

Function Get_Acceleration(o)
Acceleration(o)=Sqrt(AccX(o)^2+AccY(o)^2+AccZ(o)^2)
Endfunction

Function Set_Acceleration(o)
AccX(o)=LongForceX(o)/Mass(o)
AccY(o)=LongForceY(o)/Mass(o)
AccZ(o)=LongForceZ(o)/Mass(o)
Endfunction

Function Set_Velocity(o)
VelX(o)=VelX(o)+AccX(o) : VelX(o)=VelX(o)*TractionConstant(o)
VelY(o)=VelY(o)+AccY(o) : VelY(o)=VelY(o)*TractionConstant(o)
VelZ(o)=VelZ(o)+AccZ(o) : VelZ(o)=VelZ(o)*TractionConstant(o)
Endfunction

Function Get_Angular_Acceleration(o)
AngAccX(o)=AngForceX(o)/Mass(o)
AngAccY(o)=AngForceY(o)/Mass(o)
AngAccZ(o)=AngForceZ(o)/Mass(o)
Endfunction

Function Get_Angular_Velocity(o)
AngVelX(o)=AngVelX(o)+AngAccX(o) : AngVelX(o)=AngVelX(o)*AngularTractionConstant(o)
AngVelY(o)=AngVelY(o)+AngAccY(o) : AngVelY(o)=AngVelY(o)*AngularTractionConstant(o)
AngVelZ(o)=AngVelZ(o)+AngAccZ(o) : AngVelZ(o)=AngVelZ(o)*AngularTractionConstant(o)
Endfunction

Function Update_Coordinates_Straight(o)
PositionX(o)=PositionX(o)+(DirX(o)*Movespeed(o))
PositionY(o)=PositionY(o)+(DirY(o)*Movespeed(o))
PositionZ(o)=PositionZ(o)+(DirZ(o)*Movespeed(o))
Endfunction

Function Update_Coordinates_Slide(o)
PositionX(o)=PositionX(o)+VelX(o)
PositionY(o)=PositionY(o)+VelY(o)
PositionZ(o)=PositionZ(o)+VelZ(o)
Endfunction

Function Update_Angles(o)
AngleX(o)=Wrapvalue(AngleX(o)+AngVelX(o))
AngleY(o)=Wrapvalue(AngleY(o)+AngVelY(o))
AngleZ(o)=Wrapvalue(AngleZ(o)+AngVelZ(o))
Endfunction

Function Get_Lateral_Forces(o)
LatForce(o)=abs(AngVelY(o))*(MoveSpeed(o)/Mass(o))
If AngVelY(o)>0.0
   PushX(o)=(sin(Wrapvalue(AngleY(o)-90))*LatForce(o))
   PushZ(o)=(cos(Wrapvalue(AngleY(o)-90))*LatForce(o))
endif
If AngVelY(o)<0.0
   PushX(o)=(sin(Wrapvalue(AngleY(o)+90))*LatForce(o))
   PushZ(o)=(cos(Wrapvalue(AngleY(o)+90))*LatForce(o))
endif
Endfunction

Function Calculate_Offsets(o,sax#,say#,out#)
OffsetX(o)=(Sin(Wrapvalue(AngleY(o)+say#))*Cos(Wrapvalue(AngleX(o)+sax#)))*Out#
OffsetY(o)=-Sin(Wrapvalue(AngleX(o)+sax#))*600
OffsetZ(o)=(Cos(Wrapvalue(AngleY(o)+say#))*Cos(Wrapvalue(AngleX(o)+sax#)))*out#
Endfunction

Function Update_Object_Position(o)
Position Object o,PositionX(o),PositionY(o),PositionZ(o)
Endfunction

Function Update_Position_With_Push_Forces(o,a,b,c)
Position object o,Positionx(o)+(PushX(o)*a),Positiony(o)+(PushY(o)*b),Positionz(o)+(PushZ(o)*c)
Endfunction

Function Update_Position_With_Offsets(o,a,b,c)
Position Object o,Positionx(o)+(offsetx(o)*a),Positiony(o)+(offsety(o)*b),PositionZ(o)+(offsetZ(o)*c)
Endfunction

Function Fully_Update_object_Position(o)
Position Object o,PositionX(o)+PushX(o)+OffsetX(o),PositionY(o)+PushY(o)+OffsetY(o),PositionZ(o)+PushZ(o)+OffsetZ(o)
Endfunction

Function Update_object_Angles(o)
Rotate object o,AngleX(o),AngleY(o),AngleZ(o)
Endfunction

Function Standard_Update(o)
Get_Old_Position(o)
Get_Object_Direction(o)
Get_Traction_Force(o)
Get_Drag_Force(o)
Get_Force(o)
Get_Speed(o)
Get_Acceleration(o)
Set_Acceleration(o)
Set_Velocity(o)
Get_Angular_Acceleration(o)
Get_Angular_Velocity(o)
Update_Coordinates_Straight(o)
Update_Angles(o)
Get_Lateral_Forces(o)
Fully_Update_Object_Position(o)
Update_Object_Angles(o)
Endfunction

`Standard setup
Sync on : Sync rate 60
Set display mode 1024,768,16 : Hide mouse
Autocam off : Set camera range 1,0x7fffffff

`Load Physics Functions
`#include "ZPhysics.dba"

`Start World Physics
Start_World_Physics()
Set_Standard_Gravity()
Set_Standard_Drag()

`Start Object Physics
wheels=4
Start_Object_Physics(1)
For w=1 to wheels+1
Start_Object_Physics(w)
   Set_Move_Constant(w,250000)
   Set_Traction_Constant(w,0.98)
   Set_Angular_Traction_Constant(w,0.98)
   Set_Body_Mass(w,10.0)
   Weigh_Object(w)
Next w
turn#=25.0

`Create Player
for w=1 to wheels
   make object box w,10,30,30
next w
make object box 5,50,30,80 : color object 5,rgb(255,0,0)

`Create world
create bitmap 1,128,128 : set current bitmap 1
cls rgb(0,128,0)
for d=1 To 10000
ink rgb(0,rnd(255),rnd(50)),0
dot rnd(128),rnd(128)
next d
blur bitmap 1,4
get image 10,0,0,128,128
set current bitmap 0 : delete bitmap 1

`Create and Mold
make matrix 1,20000,20000,50,50
randomize matrix 1,50
prepare matrix texture 1,10,1,1
update matrix 1

`**Main loop**
Do

`Controls
if upkey()=1 or keystate(17)=1 then set_linear_force(5,dirX(5)*10,dirY(5),DirZ(5)*10)
if leftkey()=1
   set_angular_force(1,0,-(Movespeed(1)/turn#),0)
   set_angular_force(2,0,-(Movespeed(1)/turn#),0)
endif
if rightkey()=1
   set_angular_force(1,0,(Movespeed(1)/turn#),0)
   set_angular_force(2,0,(Movespeed(1)/turn#),0)
endif
if upkey()=0 and keystate(17)=0 then set_linear_force(5,0,0,0)
if leftkey()=0 and rightkey()=0 then set_angular_force(1,0,0,0) : set_angular_force(2,0,0,0)
Set_Next_Angle_Y(5,1,10.0)
For w=3 to 4
   Set_Next_Angle_Y(w,5,10.0)
Next w

`Update Physics
for w=1 to wheels
   Set_Linear_Velocity(w,VelX(5),VelY(5),VelZ(5))
   Position_At_Object(w,5)
   Place_On_Matrix(w,1,Object Size Y(w)/2)
next w
Place_On_Matrix(5,1,(Object Size Y(5)/2)+(Object Size Y(5)/5))
Calculate_Offsets_AB(1,5,0,315,60) : Calculate_Offsets_AB(2,5,0,45,60)
Calculate_Offsets_AB(3,5,0,135,60) : Calculate_Offsets_AB(4,5,0,215,60)
for w=1 to 5
Get_Old_Position(w)
Get_Object_Direction(w)
Get_Traction_Force(w)
Get_Drag_Force(w)
Get_Force(w)
Get_Speed(w)
Get_Acceleration(w)
Set_Acceleration(w)
Set_Velocity(w)
Get_Angular_Acceleration(w)
Get_Angular_Velocity(w)
Update_Coordinates_Slide(w)
Update_Angles(w)
Get_Lateral_Forces(w)
Fully_Update_Object_Position(w)
Update_Object_Rotation(w)
Get_Movement_Angle(w)
next w
if upkey()=0 and keystate(17)=0
   if abs(MoveAngY(1)-AngleY(1)) < abs((MoveAngY(1)+180)-AngleY(1))
      angleY(1)=curveangle(MoveAngY(1),angleY(1),10.0)
      angleY(2)=curveangle(MoveAngY(1),angleY(2),10.0)
   else
      angleY(1)=curveangle(MoveAngY(1)+180,angleY(1),10.0)
      angleY(2)=curveangle(MoveAngY(1)+180,angleY(2),10.0)
   endif
endif


`Camera
ca#=curveangle(angleY(5),ca#,20.0)
cx#=positionX(5)-sin(ca#)*600
cz#=positionZ(5)-cos(ca#)*600
cy#=positionY(5)+300
position camera cx#,cy#,cz#
point camera positionX(5),positionY(5),positionZ(5)

`Info
ink rgb(255,255,255),0
Update_Speed(5) : Get_Movement_Time(5) : Get_Movement_Distance(5)
Get_Kinetic_Energy(5) : Get_Work_Done(5) : Get_Excerted_Power(5)
Text 10,10,"Kinetic Energy: "+str$(KineticEnergy(5))
Text 10,30,"Work Done: "+str$(WorkDone(5))
Text 10,50,"Excerted Power: "+str$(Power(5))
Text 1020-text width(str$(screen fps())),740,str$(screen fps())

`**End loop**
Sync
Loop

`**Functions**
Function Start_World_Physics()
Dim Gravity(0) as float
Dim Drag(0) as Float
Dim Windx(0) as float
Dim Windy(0) as float
Dim Windz(0) as float
Endfunction

Function Start_Object_Physics(o)
Dim Positionx(o) as Float : Dim Positiony(o) as Float : Dim Positionz(o) as Float
Dim OldPosX(o) as Float : Dim OldPosY(o) as Float : Dim OldPosZ(o) as Float
Dim Anglex(o) as Float : Dim Angley(o) as Float : Dim Anglez(o) as Float
Dim MoveAngX(o) as Float : Dim MoveAngY(o) as Float : Dim MoveAngZ(o) as Float
Dim LongForceX(o) as Float : Dim LongForceY(o) as Float : Dim LongForceZ(o) as Float
Dim TracForceX(o) as Float : Dim TracForceY(o) as Float : Dim TracForceZ(o) as Float
Dim DragForceX(o) as Float : Dim DragForceY(o) as Float : Dim DragForceZ(o) as Float
Dim AngForceX(o) as Float : Dim AngForceY(o) as Float : Dim AngForceZ(o) as Float
Dim AngVelX(o) as Float : Dim AngVelY(o) as Float : Dim AngVelZ(o) as Float
Dim AngAccX(o) as Float : Dim AngAccY(o) as Float : Dim AngAccZ(o) as Float
Dim VelX(o) as Float : Dim VelY(o) as Float : Dim VelZ(o) as Float
Dim AccX(o) as Float : Dim AccY(o) as Float : Dim AccZ(o) as Float
Dim DirX(o) as Float : Dim DirY(o) as Float : Dim DirZ(o) as Float
Dim OffsetX(o) as Float : Dim OffsetY(o) as Float : Dim OffsetZ(o) as Float
Dim MoveTime(o) as Float : Dim MoveDist(o) as Float : Dim Power(o) as Float
Dim Movespeed(o) as Float : Dim MoveConstant(o) as Float : Dim Acceleration(o) as Float
Dim TractionConstant(o) as Float : Dim AngularTractionConstant(o) as Float
Dim Mass(o) as Float : Dim Weight(o) as Float
Dim KineticEnergy(o) as Float: Dim WorkDone(o) as Float
Dim LatForce(o) as Float : Dim PushX(o) as Float : Dim PushY(o) as Float : Dim PushZ(o) as Float
Dim Force(o) as Float
Endfunction

Function Set_Standard_Gravity()
Gravity(0)=9.8
Endfunction

Function Set_Standard_Drag()
Drag(0)=0.47
Endfunction

Function Set_Gravity(g#)
Gravity(0)=g#
Endfunction

Function Set_Drag(d#)
Drag(0)=d#
Endfunction

Function Set_Move_Constant(o,c#)
MoveConstant(o)=c#
Endfunction

Function Set_Traction_Constant(o,c#)
TractionConstant(o)=c#
Endfunction

Function Set_Angular_Traction_Constant(o,c#)
AngularTractionConstant(o)=c#
Endfunction

Function Set_Body_Mass(o,m#)
Mass(o)=m#
Endfunction

Function Weigh_Object(o)
Weight(o)=Mass(o)*Gravity(0)
Endfunction

Function Set_Linear_Force(o,x,y,z)
LongForceX(o)=x : LongForceY(o)=y : LongForceZ(o)=z
Endfunction

Function Set_Linear_Velocity(o,x,y,z)
VelX(o)=x : VelY(o)=y : VelZ(o)=z
Endfunction

Function Set_Angular_Force(o,x,y,z)
AngForceX(o)=x : AngForceY(o)=y : AngForceZ(o)=z
Endfunction

Function Set_Angular_Velocity(o,x,y,z)
AngVelX(o)=x : AngVelY(o)=y : AngVelZ(o)=z
Endfunction

Function Move_Object(o,spd)
VelX(o)=dirX(o)*spd
VelY(o)=dirY(o)*spd
VelZ(o)=dirZ(o)*spd
Endfunction

Function Push_Object(o,spd)
LongForceX(o)=dirX(o)*spd
LongForceY(o)=dirY(o)*spd
LongForceZ(o)=dirZ(o)*spd
Endfunction

Function Control_Object_Using_Forces(o,push,turn)
if upkey()=1 then Push_Object(o,push)
if downkey()=1 then Push_Object(o,-push)
if leftkey()=1 then set_angular_Force(o,0,-(MoveSpeed(o)/turn),0)
if rightkey()=1 then set_angular_Force(o,0,(MoveSpeed(o)/turn),0)
if upkey()=0 and downkey()=0 then Push_Object(o,0)
if leftkey()=0 and rightkey()=0 then set_angular_force(o,0,0,0)
Standard_Update(o)
Endfunction

Function Control_Object_Using_Velocities(o,spd,turn)
if upkey()=1 then Move_Object(o,spd)
if downkey()=1 then Move_Object(o,-spd)
if leftkey()=1 then set_angular_velocity(o,0,-turn,0)
if rightkey()=1 then set_angular_velocity(o,0,turn,0)
if upkey()=0 and downkey()=0 then Move_Object(o,0)
if leftkey()=0 and rightkey()=0 then set_angular_velocity(o,0,0,0)
Standard_Update(o)
Endfunction

Function Get_Movement_Angle(o)
MoveAngX(o)=AtanFull(VelY(o),VelZ(o))
MoveAngY(o)=AtanFull(VelX(o),VelZ(o))
MoveAngZ(o)=AtanFull(VelX(o),VelY(o))
Endfunction

Function Set_Next_Angle_X(a,b,smth#)
AngleX(a)=Curveangle(AngleX(b),AngleX(a),smth#)
Endfunction

Function Set_Next_Angle_Y(a,b,smth#)
AngleY(a)=Curveangle(AngleY(b),AngleY(a),smth#)
Endfunction

Function Set_Next_Angle_Z(a,b,smth#)
AngleZ(a)=Curveangle(AngleZ(b),AngleZ(a),smth#)
Endfunction

Function Position_At_Object(a,b)
PositionX(a)=PositionX(b)
PositionY(a)=PositionY(b)
PositionZ(a)=PositionZ(b)
Endfunction

Function Get_Old_Position(o)
OldPosX(o)=PositionX(o)-VelX(o)
OldPosY(o)=PositionY(o)-VelY(o)
OldPosZ(o)=PositionZ(o)-VelZ(o)
Endfunction

Function Get_Object_Direction(o)
DirX(o)=sin(AngleY(o))*cos(AngleX(o))
DirY(o)=-sin(AngleX(o))
DirZ(o)=cos(AngleY(o))*cos(AngleX(o))
Endfunction

Function Get_Traction_Force(o)
TracForceX(o)=DirX(o)*MoveConstant(o)
TracForceY(o)=DirY(o)*MoveConstant(o)
TracForceZ(o)=DirZ(o)*MoveConstant(o)
Endfunction

Function Get_Drag_Force(o)
DragForceX(o)=-Drag(0)*VelX(o)*MoveSpeed(o)
DragForceY(o)=-Drag(0)*VelY(o)*MoveSpeed(o)
DragForceZ(o)=-Drag(0)*VelZ(o)*MoveSpeed(o)
Endfunction

Function Get_Force(o)
Force(o)=Sqrt(LongForceX(o)^2+LongForceY(o)^2+LongForceZ(o)^2)
Endfunction

Function Get_Speed(o)
MoveSpeed(o)=sqrt(VelX(o)^2 + VelY(o)^2 + VelZ(o)^2)
Endfunction

Function Get_Acceleration(o)
Acceleration(o)=Sqrt(AccX(o)^2+AccY(o)^2+AccZ(o)^2)
Endfunction

Function Set_Acceleration(o)
AccX(o)=LongForceX(o)/Mass(o)
AccY(o)=LongForceY(o)/Mass(o)
AccZ(o)=LongForceZ(o)/Mass(o)
Endfunction

Function Set_Velocity(o)
VelX(o)=VelX(o)+AccX(o) : VelX(o)=VelX(o)*TractionConstant(o)
VelY(o)=VelY(o)+AccY(o) : VelY(o)=VelY(o)*TractionConstant(o)
VelZ(o)=VelZ(o)+AccZ(o) : VelZ(o)=VelZ(o)*TractionConstant(o)
Endfunction

Function Get_Angular_Acceleration(o)
AngAccX(o)=AngForceX(o)/Mass(o)
AngAccY(o)=AngForceY(o)/Mass(o)
AngAccZ(o)=AngForceZ(o)/Mass(o)
Endfunction

Function Get_Angular_Velocity(o)
AngVelX(o)=AngVelX(o)+AngAccX(o) : AngVelX(o)=AngVelX(o)*AngularTractionConstant(o)
AngVelY(o)=AngVelY(o)+AngAccY(o) : AngVelY(o)=AngVelY(o)*AngularTractionConstant(o)
AngVelZ(o)=AngVelZ(o)+AngAccZ(o) : AngVelZ(o)=AngVelZ(o)*AngularTractionConstant(o)
Endfunction

Function Update_Coordinates_Straight(o)
PositionX(o)=PositionX(o)+(DirX(o)*Movespeed(o))
PositionY(o)=PositionY(o)+(DirY(o)*Movespeed(o))
PositionZ(o)=PositionZ(o)+(DirZ(o)*Movespeed(o))
Endfunction

Function Update_Coordinates_Slide(o)
PositionX(o)=PositionX(o)+VelX(o)
PositionY(o)=PositionY(o)+VelY(o)
PositionZ(o)=PositionZ(o)+VelZ(o)
Endfunction

Function Update_Angles(o)
AngleX(o)=Wrapvalue(AngleX(o)+AngVelX(o))
AngleY(o)=Wrapvalue(AngleY(o)+AngVelY(o))
AngleZ(o)=Wrapvalue(AngleZ(o)+AngVelZ(o))
Endfunction

Function Get_Lateral_Forces(o)
LatForce(o)=abs(AngVelY(o))*(MoveSpeed(o)/Mass(o))
If AngVelY(o)>0.0
   PushX(o)=(sin(Wrapvalue(AngleY(o)-90))*LatForce(o))
   PushZ(o)=(cos(Wrapvalue(AngleY(o)-90))*LatForce(o))
endif
If AngVelY(o)<0.0
   PushX(o)=(sin(Wrapvalue(AngleY(o)+90))*LatForce(o))
   PushZ(o)=(cos(Wrapvalue(AngleY(o)+90))*LatForce(o))
endif
Endfunction

Function Calculate_Offsets(o,sax#,say#,out#)
OffsetX(o)=(Sin(Wrapvalue(AngleY(o)+say#))*Cos(Wrapvalue(AngleX(o)+sax#)))*Out#
OffsetY(o)=-Sin(Wrapvalue(AngleX(o)+sax#))*out#
OffsetZ(o)=(Cos(Wrapvalue(AngleY(o)+say#))*Cos(Wrapvalue(AngleX(o)+sax#)))*out#
Endfunction

Function Calculate_Offsets_AB(a,b,sax#,say#,out#)
OffsetX(a)=(Sin(Wrapvalue(AngleY(b)+say#))*Cos(Wrapvalue(AngleX(b)+sax#)))*out#
OffsetY(a)=-Sin(Wrapvalue(AngleX(b)+sax#))*out#
OffsetZ(a)=(Cos(Wrapvalue(AngleY(b)+say#))*Cos(Wrapvalue(AngleX(b)+sax#)))*out#
Endfunction

Function Place_On_Matrix(o,m,h#)
PositionY(o)=Get ground height(m,positionX(o),PositionZ(o))+h#
Endfunction

Function Update_Object_Position(o)
Position Object o,PositionX(o),PositionY(o),PositionZ(o)
Endfunction

Function Update_Position_With_Push_Forces(o,a,b,c)
Position object o,Positionx(o)+(PushX(o)*a),Positiony(o)+(PushY(o)*b),Positionz(o)+(PushZ(o)*c)
Endfunction

Function Update_Position_With_Offsets(o,a,b,c)
Position Object o,Positionx(o)+(offsetx(o)*a),Positiony(o)+(offsety(o)*b),PositionZ(o)+(offsetZ(o)*c)
Endfunction

Function Fully_Update_object_Position(o)
Position Object o,PositionX(o)+PushX(o)+OffsetX(o),PositionY(o)+PushY(o)+OffsetY(o),PositionZ(o)+PushZ(o)+OffsetZ(o)
Endfunction

Function Update_object_Rotation(o)
Rotate object o,AngleX(o),AngleY(o),AngleZ(o)
Endfunction

Function Standard_Update(o)
Get_Old_Position(o)
Get_Object_Direction(o)
Get_Traction_Force(o)
Get_Drag_Force(o)
Get_Force(o)
Get_Speed(o)
Get_Acceleration(o)
Set_Acceleration(o)
Set_Velocity(o)
Get_Angular_Acceleration(o)
Get_Angular_Velocity(o)
Update_Coordinates_Straight(o)
Update_Angles(o)
Get_Lateral_Forces(o)
Fully_Update_Object_Position(o)
Update_Object_Rotation(o)
Endfunction

Function Update_Speed(o)
if MoveSpeed(o)<0.9 then MoveSpeed(o)=MoveSpeed(o)*0.0
Endfunction

Function Get_Movement_Time(o)
If MoveSpeed(o)<>0 then MoveTime(o)=Movetime(o)+1.0 else MoveTime(o)=0.0
Endfunction

Function Get_Movement_Distance(o)
MoveDist(o)=MoveSpeed(o)*MoveTime(o)
Endfunction

Function Get_Kinetic_Energy(o)
KineticEnergy(o)=(Mass(o)/2)*(Movespeed(o)^2)
Endfunction

Function Get_Work_Done(o)
WorkDone(o)=Force(o)*MoveDist(o)
Endfunction

Function Get_Excerted_Power(o)
If MoveTime(o)>0.0 then Power(o)=KineticEnergy(o)/MoveTime(o) else Power(o)=0.0
Endfunction

SYNC ON : SYNC RATE 60 : HIDE MOUSE
SET DISPLAY MODE 1024,768,32

SET TEXT SIZE 28 : SET TEXT FONT "Impact"
SYNC : CENTER TEXT 512,350,"Loading..." : SYNC

INK RGB(200,200,200),0
BOX 1,1,128,128
FOR I=0 TO 2000
   COLR=RND(255)
   INK RGB(COLR,COLR,COLR),0
   DOT RND(128),RND(128)
NEXT I
GET IMAGE 1,1,1,128,128

CREATE BITMAP 1,128,128
SET CURRENT BITMAP 1
FOR I=1 TO 16
   FOR J=1 TO 16
      IF RND(1)=1
         INK RGB(255,0,0),0
      ELSE
         INK RGB(255,255,0),0
      ENDIF
      BOX ((J*8)-8),((I*8)-8),(J*8),(I*8)
   NEXT J
NEXT I
BLUR BITMAP 1,6
GET IMAGE 2,1,1,128,128
SET CURRENT BITMAP 0
DELETE BITMAP 1

MAKE MATRIX 1,5000,5000,5,5
PREPARE MATRIX TEXTURE 1,1,1,1
SET MATRIX TEXTURE 1,1,1
FOR I=1 TO 4
   SET MATRIX HEIGHT 1,I,1,1000
   SET MATRIX HEIGHT 1,I,4,1000
NEXT I
FOR I=2 TO 3
   SET MATRIX HEIGHT 1,1,I,1000
   SET MATRIX HEIGHT 1,2,I,1000
   SET MATRIX HEIGHT 1,3,I,1000
   SET MATRIX HEIGHT 1,4,I,1000
NEXT I
SET MATRIX HEIGHT 1,2,2,0
SET MATRIX HEIGHT 1,3,3,0
SET MATRIX HEIGHT 1,1,4,0
SET MATRIX HEIGHT 1,4,1,0
UPDATE MATRIX 1

LAVA=2
MAKE MATRIX LAVA,7500,7500,50,50
PREPARE MATRIX TEXTURE LAVA,2,1,1
SET MATRIX TEXTURE LAVA,1,1
SET MATRIX LAVA,0,0,1,1,1,1,1

PLR=1
MAKE OBJECT SPHERE PLR,50,12,12
SET OBJECT COLLISION ON PLR

FOG ON : FOG COLOR RGB(255,60,0) : FOG DISTANCE 1500
BACKDROP ON : COLOR BACKDROP RGB(255,60,0)

START:
PLRA=90 : PLRY#=GET GROUND HEIGHT(1,1000,1000)+25
DIM LAVASIN(50,50)
COLOR OBJECT PLR,RGBB(200)
POSITION OBJECT PLR,1000,PLRY#,1000
YROTATE OBJECT PLR,PLRA

FOR I=1 TO 50
   FOR J=1 TO 50
      LAVASIN(I,J)=RND(360)
   NEXT J
NEXT I

POSITION MATRIX LAVA,-1250,0,-1250

DO
   IF UPKEY()=1 THEN MOVE OBJECT PLR,20
   IF LEFTKEY()=1 THEN PLRA=PLRA-2
   IF RIGHTKEY()=1 THEN PLRA=PLRA+2
   YROTATE OBJECT PLR,PLRA
   PLRY#=GET GROUND HEIGHT(1,OBJECT POSITION X(PLR),OBJECT POSITION Z(PLR))+25
   POSITION OBJECT PLR,OBJECT POSITION X(PLR),PLRY#,OBJECT POSITION Z(PLR)
   SET CAMERA TO OBJECT ORIENTATION PLR
   POSITION CAMERA OBJECT POSITION X(PLR),PLRY#+100,OBJECT POSITION Z(PLR)
   MOVE CAMERA -175 : POINT CAMERA OBJECT POSITION X(PLR),PLRY#+15,OBJECT POSITION Z(PLR)

   POSITION MATRIX 1,RND(5),RND(5),RND(5)

   POSITION MATRIX LAVA,-1250,MATRIX POSITION Y(LAVA)+1,-1250
   IF OBJECT POSITION Y(PLR)<=MATRIX POSITION Y(LAVA) THEN GOSUB DEATH

   FOR I=1 TO 50
      FOR J=1 TO 50
         INC LAVASIN(I,J),5
         IF LAVASIN(I,J)>360 THEN LAVASIN(I,J)=0
         SET MATRIX HEIGHT LAVA,I,J,SIN(LAVASIN(I,J))*50
      NEXT J
   NEXT I
   UPDATE MATRIX LAVA

   SYNC
LOOP

DEATH:
   INK RGB(255,255,255),0
   COLOR OBJECT PLR,0 : SYNC
   CENTER TEXT 512,350,"You Died" : SYNC
   SLEEP 1000 : CLS
   GOSUB START

Rem Project: Gt_Racer
Rem Created: 27/12/2005 12:05:37

Rem ***** Main Source File *****
rem
sync on : sync rate 60
`load object "race track.x", 2
`load image "gt4 works.jpg", 2 , 1
make matrix 1, 100, 100, 10, 10

`position object 2, 0, -3, 0
`load object "gt4.x", 1
make object cube 1, 10
show object 1
`texture object 1, 2
position object 1, 0, -3, 5
make camera 1

do
x#=object position x(1)
y#=object position y(1)
z#=object position z(1)
oldx#=x#
oldy#=y#
oldz#=z#
if upkey()=1 then xspeed#=sin(ay#)*30 : zspeed#=cos(ay#)*30
if downkey()=1 then xspeed#=sin(ay#)*-30 : zspeed#=cos(ay#)*-30
if rightkey()=1 then ay#=wrapvalue(ay#+2)
if leftkey()=1 then ay#=wrapvalue(ay#-2)
xspeed#=xspeed#*0.97
zspeed#=zspeed#*0.97
x#=x#+xspeed#
z#=z#+zspeed#
position object 1,x#,y#,z#
yrotate object 1,ay#
speed#=sqrt((x#-oldx#)^2+(y#-oldy#)^2+(z#-oldz#)^2)
Create_Speedometer(speed#,100,100,4,20,rgb(255,0,255),rgb(0,255,0))
ca#=curveangle(ay#,ca#,20.0)
cx#=x#-sin(ca#)*300
cz#=z#-cos(ca#)*300
cy#=curvevalue(y#+100,cy#,20.0)
position camera cx#,cy#,cz#
point camera x#,y#,z#
position camera 1,cx#,cy#,cz#
yrotate camera 1,wrapvalue(camera angle y(0)+180)
set camera to image 1,10,128,64
sprite 10,screen width()/2-64,30,10
sync
loop

Function create_speedometer(movespeed#,speedox,speedoy,speedow,speedl,col1,col2)

ink col1,0
line speedox,speedoy,speedox+sin(wrapvalue(int(360-movespeed#)+180))*(speedow/2),speedoy+cos(wrapvalue(int(360-movespeed#)+180))*(speedow/2)
line speedox,speedoy,speedox-sin(wrapvalue(int(360-movespeed#)+180))*(speedow/2),speedoy-cos(wrapvalue(int(360-movespeed#)+180))*(speedow/2)
line speedox+sin(wrapvalue(int(360-movespeed#)+180))*(speedow/2),speedoy+cos(wrapvalue(int(360-movespeed#)+180))*(speedow/2),speedox+sin(wrapvalue(int(360-movespeed#)+270))*speedl,speedoy+cos(wrapvalue(int(360-movespeed#)+270))*speedl
line speedox-sin(wrapvalue(int(360-movespeed#)+180))*(speedow/2),speedoy-cos(wrapvalue(int(360-movespeed#)+180))*(speedow/2),speedox+sin(wrapvalue(int(360-movespeed#)+270))*speedl,speedoy+cos(wrapvalue(int(360-movespeed#)+270))*speedl
ink col2,0
text speedox-(speedl*1.5),speedoy,"0"
text speedox,speedoy-(speedl*1.5),"90"
text speedox+(speedl*1.5),speedoy,"180"
text speedox,speedoy+(speedl*1.5),"270"

endfunction

`Racing Challenge by Ric

detail=40

`display settings
if check display mode(1024,768,32)=1 then set display mode 1024,768,32
sync on
sync rate 50
hide mouse

autocam off
hide light 0
color backdrop 0
set ambient light 0
make light 1
set light range 1,3
`set initial value for image numbers
global image_seed
image_seed=1000000

create bitmap 1,screen width(),screen height()

position light 1,0,0,-1
set light range 1,2

make object sphere 1,1,100,100
set object specular power 1,20
set object specular 1,rgb(255,255,255)
position camera 0,0,-2
sync
sun=free_image()
get image sun,screen width()/4,screen height()/4,screen width()*.75,screen height()*.75,1
set object light 1,0
sync
flare1=free_image()
get image flare1,screen width()/4,screen height()/4,screen width()*.75,screen height()*.75,1
delete object 1

delete bitmap 1
text 0,0,"Please wait ......"
sync

`camera and light settings

set camera range 0.01,10000
position camera 64,0,64
color backdrop 0
set ambient light 20


`set size of terrain
global matrixsize
matrixsize=128

`fog settings
fog on
fog color rgb(200,200,200)
fog distance matrixsize

`set initial value for object numbers
global object_seed
object_seed=1000000



`position light for world illumination
position light 1,200,30,200
set light range 1,100000
make object cube 2000,1
hide object 2000
position object 2000,light position x(1),light position y(1),light position z(1)

global skyimage
create_sky()


global height#
global rows
global columns
global terrainimage
`set size of terrain
rows=detail
columns=detail
global peaksoff
global peaksoff2
global dim peakx(100)
global dim peakz(100)
global dim peakh(100)
global dim peakx2(100)
global dim peakz2(100)
global dim peakh2(100)
terraform()
create_trench()
Matrix_to_Object(1000,1,128,128,rows,rows,rows,rows,rows,rows)


ink rgb(255,255,255),0

gosub calculate_shadows
delete object 1000
set current bitmap 0

global sphere1
global sphere2
global sphere3
create_skysphere()

sunplane=free_object()
make object plain sunplane,160*4/3.0,160
set object fog sunplane,0
texture object sunplane,sun
set object light sunplane,0
ghost object on sunplane,2
set object transparency sunplane,2
flareplane1=free_object()
make object plain flareplane1,0.3*4/3.0,0.3
texture object flareplane1,flare1
set object light flareplane1,0
set alpha mapping on flareplane1,20
position object flareplane1,0,-0.5,0
`reference object for distant light source
dummy=free_object()
make object cube dummy,1
hide object dummy
point object dummy,light position x(1),light position y(1),light position z(1)
move object dummy,10000



position camera 64,height#+2+get ground height(1,64,64),64

`boaty and watery stuff


`set mouse position to centre
position mouse 0.5*screen width(),0

 camdist#=4 `cam distance
 camheight#=2 `cam height
 smoothness=50 `chasecam smoothness


 `wavegenerator1 x,z coordinates:
 a=15
 b=10

 `wavegenerator2 x,z coordinates:
 c=35
 d=40

 `wave characteristics
 elasticity#=0.001 `determines wave speed
 damping#=0.992 `determines the decay of the wave height
 amplitude#=2.5  `the wave height at the generator points
 frequency#=2.0   `how quickly the generator points move up and down
 current#=0.01  `how much the wave slope affects the jetski

`******************

`***jetski motion properties***

friction#=0.975 `1=no friction - slides all over the place
               `0=big friction - like moving through concrete

power#=0.007    `controls the acceleration of the jetski


dim tile(rows+10,columns+10)
dim v#(rows*2,columns*2)

make matrix 2,rows*4,columns*4,rows,columns


`give every tile of the matrix a coordinate - tile(x,z)
for x=1 to rows
for z=1 to columns
tile(x,z)=tile
tile=tile+1
next z
next x

`fix the wave generator points
dim fixed(rows*columns*2)
fixed(tile(a,b))=1
fixed(tile(c,d))=1

create bitmap 1,10,10
ink rgb(0,0,255),0
box 0,0,10,10
`ink rgb(0,0,180),0
`box 0,0,9,9
get image 1,0,0,10,10,1
prepare matrix texture 2,1,1,1
ghost matrix on 2,3


`Make a simple jetski - object number 1
  make object cube 1,0.1
  make object box 2,2,1,4
  make mesh from object 1,2
  delete object 2
  make object cone 3,2
  scale object 3,100,200,100
  xrotate object 3,90
  make mesh from object 2,3
  delete object 3
  make object sphere 4,2
  make mesh from object 3,4
  delete object 4
  add limb 1,1,1
  add limb 1,2,2
  offset limb 1,2,0,0,2
  add limb 1,3,3
  ink rgb(250,0,50),0
  box 0,0,10,10
  get image 2,0,0,10,10
  set current bitmap 0
  texture object 1,2
  ink rgb(255,255,255),0
  scale object 1,20,20,20

`*****************************************************



`*******initial positions and speeds *********


number=50  `number of floating boxes
dim x#(number)
dim y#(number)
dim z#(number)
dim frontspeed#(number)
dim sidespeed#(number)

x#(1)=64  `initial x position of jetski
z#(1)=96  `initial z position of jetski

for object=2 to number

make object cube object,5
hide object object

x#(object)=rnd(500)  `initial x position of boxes
z#(object)=rnd(500)  `initial z position of boxes

next object

`**********************************


`main program loop
do

position object sunplane,camera position x(),camera position y(),camera position z()
point object sunplane,object position x(dummy),object position y(dummy),object position z(dummy)
move object sunplane,128*7
point object sunplane,camera position x(),camera position y(),camera position z()

gosub wave

gosub control

gosub get_wave_slope

gosub chasecam

text 0,0,str$(screen fps())

sync
loop

`****************************************

`functions

function terraform

`create terrain
make matrix 1,matrixsize,matrixsize,columns,rows
set matrix 1,0,0,0,2,1,1,0

for peak=1 to 50
peakh(peak)=rnd(5)+5
peakx(peak)=rnd(rows-8)+4
peakz(peak)=rnd(rows-8)+4
set matrix height 1,peakx(peak),peakz(peak),peakh(peak)
next peak
update matrix 1


for times=1 to 120
if times>=90 then peaksoff=1
smooth_matrix(1,rows,rows)
next times
peaksoff=0


endfunction

function create_trench()

data$="663663222222222222121477441114788888889698744444478998888963336666"
x=20
z=30
length=len(data$)
for n=1 to length
check$=right$(left$(data$,n),1)
if check$="1" then dec x:dec z
if check$="2" then dec z
if check$="3" then inc x:dec z
if check$="4" then dec x
if check$="6" then inc x
if check$="7" then inc z:dec x
if check$="8" then inc z
if check$="9" then inc z:inc x

for u=-1 to 1
for v=-1 to 1
set matrix height 1,x+u,z+v,-2
next v
next u
next n


`fill with water
water=free_object()
make object plain water,4000,4000
xrotate object water,-90
ghost object on water,3
set object fog water,0
position object water,0,-1.8,0
color object water,rgb(0,0,255)


water2=free_object()
make object plain water2,4000,4000
xrotate object water2,-90
`ghost object on water2,3
set object fog water2,0
position object water2,0,-2.85,0
color object water2,rgb(0,0,255)


endfunction


function smooth_matrix(matnum,tilex,tilez)

for x=2 to tilex-2
for z=2 to tilez-2
b=-get matrix height(matnum,x,z)+get matrix height(matnum,x,z+1)
d=-get matrix height(matnum,x,z)+get matrix height(matnum,x+1,z)
f=-get matrix height(matnum,x,z)+get matrix height(matnum,x,z-1)
h=-get matrix height(matnum,x,z)+get matrix height(matnum,x-1,z)
total=b+d+f+h

set matrix height matnum,x,z,get matrix height(matnum,x,z)+total*0.1

next z
next x



for peak=1 to 100
`peakh(peak)=peakh(peak)*0.999
if matnum=1
  if peaksoff=0 then set matrix height 1,peakx(peak),peakz(peak),peakh(peak)
endif
if matnum=2
  if peaksoff=0 then set matrix height 2,peakx2(peak),peakz2(peak),peakh2(peak)
endif
next peak

update matrix matnum
endfunction

function free_object

object=object_seed
repeat
inc object
until object exist(object)=0

endfunction object

function check_free_object_number

object=object_seed
repeat
inc object
until object exist(object)=0

endfunction object

function free_image

image=image_seed
repeat
inc image
until image exist(image)=0

endfunction image



function create_skysphere

`get image for sphere by taking snapshot of terrain
terrainimage=free_image()
set camera fov 90
fog off
position camera matrixsize/2.0,3,-matrixsize*0.3
sync
get image terrainimage,0,0,screen width(),screen height()
fog on
fog color rgb(200,200,200)
fog distance matrixsize
set camera fov 60


sphere2=free_object()
make object sphere sphere2,matrixsize*10
set object collision off sphere2
xrotate object sphere2,180
`yrotate object sphere2,30
set object cull sphere2,0
texture object sphere2,terrainimage
set object texture sphere2,2,0
scale object texture sphere2,8,2.6
set object fog sphere2,0
set object transparency sphere2,2
ghost object on sphere2,5
set alpha mapping on sphere2,20
set object emissive sphere2,rgb(80,60,100)
`set object specular power sphere2,100
`fade object sphere2,200
`set object light sphere2,0
position object sphere2,matrixsize/2,-135,matrixsize/2

`sphere 3 (sky) - a third, outer sphere for the sky
sphere3=free_object()
make object sphere sphere3,matrixsize*20
set object collision off sphere3
set object cull sphere3,0
texture object sphere3,skyimage
scale object texture sphere3,1,1.8
set object light sphere3,0
set object fog sphere3,0
position object sphere3,matrixsize/2,0,matrixsize/2

endfunction

function create_sky()
`creates a dark to light bluish gradient
cls
for n=0 to 250
  ink rgb(n,n-20,n-40),0
  line 0,n,250,n
next n

skyimage=free_image()
get image skyimage,0,0,250,250

endfunction







function move_camera

control camera using arrowkeys 0,0.2,3
position camera camera position x(),height#+.2+get ground height(1,camera position x(),camera position z()),camera position z()
if inkey$()="a" then inc height#,1
if inkey$()="z" then dec height#,1
endfunction

calculate_shadows:
undim pixelshaded1(128,128)
undim pixelshaded2(128,128)
dim pixelshaded1(128,128)
dim pixelshaded2(128,128)
`if light1>0
create bitmap 1,128,128
set current bitmap 1
ink rgb(20,155,0),0
box 0,0,128,128
for n=0 to 5000
tone=rnd(100)+150
ink rgb(tone/10,tone+15,tone/15),0
dot rnd(128),rnd(128)
next n
ink rgb(10,10,10),0
for y=0 to 128
for x=0 to 128


`shadows cast by terrain - matrix converted to object 1000
ray2#=intersect object(1000,x,get ground height(1,x,y)+0.1,y,light position x(1),light position y(1),light position z(1))
if ray2#=0 then ray2#=1000
if ray2#<0 then ray2#=1
if ray2#<=100
  tone=40+ray2#*8

  `if pixelshaded2(x,y)=1 or pixelshaded1(x,y)=1 then tone=tone/1.5
  if tone>155 then tone=155
  ink rgb(20,tone,0),0
  dot x,128-y
  `pixelshaded1(x,y)=1
endif


next x
next y
blur bitmap 1,3
shadowmap=5000000
get image shadowmap,0,0,128,128

delete bitmap 1
`set light mapping on 2,shadowmap
 `(light1 exist)

`texture matrix
prepare matrix texture 1,shadowmap,rows,columns
tile=1
for x=rows-1 to 0 step -1
for z=0 to columns-1
set matrix tile 1,z,x,tile
inc tile
next z
next x
`endif
update matrix 1

return

function Matrix_to_Object( object, matrixnum, matxsize#, matzsize#, matxsegs, matzsegs, tilex, tilez, limbx, limbz )
   `safety in case limbx and limbz is smaller than the texture x and z
   a = 0
   if limbx < tilex then limbx = tilex : a = 1
   if limbz < tilez then limbz = tilez : a = 1
   if a > 0
      sync
      print "limbx or limbz was less than tilex or tilez : Error corrected"
      print "Press any key to continue."
      sync
      wait key
   endif
   `safety in case matxsegs and matzsegs is not evenly divisible by tilex and tilez
   if matxsegs mod tilex > 0
      EXIT PROMPT "Number of matrix xsegs not evenly divisible by tilex", "Texture tilex error"
      end
   endif
   if matzsegs mod tilez > 0
      EXIT PROMPT "Number of matrix zsegs not evenly divisible by tilez", "Texture tilez error"
      end
   endif
   `safety in case matxsegs and matzsegs is not evenly divisible by tilex and tilez
   if matxsegs mod limbx > 0
      EXIT PROMPT "Number of matrix xsegs not evenly divisible by limbx", "Texture tilex error"
      end
   endif
   if matzsegs mod limbz > 0
      EXIT PROMPT "Number of matrix zsegs not evenly divisible by limbz", "Texture tilez error"
      end
   endif

   `calc the number of polys in the matrix
   num_mat_polys = (matxsegs*matzsegs)*2
   `make array to store vert info
   `storage indexs = poly number, vert number, x/y/z vert pos/norm/uv
   `polys are numbered starting from left to right, front to
   `back, one row at a time (with 2 polys per tile)
   dim vert_store#(num_mat_polys,2,7)

   `calc the width and depth of each tile on the matrix
   `formula size#/number
   mat_x_wide# = matxsize#/matxsegs
   mat_z_deep# = matzsize#/matzsegs


`******************************************************************************
   `get vert position data
   `placeholder variable for poly numbers
   a = 1
   `one row at a time (front to back)
   for j = 0 to matzsegs - 1
      `one tile at a time (left to right)
      for i = 0 to matxsegs - 1

         `tile top left poly info

         `bottom left vert (x,y,z)
         vert_store#(a,0,0) = i*mat_x_wide#
         vert_store#(a,0,1) = GET MATRIX HEIGHT(matrixnum, i, j)
         vert_store#(a,0,2) = j*mat_z_deep#
         `top left vert (x,y,z)
         vert_store#(a,1,0) = i*mat_x_wide#
         vert_store#(a,1,1) = GET MATRIX HEIGHT(matrixnum, i, j+1)
         vert_store#(a,1,2) = (j+1)*mat_z_deep#
         `top right vert (x,y,z)
         vert_store#(a,2,0) = (i+1)*mat_x_wide#
         vert_store#(a,2,1) = GET MATRIX HEIGHT(matrixnum, i+1, j+1)
         vert_store#(a,2,2) = (j+1)*mat_z_deep#
         inc a, 1

         `tile bottom right poly info

         `bottom left vert (x,y,z)
         vert_store#(a,0,0) = i*mat_x_wide#
         vert_store#(a,0,1) = GET MATRIX HEIGHT(matrixnum, i, j)
         vert_store#(a,0,2) = j*mat_z_deep#
         `top right vert (x,y,z)
         vert_store#(a,1,0) = (i+1)*mat_x_wide#
         vert_store#(a,1,1) = GET MATRIX HEIGHT(matrixnum, i+1, j+1)
         vert_store#(a,1,2) = (j+1)*mat_z_deep#
         `bottom right vert (x,y,z)
         vert_store#(a,2,0) = (i+1)*mat_x_wide#
         vert_store#(a,2,1) = GET MATRIX HEIGHT(matrixnum, i+1, j)
         vert_store#(a,2,2) = j*mat_z_deep#
         inc a, 1
      next i
   next j
`******************************************************************************


`******************************************************************************
   `calc normals for polys
   `Thanks to ADR for posting this code on the DBP forums :)
   for i = 1 to num_mat_polys
      `acuire vert positions
      P1X# = vert_store#(i,0,0)
      P1Y# = vert_store#(i,0,1)
      P1Z# = vert_store#(i,0,2)
      P2X# = vert_store#(i,1,0)
      P2Y# = vert_store#(i,1,1)
      P2Z# = vert_store#(i,1,2)
      P3X# = vert_store#(i,2,0)
      P3Y# = vert_store#(i,2,1)
      P3Z# = vert_store#(i,2,2)

      null = make vector3(1)
      null = make vector3(2)
      null = make vector3(3)

      ` -- calculate the two directional vectors for the adj and opp edges...
      set vector3 1, P1X#, P1Y#, P1Z#
      set vector3 2, P2X#, P2Y#, P2Z#
      set vector3 3, P3X#, P3Y#, P3Z#
      subtract vector3 2, 2, 1
      subtract vector3 3, 3, 1   ` -- vector 3 and 1 are now directional vectors
      normalize vector3 2,2      ` -- normalize em
      normalize vector3 3,3
      cross product vector3 1, 2,3  ` -- use the origin vector (1) to store the face normal
      normalize vector3 1,1

      `save normals (all 3 verts have same normals)
      vert_store#(i,0,3) = x vector3(1)
      vert_store#(i,0,4) = y vector3(1)
      vert_store#(i,0,5) = z vector3(1)
      vert_store#(i,1,3) = vert_store#(i,0,3)
      vert_store#(i,1,4) = vert_store#(i,0,4)
      vert_store#(i,1,5) = vert_store#(i,0,5)
      vert_store#(i,2,3) = vert_store#(i,0,3)
      vert_store#(i,2,4) = vert_store#(i,0,4)
      vert_store#(i,2,5) = vert_store#(i,0,5)

      null = delete vector3(1)
      null = delete vector3(2)
      null = delete vector3(3)
   next i
`******************************************************************************


`******************************************************************************
   `calc UV data for polys

   `save current x tile number
   xtiles = 1
   `save current z tile number
   ztiles = 1
   `calc how much to step each u data per tile
   stepu# = (1.0/tilex)
   `calc how much to step each v data per tile
   stepv# = (1.0/tilez)
   `set base u data for new set of tiles
   baseu# = 0
   `set base v data for new set of tiles
   basev# = 1-stepv#
   `poly number placeholder variable
   i = 1
   `from front to back
   for k = 1 to matzsegs
      `reset the number of x tiles to 1 and the u base to 0
      `at the beginning of each row
      xtiles = 1
      baseu# = 0
      `from left to right
      for l = 1 to matxsegs
         `2 polys per tile
         for m = 1 to 2
            `write all 3 verts of each matrix poly
            for j = 0 to 2
               `select which formula to apply depending on polygon side
               `and vert number (0 to 2)
               `j selects the vert number
               `m selects the polygon side (back/left or front right : 1 or 2)
               select j
                  case 0
                     if m = 1
                        testu# = baseu#
                        testv# = basev#+stepv#
                     else
                        testu# = baseu#
                        testv# = basev#+stepv#
                     endif
                  endcase
                  case 1
                     if m = 1
                        testu# = baseu#
                        testv# = basev#
                     else
                        testu# = baseu#+stepu#
                        testv# = basev#
                     endif
                  endcase
                  case 2
                     if m = 1
                        testu# = baseu#+stepu#
                        testv# = basev#
                     else
                        testu# = baseu#+stepu#
                        testv# = basev#+stepv#
                     endif
                  endcase
               endselect
               `store calculated data for each vert of each poly
               `u data
               vert_store#(i,j,6) = testu#
               `v data
               vert_store#(i,j,7) = testv#
            next j
            inc i, 1
         `next polygon side of this tile
         next m
         `update u data and xtiles place holder
         inc baseu#, stepu#
         inc xtiles, 1
         `reset data when texture tile width has been reached
         if xtiles > tilex
            baseu# = 0
            xtiles = 1
         endif
      `next x tile
      next l
      `update v data and ztiles place holder
      dec basev#, stepv#
      inc ztiles, 1
      `reset data when texture depth has been reached
      if ztiles > tilez
         basev# = 1-stepv#
         ztiles = 1
      endif
   `next z tile
   next k
`******************************************************************************


`******************************************************************************
   `make object from matrix verts

   `calc number of memblocks needed
   tempx = matxsegs/limbx
   tempz = matzsegs/limbz
   tempmem = tempx * tempz
   `make arrays to store memblock positions and memblock numbers
   dim membhold(tempx, tempz)
   dim mempos(tempmem)
   `enter memblock numbers for each texture tile based on the x and z limbs to use
   a = 1
   for i = 1 to tempz
      for j = 1 to tempx
         membhold(j, i) = a
         inc a, 1
      next j
   next i
   `calc the memblock size
   `formula is 12 byte header + ((( 32 bytes per vert * limb x segs)*(limb z segs * 2 polys))* 3 verts per poly)
   mat_mem_size = 12+(((32*limbx)*(limbz*2))*3)
   `make memblocks and write headers for each memblock and set the beginning
   `position for vert data to 12
   for i = 1 to tempmem
      make memblock i, mat_mem_size
      `write objects fvf format as 274
      write memblock dword i, 0, 274
      `write bytes per vert as 32 (8 floats xpos#,ypos#,zpos#,xnorm#,ynorm#,znorm#,u#,v#)
      write memblock dword i, 4, 32
      `write number of verts in matrix (polys*3)
      write memblock dword i, 8, ((limbx*limbz)*2)*3
      `set beginning position for vert data after header(0-8)+4 = 12
      mempos(i) = 12
   next i

   `image tile, row, and memblock placeholders
   xtile = 1
   ztile = 1
   xrow = 1
   xmem = 1
   zmem = 1
   `polygon placeholders (b=2 is 1 tile complete)
   b = 1
   `enter all polys' verts to memblock
   for i = 1 to num_mat_polys
      `select memblock number to write to use current x and z limb tile to pick
      tempmemnum = membhold(xmem, zmem)
      `write all 3 verts of each matrix poly
      for j = 0 to 2
         `xpos
         write memblock float tempmemnum, mempos(tempmemnum), vert_store#(i,j,0)
         `increment the current memblock position
         inc mempos(tempmemnum), 4
         `ypos
         write memblock float tempmemnum, mempos(tempmemnum), vert_store#(i,j,1)
         inc mempos(tempmemnum), 4
         `zpos
         write memblock float tempmemnum, mempos(tempmemnum), vert_store#(i,j,2)
         inc mempos(tempmemnum), 4
         `xnorm
         write memblock float tempmemnum, mempos(tempmemnum), vert_store#(i,j,3)
         inc mempos(tempmemnum), 4
         `ynorm
         write memblock float tempmemnum, mempos(tempmemnum), vert_store#(i,j,4)
         inc mempos(tempmemnum), 4
         `znorm
         write memblock float tempmemnum, mempos(tempmemnum), vert_store#(i,j,5)
         inc mempos(tempmemnum), 4
         `u data
         write memblock float tempmemnum, mempos(tempmemnum), vert_store#(i,j,6)
         inc mempos(tempmemnum), 4
         `v data
         write memblock float tempmemnum, mempos(tempmemnum), vert_store#(i,j,7)
         inc mempos(tempmemnum), 4
      next j
      `after each poly increase b by 1
      inc b, 1
      `after 2 polys have been completed 1 tile has been entered
      if b > 2
         `reset poly count to 1
         b = 1
         `move to next x tile
         inc xtile, 1
         `if at the next limb
         if xtile > limbx
            `reseet the x tile to 1
            xtile = 1
            `update the tile count to know when to move to the next row
            inc xrow, limbx
            `inc the xmem selector placeholder
            inc xmem, 1
            `if at the last xmem wide then start back at the begining
            if xmem > tempx
               xmem = 1
            endif
         endif
         `when at the end of the row
         if xrow > matxsegs
            `go back to the far left
            xrow = 1
            `update the number of rows done
            inc ztile, 1
            `when the number of rows done is > the preset depth
            if ztile > limbz
               `reset row count back to 1
               ztile = 1
               `move to the next zmem selector placeholder
               inc zmem, 1
               `if row count overlaps predetermined number
               if zmem > tempz
                  `wrap to 1
                  zmem = 1
               endif
            endif
         endif
      endif
   next i
   `undim arrays when done with them
   undim vert_store#(0,0,0)
   undim membhold(0, 0)
   undim mempos(0)
   `set the mesh number to use (can be replaced with a findfreemesh function)
   mesh = 1
   `make object and limbs from data setup
   for i = 1 to tempmem
      `make temp mesh from info
      make mesh from memblock mesh, i
      `delete temp memblock
      delete memblock i
      `if it is the first memblock mesh, use it as the base object
      if i = 1
         `make temp object for NGC
         make object object, mesh, 0
      `else turn each memblock mesh into limbs
      else
         add limb object, i-1, mesh
      endif
      `delete temp mesh
      delete mesh mesh
   next i
`******************************************************************************
endfunction



function check_free_image_number()

`this functionworks in the same way as the check_free_object function,
`except it works for images.

image=image_seed
repeat
inc image
until image exist(image)=0

endfunction image


wave:

`****move the generator points up and down****

y#=amplitude#*sin(theta#)
theta#=theta#+frequency#

if fixed(tile(a,b))=1
set matrix height 2,a,b,y#
endif
if fixed(tile(c,d))=1
set matrix height 2,c,d,y#
endif

`********************************************

`****make wave propogate through matrix******

for x=1+1 to rows-1
for z=1+1 to columns-1

`check x+1 height
if x<rows
distxp1#=get matrix height(2,x+1,z)-get matrix height (2,x,z)
endif


`check x-1 height
if x>1
distxm1#=get matrix height(2,x-1,z)-get matrix height (2,x,z)
endif

`check z+1 height
if z<columns
distzp1#=get matrix height(2,x,z+1)-get matrix height (2,x,z)
endif

`check z-1 height
if z>1
distzm1#=get matrix height(2,x,z-1)-get matrix height (2,x,z)
endif

`calculate vector some of heights adjacent to tile
`and make this proportional to the tile's acceleration
vectorsum#=distxp1#+distxm1#+distzp1#+distzm1#
a#=vectorsum#*elasticity#
`increase tile's velocity by it's acceleration amount
v#(tile(x,z))=v#(tile(x,z))+a#
if fixed(tile(x,z))=0
`reposition tile
set matrix height 2,x,z,get matrix height(2,x,z)+v#(tile(x,z))
v#(tile(x,z))=v#(tile(x,z))*damping#
endif

next z
next x

update matrix 2

`****************************************************

return


get_wave_slope:

for object=1 to number

   y#(object)=get ground height(2,x#(object),z#(object))
   yfront#=get ground height(2,newxvalue(x#(object),angle#,1),newzvalue(z#(object),angle#,1))
   yside#=get ground height(2,newxvalue(x#(object),angle#+90,1),newzvalue(z#(object),angle#+90,1))

   frontspeed#(object)=frontspeed#(object)+(y#(object)-yfront#)*current#
   sidespeed#(object)=sidespeed#(object)+(y#(object)-yside#)*current#

   x#(object)=newxvalue(x#(object),angle#,frontspeed#(object))
   z#(object)=newzvalue(z#(object),angle#,frontspeed#(object))

   x#(object)=newxvalue(x#(object),angle#+90,sidespeed#(object))
   z#(object)=newzvalue(z#(object),angle#+90,sidespeed#(object))

position object object,x#(object),y#(object),z#(object)

next object

return

control:

`***steering****


if mousex()>0.5*screen width() then polarity=1
if mousex()<0.5*screen width() then polarity=-1
`the square function is so that the steering is less sensitive
`in the middle and more sensitive for extreme turns
`note: squaring makes the mouse position always positive,
`hence the need for a polarity.
turn#=(((mousex()-(0.5*screen width()))^2)/10000)*polarity
if turn#>1 then turn#=1
if turn#<-1 then turn#=-1
yrotate object 1,object angle y(1)+turn#
direction#=wrapvalue(direction#+turn#)
rotate limb 1,1,0,0,-turn#*20 `roll the jetski body according to turn amount

`***************

`*****speed*****

if upkey()=1 then frontspeed#(1)=frontspeed#(1)+(power#*cos(direction#)):sidespeed#(1)=sidespeed#(1)+(power#*sin(direction#))
if downkey()=1 then frontspeed#(1)=frontspeed#(1)-(power#*cos(direction#)):sidespeed#(1)=sidespeed#(1)-(power#*sin(direction#))

for object=1 to number
frontspeed#(object)=frontspeed#(object)*friction#
sidespeed#(object)=sidespeed#(object)*friction#
next object

`***************

return

chasecam:

yrotate camera curveangle(direction#,camera angle y(),smoothness)
position camera object position x(1),object position y(1)+camheight#,object position z(1)
move camera -camdist#
point camera object position x(1),object position y(1)+camheight#,object position z(1)

return

`********have a nice day*************

IF (boathitsmountainslol) THEN BoatAngle=360-BoatAngle

Sync On:Sync Rate 60:Autocam Off
Set Camera Range .01,100
Global GSize as Integer
GSize=100
Dim GlobalHeights(GSize,GSize)
Global Grav as Float
Grav=-.01
Global stearing as float
Global move# as float
Global acc# as float
Global Speed as float
Global friction# as float
Global heavy as float
Global Choing as Integer

Null=Make Vector3(1)
Null=Make Vector3(2)
Null=Make Vector3(3)

Type Vector
   Object as DWord
   x as float
   y as float
   z as float
Endtype

Dim ObjectPhysics(3) as Vector



Print "Drawing Grass..."
Sync
Grass()
Print "Global Heights..."
Sync
Establish_Globals()

Make_Terrain(64,64)
Print "Terrain Done..."
Sync


Position Camera 0,0,0
x#=Camera Position X()
z#=Camera Position Z()
y#=Get Terrain Ground Height(1,x#,z#)+150
Position Camera x#,y#,z#

Make_Vehicle()


lx=1:ly=1
Do
   turn#=(RightKey()-LeftKey())*1.0
   stearing = stearing + turn#-((stearing>0)-(stearing<0))*.3
   move#=(DownKey()-Upkey())*.002
   acc#=acc#+move#
   heavy=0
   If Choing=0
      If Lower$(Inkey$())="f"
         Move Object Up 21,.4
         Move Object Up 23,.4
         Choing=100
      Endif
      n#=0
      If SpaceKey()
         For i = 21 to 24
            Move Object Up i,.4+n#
            n#=n#+.02
         Next i
         Choing=100
      Endif
   Else
      Choing=Choing-1
   Endif
   For i = 21 to 22
      if stearing <-30 then stearing = -30
      if stearing >30 then stearing = 30
      gy#=Get Terrain Ground Height(1,Object Position X(i),Object Position Z(i))+Object Size Y(i)*.5
      if gy#>Object Position Y(i)
         Move Object i,acc#
         friction#=.35
      Else
         friction#=.999
         heavy=heavy-.0025
         if gy#<Object Position Y(i)-1.0 Then acc#=acc#-(.001*acc#)
      Endif
      Physics(i)
      Center Text Object Screen X(i),Object Screen Y(i)-20,Str$(i)
   Next i
   for i = 23 to 24
      gy#=Get Terrain Ground Height(1,Object Position X(i),Object Position Z(i))+Object Size Y(i)*.5
      if gy#>Object Position Y(i)
         Move Object i,acc#*.5
         friction#=.45
      Else
         friction#=.998
         heavy=heavy-.0025
         if gy#<Object Position Y(i)-1.0 Then acc#=acc#-(.001*acc#)
      Endif
      Physics(i)
      Center Text Object Screen X(i),Object Screen Y(i)-20,Str$(i)
   next i
   Space_Wheels(21,22,23,24,25,26,27)
   Rotate Camera 0,0,0
   Pitch Camera Down 6
   Move Camera -18
   Set Cursor 0,0
   Print "FPS: ";Screen FPS()
   Print "Speed: ";Int(-acc#*100)
   Sync
Loop

End

Function Make_Terrain(x,y)
   Print "Heightmap..."
   Sync
   Create_Heightmap(x,y)
   Print "Make Terrain..."
   Sync
   If Object Exist(1) then Delete Object 1
   Make Object Terrain 1
   Set Terrain Heightmap 1,"Heightmap.bmp"
   Set Terrain Texture 1,3,4
   Set Terrain Tiling 1,2
   Set Terrain Scale 1,5,.3,5
   Build Terrain 1
   Set Object Collision Off 1
Endfunction


Function Establish_Globals()
   GlobalHeights(1,1)=Rnd(50)+100
   For i=2 to GSize
      GlobalHeights(i,1)=rnd(50)+100
      adj=GlobalHeights(i-1,1)+rnd(13)-7
      GlobalHeights(i,1)=adj
   Next i
   For i=2 to GSize
      GlobalHeights(1,i)=rnd(50)+100
      adj=GlobalHeights(1,i-1)+rnd(13)-7
      GlobalHeights(1,i)=adj
   Next i
   For j = 2 to GSize
      For i= 2 to GSize
         adj=GlobalHeights(i-1,j-1)
         adj=adj+GlobalHeights(i,j-1)
         adj=adj+GlobalHeights(i-1,j)
         adj=adj/3
         nh=rnd(13)-7
         gh=adj+nh
         if gh<0 then gh=0
         if gh>200 then gh=200
         GlobalHeights(i,j)=gh
      Next i
   Next j
   Sleep 1000
Endfunction

Function Global_Position(n)
   Do
      if n>0 and n<GSize then Exit
      if n<1 then n=n+GSize
      if n>GSize then n=n-GSize
   Loop
Endfunction n

Function Splotch(h,s)
   Lock Pixels
   Ink 0,0
   Box 200,200,230+s,230+s
   c=rnd(7)+h-3
   Ink Rgb(c,c,c),0
   r=rnd(s/2)+s/4
   For i# = 0.0 to 360.0 step 2
      x=sin(i#)*r+215+s/2
      y=cos(i#)*r+215+s/2
      line 207+s/2,207+s/2,x,y
      r=r+(rnd(2)-1)
      if r>s+15 then r=s+15
      m=s-15:If m<2 then m=2
      if r<m then r=m
   Next i#
   Unlock Pixels
   Get Image 1,200,200,230+s,230+s
Endfunction

Function Create_Heightmap(x,y)
   Ink rgb(128,128,128),0
   Box 0,0,199,199
   For ny=0 to 5
      For nx=0 to 5
         Randomize Timer()
         TileHeight=10+rnd(6)+rnd(6)+rnd(6)
         ox=32+(nx*32):oy=32+(ny*32)
         mx=32+((nx+1)*32):my=32+((ny+1)*32)
         Splotch(TileHeight,32)
         Sprite 1,200,200,1
         Set Sprite 1,0,1
         Set Sprite Alpha 1,200
         Paste Sprite 1,ox+rnd(32)-16,oy+rnd(32)-16
         Box ox,oy,mx,my
      Next nx
   Next ny
   For ny=0 to 5
      For nx=0 to 5
         Randomize Timer()
         ox=32+(nx*32):oy=32+(ny*32)
         For i = 1 to rnd(3)+2
            sx=ox+rnd(32)-16:sy=oy+rnd(32)-16
            h=point(sx,sy)+Rnd(5)-Rnd(5)
            Splotch(h,15)
            Sprite 1,200,200,1
            Set Sprite 1,0,1
            Set Sprite Alpha 1,128
            Paste Sprite 1,sx,sy
         Next i
         Randomize Timer()
         ox=32+(nx*32):oy=32+(ny*32)
         for noise=1 to rnd(50)+rnd(50)+rnd(50)+50
            noisex=rnd(32)+ox
            noisey=rnd(32)+oy
            c=point(noisex,noisey)
            c=rgbr(c)+(rnd(5)-2)
            if c<0 then c=0
            if c>255 then c=255
            ink rgb(c,c,c),0
            dot noisex,noisey
         Next noise
      Next nx
   Next ny
   if file exist("Heightmap.bmp") Then Delete File "Heightmap.bmp"
   Get Image 2,32,32,32+x*5,32+y*5,0
   Save Image "Heightmap.bmp",2
   If Sprite Exist(1) Then Delete Sprite 1
   Ink Rgb(255,255,255),0
   Print "Heightmap Done"
   Sync
Endfunction


Function Grass()
   If Bitmap Exist(1) Then Delete Bitmap 1
   Create Bitmap 1,512,512
   Set Current Bitmap 1
   Ink Rgb(0,90,0),0
   Box 0,0,511,511
   g=0
   Do
      Splotch(180,15)
      Sprite 1,200,200,1
      Size Sprite 1,rnd(3)+1,Rnd(3)+1
      Set Sprite Alpha 1,128
      Set Sprite Diffuse 1,0,255,0
      Paste Sprite 1,rnd(67)-2,rnd(67)-2
      g=g+1
      If rnd(500)+500<g Then Exit
   Loop
   Get Image 3,1,1,64,64,0
   Ink Rgb(0,90,0),0
   Box 0,0,511,511
   g=0
   Do
      Splotch(rnd(80)+130,30)
      Sprite 1,200,200,1
      Size Sprite 1,rnd(10)+10,Rnd(10)+10
      Set Sprite Alpha 1,64
      Set Sprite Diffuse 1,0,255,0
      Paste Sprite 1,rnd(75)-10,rnd(75)-10
      g=g+1
      If rnd(1000)+200<g Then Exit
   Loop
   Get Image 4,1,1,64,64,0
   Set Current Bitmap 0
   Delete Bitmap 1
   Delete Sprite 1
Endfunction


Function Physics(ObjectID)
   Flag=0
   For i = 0 to Array Count(ObjectPhysics())
      if ObjectPhysics(i).Object=ObjectID
         Flag=1
         Exit
      Endif
   Next i
   If Flag = 0 Then ExitFunction
   x#=Object Position X(ObjectID):y#=Object Position Y(ObjectID):z#=Object Position Z(ObjectID)
   Position Object ObjectID,x#+X Vector3(ObjectID),y#+Y Vector3(ObjectID)+Grav+heavy,z#+Z Vector3(ObjectID)
   If Object Position X(ObjectID)>Object Size X(1)-1
      Position Object ObjectID,Object Size X(1)-1,Object Position Y(ObjectID),Object Position Z(ObjectID)
   Endif
   If Object Position X(ObjectID)<1
      Position Object ObjectID,1,Object Position Y(ObjectID),Object POsition Z(ObjectID)
   Endif
   If Object Position Z(ObjectID)>Object Size Z(1)-1
      Position Object ObjectID,Object Position X(ObjectID),Object Position Y(ObjectID),Object Size Z(1)-1
   Endif
   If Object Position Z(ObjectID)<1
      Position Object ObjectID,Object Position X(ObjectID),Object Position Y(ObjectID),1
   Endif
   gy#=Get Terrain Ground Height(1,Object Position X(ObjectID),Object Position Z(ObjectID))+Object Size Y(ObjectID)*.5
   dy#=gy#-Object Position Y(ObjectID)
   bounce=0
   if dy#>0
      Position Object ObjectID,Object Position X(ObjectID),Object Position Y(ObjectID)+dy#*.85,Object Position Z(ObjectID)
      bounce=1
   Endif

   Set Vector3 ObjectID,Object Position X(ObjectID)-ObjectPhysics(i).x,Object Position Y(ObjectID)-ObjectPhysics(i).y,Object Position Z(ObjectID)-ObjectPhysics(i).z
   Speed=Length Vector3(ObjectID)
   If Speed>1
      Multiply Vector3 ObjectID,1.0/Length Vector3(ObjectID)
   Endif
   If Bounce=1
      Set Vector3 ObjectID,X Vector3(ObjectID),Y Vector3(ObjectID)*.1,Z Vector3(ObjectID)
   Endif
   Set Vector3 ObjectID,X Vector3(ObjectID)*friction#,Y Vector3(ObjectID)*friction#,Z Vector3(ObjectID)*friction#
   ObjectPhysics(i).x=Object Position X(ObjectID)
   ObjectPhysics(i).y=Object Position Y(ObjectID)
   ObjectPhysics(i).z=Object Position Z(ObjectID)
Endfunction

Function Space_Wheels(Wheel1,Wheel2,Wheel3,Wheel4,Axel1,Axel2,Axel3)
   For run = 1 to 2
      ax#=(Object Position X(Wheel1)+Object Position X(Wheel2)+Object Position X(Wheel3)+Object Position X(Wheel4))/4.0
      ay#=(Object Position Y(Wheel1)+Object Position Y(Wheel2)+Object Position Y(Wheel3)+Object Position Y(Wheel4))/4.0
      az#=(Object Position Z(Wheel1)+Object Position Z(Wheel2)+Object Position Z(Wheel3)+Object Position Z(Wheel4))/4.0

      lay#=ay#:low=0
      If Object Position Y(Wheel1)<lay#
         lay#=Object Position Y(Wheel1)
         low=Wheel1
      Endif
      If Object Position Y(Wheel2)<lay#
         lay#=Object Position Y(Wheel2)
         low=Wheel2
      Endif
      If Object Position Y(Wheel3)<lay#
         lay#=Object Position Y(Wheel3)
         low=Wheel3
      Endif
      If Object Position Y(Wheel4)<lay#
         lay#=Object Position Y(Wheel4)
         low=Wheel4
      Endif
      If ay#-lay#>.5 and low<>0
         Set Vector3 1,Object Position X(low),lay#,Object Position Z(low)
         Normalize Vector3 1,1
         ax#=ax#+X Vector3(1)*.005:az#=az#+Z Vector3(1)*.005
      Endif

      fcx#=(Object Position X(Wheel1)+Object Position X(Wheel2))/2.0
      fcy#=(Object Position Y(Wheel1)+Object Position Y(Wheel2))/2.0
      fcz#=(Object Position Z(Wheel1)+Object Position Z(Wheel2))/2.0

      Set Vector3 1,Object Position X(Wheel2)-Object Position X(Wheel1),Object Position Y(Wheel2)-Object Position Y(Wheel1),Object Position Z(Wheel2)-Object Position Z(Wheel1)
      Position Object Axel1,fcx#,fcy#,fcz#
      Point Object Axel1,fcx#+X Vector3(1),fcy#+Y Vector3(1),fcz#+Z Vector3(1)
      Position Object Wheel1,Limb Position X(Axel1,2),Limb Position Y(Axel1,2),Limb Position Z(Axel1,2)
      Position Object Wheel2,Limb Position X(Axel1,1),Limb Position Y(Axel1,1),Limb Position Z(Axel1,1)

      rcx#=(Object Position X(Wheel3)+Object Position X(Wheel4))/2.0
      rcy#=(Object Position Y(Wheel3)+Object Position Y(Wheel4))/2.0
      rcz#=(Object Position Z(Wheel3)+Object Position Z(Wheel4))/2.0

      Set Vector3 1,Object Position X(Wheel4)-Object Position X(Wheel3),Object Position Y(Wheel4)-Object Position Y(Wheel3),Object Position Z(Wheel4)-Object Position Z(Wheel3)
      Position Object Axel1,rcx#,rcy#,rcz#
      Point Object Axel1,rcx#+X Vector3(1),rcy#+Y Vector3(1),rcz#+Z Vector3(1)
      Position Object Wheel3,Limb Position X(Axel1,2),Limb Position Y(Axel1,2),Limb Position Z(Axel1,2)
      Position Object Wheel4,Limb Position X(Axel1,1),Limb Position Y(Axel1,1),Limb Position Z(Axel1,1)

      Set Vector3 1,Object Position X(Wheel4)-Object Position X(Wheel1),Object Position Y(Wheel4)-Object Position Y(Wheel1),Object Position Z(Wheel4)-Object Position Z(Wheel1)
      Position Object Axel2,ax#,ay#,az#
      Point Object Axel2,ax#+X Vector3(1),ay#+Y Vector3(1),az#+Z Vector3(1)
      Position Object Wheel1,Limb Position X(Axel2,2),Limb Position Y(Axel2,2),Limb Position Z(Axel2,2)
      Position Object Wheel4,Limb Position X(Axel2,1),Limb Position Y(Axel2,1),Limb Position Z(Axel2,1)

      Set Vector3 1,Object Position X(Wheel3)-Object Position X(Wheel2),Object Position Y(Wheel3)-Object Position Y(Wheel2),Object Position Z(Wheel3)-Object Position Z(Wheel2)
      Position Object Axel3,ax#,ay#,az#
      Point Object Axel3,ax#+X Vector3(1),ay#+Y Vector3(1),az#+Z Vector3(1)
      Position Object Wheel2,Limb Position X(Axel3,2),Limb Position Y(Axel3,2),Limb Position Z(Axel3,2)
      Position Object Wheel3,Limb Position X(Axel3,1),Limb Position Y(Axel3,1),Limb Position Z(Axel3,1)
   Next run

   Point Object Wheel1,fcx#,fcy#,fcz#
   Turn Object Left Wheel1,90.0-stearing

   Point Object Wheel2,fcx#,fcy#,fcz#
   Turn Object Right Wheel2,90.0+stearing


   Point Object Wheel3,Object Position X(Wheel4),Object Position Y(Wheel4),Object Position Z(Wheel4)
   Turn Object Left Wheel3,90.0
   Point Object Wheel4,Object Position X(Wheel3),Object Position Y(Wheel3),Object Position Z(Wheel3)
   Turn Object Right Wheel4,90.0

   Position Camera ax#,ay#,az#
Endfunction

Function Make_Axel(ObjectID,length#)
   Make Object Triangle ObjectID, 0,0,0, 0,.01,0, 0,0,0
   Make Mesh From Object 1,ObjectID
   Delete Object ObjectID
   Make Object Cylinder ObjectID,.1
   Scale Object ObjectID,100,length#*1000.0,100
   Pitch Object Down ObjectID,90
   Fix Object Pivot ObjectID
   Make Mesh From Object 2,ObjectID
   Delete Object ObjectID
   Make Object ObjectID,2,0
   Add Limb ObjectID,1,1
   Offset Limb ObjectID,1,0,0,length#*.5
   Add Limb ObjectID,2,1
   Offset Limb ObjectID,2,0,0,length#*-.5
Endfunction


Function Make_Vehicle()
   For i = 21 to 24
      Make Object Cylinder i,1
      Roll Object Right i,90
      Fix Object Pivot i
      Set Object Cull i,0
      Null=Make Vector3(i)
      ObjectPhysics(i-21).Object = i
   Next i
   Position Object 21,5,20,5
   Position Object 22,2,20,5
   Position Object 23,5,20,1
   Position Object 24,2,20,1
   Make_Axel(25,3)
   Make_Axel(26,5)
   Make_Axel(27,5)
Endfunction

`darkbasicproguys racing entry
` all the setup stuff

Sync on : Sync rate 60
Set display mode 1024,768,16 : Hide mouse
autocam off `Set camera range 1,0x7fffffff

rem start


mtrxLand as integer : mtrxLand=1
txtrGrass as integer : txtrGrass=2
objCar as integer : objCar=1
hillheight=100

` make grass texture (from zotoasters entry)

create bitmap 1,128,128
set current bitmap 1
cls rgb(0,128,0)
for d=1 To 10000
ink rgb(0,rnd(255),rnd(50)),0
dot rnd(128),rnd(128)
next d
blur bitmap 1,4
get image txtrGrass,0,0,128,128
set current bitmap 0
delete bitmap 1

` make land (from zotoasters entry)

make matrix mtrxLand,20000,20000,50,50
randomize matrix mtrxLand,hillheight
prepare matrix texture mtrxLand,txtrGrass,1,1
update matrix mtrxLand
position matrix mtrxLand,0,0,0

` make car

make object box objCar,.5,.3,.8
color object objCar,rgb(2,255,255)
position object objCar,0,1,30

` camera

position camera 0,5,0

` main program loop*****

do
   control camera using arrowkeys 0,.1,1
   position camera camera position x(),5,camera position z()
   sync
loop

`Standard setup
Sync on : Sync rate 60
Set display mode 1024,768,16 : Hide mouse
Autocam off : Set camera range 1,0x7fffffff

`Load Physics Functions
`#include "ZPhysics.dba"

`Start World Physics
Start_World_Physics()
Set_World_Standards()

`Start Object Physics
wheels=4
Start_Object_Physics(1)
For w=1 to wheels+1
Start_Object_Physics(w)
   Set_Object_Constants(w,20,0.97,0.98,10.0)
   Weigh_Object(w)
Next w
turn#=25.0
fuel#=200.0

`Create Player
for w=1 to wheels
   make object box w,10,30,30
next w
make object box 5,50,30,80 : color object 5,rgb(255,0,0)

`Create world
create bitmap 1,128,128 : set current bitmap 1
cls rgb(0,128,0)
for d=1 To 10000
ink rgb(0,rnd(255),rnd(50)),0
dot rnd(128),rnd(128)
next d
blur bitmap 1,4
get image 10,0,0,128,128
set current bitmap 0 : delete bitmap 1

`Create and Mold
make matrix 1,20000,20000,50,50
randomize matrix 1,50
prepare matrix texture 1,10,1,1
update matrix 1

`**Main loop**
Do

`Controls
if upkey()=1 or keystate(17)=1 then set_linear_force(5,dirX(5)*MoveConstant(5),dirY(5),DirZ(5)*MoveConstant(5))
if leftkey()=1
   set_angular_force(1,0,-(Movespeed(1)/turn#),0)
   set_angular_force(2,0,-(Movespeed(1)/turn#),0)
endif
if rightkey()=1
   set_angular_force(1,0,(Movespeed(1)/turn#),0)
   set_angular_force(2,0,(Movespeed(1)/turn#),0)
endif
if upkey()=0 and keystate(17)=0 or fuel#=0.0 then set_linear_force(5,0,0,0)
if leftkey()=0 and rightkey()=0 then set_angular_force(1,0,0,0) : set_angular_force(2,0,0,0)
Set_Next_Angle_Y(5,1,10.0)
For w=3 to 4
   Set_Next_Angle_Y(w,5,10.0)
Next w
for w=1 to 4
   AngleX(w)=Wrapvalue(AngleX(w)+(movespeed(5)/3.14))
next w

`Update Physics
for w=1 to wheels
   Set_Linear_Velocity(w,VelX(5),VelY(5),VelZ(5))
   Position_At_Object(w,5)
   Place_On_Matrix(w,1,Object Size Y(w)/2)
next w
Place_On_Matrix(5,1,(Object Size Y(5)/2)+(Object Size Y(5)/5))
Calculate_Offsets_AB(1,5,0,315,60) : Calculate_Offsets_AB(2,5,0,45,60)
Calculate_Offsets_AB(3,5,0,135,60) : Calculate_Offsets_AB(4,5,0,215,60)
for w=1 to 5
Get_Old_Position(w)
Get_Object_Direction(w)
Get_Traction_Force(w)
Get_Drag_Force(w)
Get_Force(w)
Get_Speed(w)
Get_Acceleration(w)
Set_Acceleration(w)
Set_Velocity(w)
Get_Angular_Acceleration(w)
Get_Angular_Velocity(w)
Update_Coordinates_Slide(w)
Update_Angles(w)
Get_Lateral_Forces(w)
Fully_Update_Object_Position(w)
Update_Object_Rotation(w)
Get_Movement_Angle(w)
next w
if upkey()=0 and keystate(17)=0
   if abs(MoveAngY(1)-AngleY(1)) < abs((MoveAngY(1)+180)-AngleY(1))
      angleY(1)=curveangle(MoveAngY(1),angleY(1),10.0)
      angleY(2)=curveangle(MoveAngY(1),angleY(2),10.0)
   else
      angleY(1)=curveangle(MoveAngY(1)+180,angleY(1),10.0)
      angleY(2)=curveangle(MoveAngY(1)+180,angleY(2),10.0)
   endif
endif


`Camera
ca#=curveangle(angleY(5),ca#,20.0)
cx#=positionX(5)-sin(ca#)*600
cz#=positionZ(5)-cos(ca#)*600
cy#=positionY(5)+300
position camera cx#,cy#,cz#
point camera positionX(5),positionY(5),positionZ(5)

`Create DashBoard
dec fuel#,movespeed(5)/2000
if fuel#<=0.0 then fuel#=0.0
ink rgb(0,0,0),0 : box 0,648,1024,768
create_meter(movespeed(5),60,708,5,30,rgb(255,0,255),rgb(255,255,0),"0","90","180","270")
create_meter(fuel#,964,708,5,30,rgb(255,0,255),rgb(255,255,0),"0","90","180","270")
ink rgb(255,255,255),0
text 110,652,"Speed: "
text 914-text width("Fuel: "),652,"Fuel: "

`**End loop**
Sync
Loop


`**Functions**
Function Start_World_Physics()
Dim Gravity(0) as float
Dim Drag(0) as Float
Dim Windx(0) as float
Dim Windy(0) as float
Dim Windz(0) as float
Endfunction

Function Start_Object_Physics(o)
Dim Positionx(o) as Float : Dim Positiony(o) as Float : Dim Positionz(o) as Float
Dim OldPosX(o) as Float : Dim OldPosY(o) as Float : Dim OldPosZ(o) as Float
Dim Anglex(o) as Float : Dim Angley(o) as Float : Dim Anglez(o) as Float
Dim MoveAngX(o) as Float : Dim MoveAngY(o) as Float : Dim MoveAngZ(o) as Float
Dim LongForceX(o) as Float : Dim LongForceY(o) as Float : Dim LongForceZ(o) as Float
Dim TracForceX(o) as Float : Dim TracForceY(o) as Float : Dim TracForceZ(o) as Float
Dim DragForceX(o) as Float : Dim DragForceY(o) as Float : Dim DragForceZ(o) as Float
Dim AngForceX(o) as Float : Dim AngForceY(o) as Float : Dim AngForceZ(o) as Float
Dim AngVelX(o) as Float : Dim AngVelY(o) as Float : Dim AngVelZ(o) as Float
Dim AngAccX(o) as Float : Dim AngAccY(o) as Float : Dim AngAccZ(o) as Float
Dim VelX(o) as Float : Dim VelY(o) as Float : Dim VelZ(o) as Float
Dim AccX(o) as Float : Dim AccY(o) as Float : Dim AccZ(o) as Float
Dim DirX(o) as Float : Dim DirY(o) as Float : Dim DirZ(o) as Float
Dim OffsetX(o) as Float : Dim OffsetY(o) as Float : Dim OffsetZ(o) as Float
Dim MoveTime(o) as Float : Dim MoveDist(o) as Float : Dim Power(o) as Float
Dim Movespeed(o) as Float : Dim MoveConstant(o) as Float : Dim Acceleration(o) as Float
Dim TractionConstant(o) as Float : Dim AngularTractionConstant(o) as Float
Dim Mass(o) as Float : Dim Weight(o) as Float
Dim KineticEnergy(o) as Float: Dim WorkDone(o) as Float
Dim LatForce(o) as Float : Dim PushX(o) as Float : Dim PushY(o) as Float : Dim PushZ(o) as Float
Dim Force(o) as Float
Endfunction

Function Set_Standard_Gravity()
Gravity(0)=9.8
Endfunction

Function Set_Standard_Drag()
Drag(0)=0.43
Endfunction

Function Set_World_Standards()
Gravity(0)=9.8
Drag(0)=0.43
Endfunction

Function Set_Gravity(g#)
Gravity(0)=g#
Endfunction

Function Set_Drag(d#)
Drag(0)=d#
Endfunction

Function Set_Object_Constants(o,mv#,t#,a#,m#)
Set_Move_Constant(o,mv#)
Set_Traction_Constant(o,t#)
Set_Angular_Traction_Constant(o,a#)
Set_Body_Mass(o,m#)
Endfunction

Function Set_Move_Constant(o,c#)
MoveConstant(o)=c#
Endfunction

Function Set_Traction_Constant(o,c#)
TractionConstant(o)=c#
Endfunction

Function Set_Angular_Traction_Constant(o,c#)
AngularTractionConstant(o)=c#
Endfunction

Function Set_Body_Mass(o,m#)
Mass(o)=m#
Endfunction

Function Weigh_Object(o)
Weight(o)=Mass(o)*Gravity(0)
Endfunction

Function Set_Linear_Force(o,x,y,z)
LongForceX(o)=x : LongForceY(o)=y : LongForceZ(o)=z
Endfunction

Function Set_Linear_Velocity(o,x,y,z)
VelX(o)=x : VelY(o)=y : VelZ(o)=z
Endfunction

Function Set_Angular_Force(o,x,y,z)
AngForceX(o)=x : AngForceY(o)=y : AngForceZ(o)=z
Endfunction

Function Set_Angular_Velocity(o,x,y,z)
AngVelX(o)=x : AngVelY(o)=y : AngVelZ(o)=z
Endfunction

Function Move_Object(o,spd)
VelX(o)=dirX(o)*spd
VelY(o)=dirY(o)*spd
VelZ(o)=dirZ(o)*spd
Endfunction

Function Push_Object(o,spd)
LongForceX(o)=dirX(o)*spd
LongForceY(o)=dirY(o)*spd
LongForceZ(o)=dirZ(o)*spd
Endfunction

Function Control_Object_Using_Forces(o,push,turn)
if upkey()=1 then Push_Object(o,push)
if downkey()=1 then Push_Object(o,-push)
if leftkey()=1 then set_angular_Force(o,0,-(MoveSpeed(o)/turn),0)
if rightkey()=1 then set_angular_Force(o,0,(MoveSpeed(o)/turn),0)
if upkey()=0 and downkey()=0 then Push_Object(o,0)
if leftkey()=0 and rightkey()=0 then set_angular_force(o,0,0,0)
Standard_Update(o)
Endfunction

Function Control_Object_Using_Velocities(o,spd,turn)
if upkey()=1 then Move_Object(o,spd)
if downkey()=1 then Move_Object(o,-spd)
if leftkey()=1 then set_angular_velocity(o,0,-turn,0)
if rightkey()=1 then set_angular_velocity(o,0,turn,0)
if upkey()=0 and downkey()=0 then Move_Object(o,0)
if leftkey()=0 and rightkey()=0 then set_angular_velocity(o,0,0,0)
Standard_Update(o)
Endfunction

Function Get_Movement_Angle(o)
MoveAngX(o)=AtanFull(VelY(o),VelZ(o))
MoveAngY(o)=AtanFull(VelX(o),VelZ(o))
MoveAngZ(o)=AtanFull(VelX(o),VelY(o))
Endfunction

Function Set_Next_Angle_X(a,b,smth#)
AngleX(a)=Curveangle(AngleX(b),AngleX(a),smth#)
Endfunction

Function Set_Next_Angle_Y(a,b,smth#)
AngleY(a)=Curveangle(AngleY(b),AngleY(a),smth#)
Endfunction

Function Set_Next_Angle_Z(a,b,smth#)
AngleZ(a)=Curveangle(AngleZ(b),AngleZ(a),smth#)
Endfunction

Function Position_At_Object(a,b)
PositionX(a)=PositionX(b)
PositionY(a)=PositionY(b)
PositionZ(a)=PositionZ(b)
Endfunction

Function Get_Old_Position(o)
OldPosX(o)=PositionX(o)-VelX(o)
OldPosY(o)=PositionY(o)-VelY(o)
OldPosZ(o)=PositionZ(o)-VelZ(o)
Endfunction

Function Get_Object_Direction(o)
DirX(o)=sin(AngleY(o))*cos(AngleX(o))
DirY(o)=-sin(AngleX(o))
DirZ(o)=cos(AngleY(o))*cos(AngleX(o))
Endfunction

Function Get_Traction_Force(o)
TracForceX(o)=DirX(o)*MoveConstant(o)
TracForceY(o)=DirY(o)*MoveConstant(o)
TracForceZ(o)=DirZ(o)*MoveConstant(o)
Endfunction

Function Get_Drag_Force(o)
DragForceX(o)=-Drag(0)*VelX(o)*MoveSpeed(o)
DragForceY(o)=-Drag(0)*VelY(o)*MoveSpeed(o)
DragForceZ(o)=-Drag(0)*VelZ(o)*MoveSpeed(o)
Endfunction

Function Get_Force(o)
Force(o)=Sqrt(LongForceX(o)^2+LongForceY(o)^2+LongForceZ(o)^2)
Endfunction

Function Get_Speed(o)
MoveSpeed(o)=sqrt(VelX(o)^2 + VelY(o)^2 + VelZ(o)^2)
Endfunction

Function Get_Acceleration(o)
Acceleration(o)=Sqrt(AccX(o)^2+AccY(o)^2+AccZ(o)^2)
Endfunction

Function Set_Acceleration(o)
AccX(o)=LongForceX(o)/Mass(o)
AccY(o)=LongForceY(o)/Mass(o)
AccZ(o)=LongForceZ(o)/Mass(o)
Endfunction

Function Set_Velocity(o)
VelX(o)=VelX(o)+AccX(o) : VelX(o)=VelX(o)*TractionConstant(o)
VelY(o)=VelY(o)+AccY(o) : VelY(o)=VelY(o)*TractionConstant(o)
VelZ(o)=VelZ(o)+AccZ(o) : VelZ(o)=VelZ(o)*TractionConstant(o)
Endfunction

Function Get_Angular_Acceleration(o)
AngAccX(o)=AngForceX(o)/Mass(o)
AngAccY(o)=AngForceY(o)/Mass(o)
AngAccZ(o)=AngForceZ(o)/Mass(o)
Endfunction

Function Get_Angular_Velocity(o)
AngVelX(o)=AngVelX(o)+AngAccX(o) : AngVelX(o)=AngVelX(o)*AngularTractionConstant(o)
AngVelY(o)=AngVelY(o)+AngAccY(o) : AngVelY(o)=AngVelY(o)*AngularTractionConstant(o)
AngVelZ(o)=AngVelZ(o)+AngAccZ(o) : AngVelZ(o)=AngVelZ(o)*AngularTractionConstant(o)
Endfunction

Function Update_Coordinates_Straight(o)
PositionX(o)=PositionX(o)+(DirX(o)*Movespeed(o))
PositionY(o)=PositionY(o)+(DirY(o)*Movespeed(o))
PositionZ(o)=PositionZ(o)+(DirZ(o)*Movespeed(o))
Endfunction

Function Update_Coordinates_Slide(o)
PositionX(o)=PositionX(o)+VelX(o)
PositionY(o)=PositionY(o)+VelY(o)
PositionZ(o)=PositionZ(o)+VelZ(o)
Endfunction

Function Update_Angles(o)
AngleX(o)=Wrapvalue(AngleX(o)+AngVelX(o))
AngleY(o)=Wrapvalue(AngleY(o)+AngVelY(o))
AngleZ(o)=Wrapvalue(AngleZ(o)+AngVelZ(o))
Endfunction

Function Get_Lateral_Forces(o)
LatForce(o)=abs(AngVelY(o))*(MoveSpeed(o)/Mass(o))
If AngVelY(o)>0.0
   PushX(o)=(sin(Wrapvalue(AngleY(o)-90))*LatForce(o))
   PushZ(o)=(cos(Wrapvalue(AngleY(o)-90))*LatForce(o))
endif
If AngVelY(o)<0.0
   PushX(o)=(sin(Wrapvalue(AngleY(o)+90))*LatForce(o))
   PushZ(o)=(cos(Wrapvalue(AngleY(o)+90))*LatForce(o))
endif
Endfunction

Function Calculate_Offsets(o,sax#,say#,out#)
OffsetX(o)=(Sin(Wrapvalue(AngleY(o)+say#))*Cos(Wrapvalue(AngleX(o)+sax#)))*Out#
OffsetY(o)=-Sin(Wrapvalue(AngleX(o)+sax#))*out#
OffsetZ(o)=(Cos(Wrapvalue(AngleY(o)+say#))*Cos(Wrapvalue(AngleX(o)+sax#)))*out#
Endfunction

Function Calculate_Offsets_AB(a,b,sax#,say#,out#)
OffsetX(a)=(Sin(Wrapvalue(AngleY(b)+say#))*Cos(Wrapvalue(AngleX(b)+sax#)))*out#
OffsetY(a)=-Sin(Wrapvalue(AngleX(b)+sax#))*out#
OffsetZ(a)=(Cos(Wrapvalue(AngleY(b)+say#))*Cos(Wrapvalue(AngleX(b)+sax#)))*out#
Endfunction

Function Place_On_Matrix(o,m,h#)
PositionY(o)=Get ground height(m,positionX(o),PositionZ(o))+h#
Endfunction

Function Update_Object_Position(o)
Position Object o,PositionX(o),PositionY(o),PositionZ(o)
Endfunction

Function Update_Position_With_Push_Forces(o,a,b,c)
Position object o,Positionx(o)+(PushX(o)*a),Positiony(o)+(PushY(o)*b),Positionz(o)+(PushZ(o)*c)
Endfunction

Function Update_Position_With_Offsets(o,a,b,c)
Position Object o,Positionx(o)+(offsetx(o)*a),Positiony(o)+(offsety(o)*b),PositionZ(o)+(offsetZ(o)*c)
Endfunction

Function Fully_Update_object_Position(o)
Position Object o,PositionX(o)+PushX(o)+OffsetX(o),PositionY(o)+PushY(o)+OffsetY(o),PositionZ(o)+PushZ(o)+OffsetZ(o)
Endfunction

Function Update_object_Rotation(o)
Rotate object o,AngleX(o),AngleY(o),AngleZ(o)
Endfunction

Function Standard_Update(o)
Get_Old_Position(o)
Get_Object_Direction(o)
Get_Traction_Force(o)
Get_Drag_Force(o)
Get_Force(o)
Get_Speed(o)
Get_Acceleration(o)
Set_Acceleration(o)
Set_Velocity(o)
Get_Angular_Acceleration(o)
Get_Angular_Velocity(o)
Update_Coordinates_Straight(o)
Update_Angles(o)
Get_Lateral_Forces(o)
Fully_Update_Object_Position(o)
Update_Object_Rotation(o)
Endfunction

Function Update_Speed(o)
if MoveSpeed(o)<0.9 then MoveSpeed(o)=MoveSpeed(o)*0.0
Endfunction

Function Get_Movement_Time(o)
If MoveSpeed(o)<>0 then MoveTime(o)=Movetime(o)+1.0 else MoveTime(o)=0.0
Endfunction

Function Get_Movement_Distance(o)
MoveDist(o)=MoveSpeed(o)*MoveTime(o)
Endfunction

Function Get_Kinetic_Energy(o)
KineticEnergy(o)=(Mass(o)/2)*(Movespeed(o)^2)
Endfunction

Function Get_Work_Done(o)
WorkDone(o)=Force(o)*MoveDist(o)
Endfunction

Function Get_Excerted_Power(o)
If MoveTime(o)>0.0 then Power(o)=KineticEnergy(o)/MoveTime(o) else Power(o)=0.0
Endfunction

Function create_meter(var#,speedox,speedoy,speedow,speedl,col1,col2,m1$,m2$,m3$,m4$)
ink col1,0
circle speedox,speedoy,speedl
line speedox,speedoy,speedox+sin(wrapvalue(int(360-var#)+180))*(speedow/2),speedoy+cos(wrapvalue(int(360-var#)+180))*(speedow/2)
line speedox,speedoy,speedox-sin(wrapvalue(int(360-var#)+180))*(speedow/2),speedoy-cos(wrapvalue(int(360-var#)+180))*(speedow/2)
line speedox+sin(wrapvalue(int(360-var#)+180))*(speedow/2),speedoy+cos(wrapvalue(int(360-var#)+180))*(speedow/2),speedox+sin(wrapvalue(int(360-var#)+270))*speedl,speedoy+cos(wrapvalue(int(360-var#)+270))*speedl
line speedox-sin(wrapvalue(int(360-var#)+180))*(speedow/2),speedoy-cos(wrapvalue(int(360-var#)+180))*(speedow/2),speedox+sin(wrapvalue(int(360-var#)+270))*speedl,speedoy+cos(wrapvalue(int(360-var#)+270))*speedl
ink col2,0
center text speedox-(speedl*1.5),speedoy,m1$ : center text speedox,speedoy-(speedl*1.5),m2$
center text speedox+(speedl*1.5),speedoy,m3$ : center text speedox,speedoy+(speedl*1.5),m4$
endfunction

`Space Race
sync rate 20:hide mouse:color backdrop 0:set display mode 1024,768,32
`fog on:fog distance 5000

`make the head
make object sphere 2,6
color object 2,RGB(132,168,253)
set object ambience 2,rgb(0,0,255)

`gooey death graphics
CREATE BITMAP 1,128,128
SET CURRENT BITMAP 1
FOR I=1 TO 64
   FOR J=1 TO 64
      IF RND(1)=1
         INK RGB(0,0,0),0
      ELSE
         INK RGB(0,0,200),0
      ENDIF
      BOX ((J*2)-2),((I*2)-2),(J*2),(I*2)
   NEXT J
NEXT I
BLUR BITMAP 1,6
GET IMAGE 2,1,1,128,128
SET CURRENT BITMAP 0
DELETE BITMAP 1

LAVA=2
MAKE MATRIX LAVA,200,5000,50,50
PREPARE MATRIX TEXTURE LAVA,2,1,1
SET MATRIX TEXTURE LAVA,1,1
SET MATRIX LAVA,0,0,1,1,1,1,1

DIM LAVASIN(50,50)

FOR I=1 TO 50
   FOR J=1 TO 50
      LAVASIN(I,J)=RND(360)
   NEXT J
NEXT I

`Create all variables
t=3
size=5
color=1
partcolor=1
bodylength=7
speed=0
nextob=87

`The body part making loop
for t=3 to bodylength

   `Make a body part
   make object sphere t,size

   `Color it
   color object t,RGB((20*color)+150,color*50,color*50)
   set object ambience t,RGB((20*color)+150,color*50,color*50)
   dec size:inc color
next t

`Make the walls for the track
make object box 10,1,200,25000
color object 10,RGB(255,0,0):set object ambience 10,RGB(255,0,0)
rotate object 10,0,180,0
make object box 11,1,200,25000
color object 11,RGB(255,0,0):set object ambience 11,RGB(255,0,0)
rotate object 11,0,180,0:position object 11,object position x(10)-150,0,0
make object box 12,200,1,25000
position object 12,object position x(10)-75,object position y(10)+100,0
color object 12,RGB(255,255,255):set object ambience 12,RGB(255,255,255)

position object 2,object position x(10)-75,0,0
POSITION MATRIX LAVA,object position x(10)-150,-100,0

`Create all of the obstacles
For x=15 to 215
   make object box x,rnd(50)+50,rnd(100)+50,rnd(100)+50
   color object x,RGB(50,200,50):set object ambience x,RGB(50,200,50)  `ghost object on x
   position object x,object position x(10)-rnd(150),object position y(10)-(rnd(200)-100),nextob
   inc nextob,150
Next x

` game loop
do
   distance2=int(object position z(2))
   `Control the way the head faces by moving the mouse
   yrotate object 2,object angle y(2)+ (mousemovex()/2)
   xrotate object 2,object angle x(2)+ (mousemovey()/2)
   if object angle x(2)<-90
      xrotate object 2,object angle x(2)-object angle x(2)
      xrotate object 2,-90
   endif
   if object angle x(2)>90
      xrotate object 2,object angle x(2)-object angle x(2)
      xrotate object 2,90
   endif

   `speed controls
   If speed < 10
      If upkey() = 1 then Inc speed
   Endif
   If speed > 0
      If downkey() = 1 then Dec speed
   Endif

   ` position the bits
   t=bodylength
   for a=4 to bodylength
      position object t,object position x(t-1),object position y(t-1),object position z(t-1)
      t=t-1
   next a

   `this want's to be outside the loop because the last body part has to be attatched to the head.
   position object 3, object position x(2),object position y(2),object position z(2)

   `move the head in the direction it's facing
   move object 2,speed

   `put the camera in the right place.
   set camera to follow object position x(2),object position y(2),object position z(2),object angle y(2),40,object position y(2),10,1
   point camera object position x(2),object position y(2),object position z(2)

   `Gooey death
   FOR I=1 TO 50
      FOR J=1 TO 50
         INC LAVASIN(I,J),2
         IF LAVASIN(I,J)>360 THEN LAVASIN(I,J)=0
         SET MATRIX HEIGHT LAVA,I,J,SIN(LAVASIN(I,J))*20
      NEXT J
   NEXT I
   UPDATE MATRIX LAVA

   `collision junk
   for a=10 to 12
      If object collision(2,a) then goto OUCH
   next a
   for b=15 to 215
      if object collision(2,b) then goto OUCH
   next b
   IF OBJECT POSITION Y(2)<=MATRIX POSITION Y(LAVA) then goto OUCH

   if camera position y() > object position y(12)
      hide object 12
   else
      show object 12
   endif

   ghost matrix on 2

   if object position z(2) > matrix position z(2)+2000
      POSITION MATRIX 2,object position x(10)-150,-100,matrix position z(2)+2000
   endif

   if object position z(2) > 20000 then goto HAPPYENDING

   distance=int(object position z(2))
   inc score,(distance-distance2)*speed

   `Text stuff
   set text size 20:ink RGB(255,255,255),0
   Text 10,10,"FPS      "+str$(screen fps())
   Text 10,22,"Speed    "+str$(speed)
   Text 10,34,"Distance "+str$(distance)
   Text 10,46,"Score    "+str$(score)
loop

OUCH:
do
ink RGB(255,255,255),0:set text size 40
center text 512,324,"YOU CRASHED AND BURNED"
loop

HAPPYENDING:
do
ink RGB(255,255,255),0:set text size 40
center text 512,324,"YOU MADE IT TO THE END.... WOO.... !!"
loop

`Standard setup
Sync on : Sync rate 60
Set display mode 1024,768,16 : Hide mouse
Autocam off : Set camera range 1,0x7fffffff

`Load Physics Functions
`#include "ZPhysics.dba"

`Start World Physics
Start_World_Physics()
Set_World_Standards()

`Create Player
create Bitmap 1,128,128 : set current bitmap 1
for x=1 to 64 : for y=1 to 64
   if rnd(1)=1 then ink rgb(133,162,119),0 else ink rgb(0,128,0),0
   box ((x*2)-2),((y*2)-2),(x*2),(y*2)
next x : next y
blur bitmap 1,3 : get image 5,1,1,128,128
set current bitmap 0 : delete bitmap 1
wheels=4
for w=1 to wheels
    create_tyre(w,10,30,20.0,0.97,0.98,10.0)
next w
Create_Chassi(5,50,30,80,20.0,0.97,0.98,10.0)
texture object 5,5 : scale object texture 5,0.1,0.1
turn#=25.0
fuel#=200.0

`Create world
create bitmap 1,128,128 : set current bitmap 1
cls rgb(0,128,0)
for d=1 To 10000
ink rgb(0,rnd(255),rnd(50)),0
dot rnd(128),rnd(128)
next d
blur bitmap 1,4
get image 10,0,0,128,128
set current bitmap 0 : delete bitmap 1

`Create and Mold
make matrix 1,20000,20000,50,50
randomize matrix 1,50
prepare matrix texture 1,10,1,1
update matrix 1

`**Main loop**
Do

`Controls
if upkey()=1 or keystate(17)=1 then set_linear_force(5,dirX(5)*MoveConstant(5),dirY(5),DirZ(5)*MoveConstant(5))
if leftkey()=1
   set_angular_force(1,0,-(Movespeed(1)/turn#),0)
   set_angular_force(2,0,-(Movespeed(1)/turn#),0)
endif
if rightkey()=1
   set_angular_force(1,0,(Movespeed(1)/turn#),0)
   set_angular_force(2,0,(Movespeed(1)/turn#),0)
endif
if upkey()=0 and keystate(17)=0 or fuel#=0.0 then set_linear_force(5,0,0,0)
if leftkey()=0 and rightkey()=0 then set_angular_force(1,0,0,0) : set_angular_force(2,0,0,0)
Set_Next_Angle_Y(5,1,10.0)
For w=3 to 4
   Set_Next_Angle_Y(w,5,10.0)
Next w

`Update Physics
for w=1 to wheels
   Set_Linear_Velocity(w,VelX(5),VelY(5),VelZ(5))
   Position_At_Object(w,5)
   Place_On_Matrix(w,1,Object Size Y(w)/2)
next w
Place_On_Matrix(5,1,(Object Size Y(5)/2)+(Object Size Y(5)/5))
Calculate_Offsets_AB(1,5,0,315,60) : Calculate_Offsets_AB(2,5,0,45,60)
Calculate_Offsets_AB(3,5,0,135,60) : Calculate_Offsets_AB(4,5,0,215,60)
for w=1 to 5
Get_Old_Position(w)
Get_Object_Direction(w)
Get_Traction_Force(w)
Get_Drag_Force(w)
Get_Force(w)
Get_Speed(w)
Get_Acceleration(w)
Set_Acceleration(w)
Set_Velocity(w)
Get_Angular_Acceleration(w)
Get_Angular_Velocity(w)
Update_Coordinates_Slide(w)
Update_Angles(w)
Get_Lateral_Forces(w)
Fully_Update_Object_Position(w)
Update_Object_Rotation(w)
Get_Movement_Angle(w)
next w
if upkey()=0 and keystate(17)=0
   if abs(MoveAngY(1)-AngleY(1)) < abs((MoveAngY(1)+180)-AngleY(1))
      angleY(1)=curveangle(MoveAngY(1),angleY(1),10.0)
      angleY(2)=curveangle(MoveAngY(1),angleY(2),10.0)
   else
      angleY(1)=curveangle(MoveAngY(1)+180,angleY(1),10.0)
      angleY(2)=curveangle(MoveAngY(1)+180,angleY(2),10.0)
   endif
endif


`Camera
ca#=curveangle(angleY(5),ca#,20.0)
cx#=positionX(5)-sin(ca#)*600
cz#=positionZ(5)-cos(ca#)*600
cy#=positionY(5)+300
position camera cx#,cy#,cz#
point camera positionX(5),positionY(5),positionZ(5)

`Create DashBoard
dec fuel#,movespeed(5)/2000
if fuel#<=0.0 then fuel#=0.0
ink rgb(0,0,0),0 : box 0,648,1024,768
create_meter(movespeed(5),60,708,5,30,rgb(255,0,255),rgb(255,255,0),"0","90","180","270")
create_meter(fuel#,964,708,5,30,rgb(255,0,255),rgb(255,255,0),"0","90","180","270")
ink rgb(255,255,255),0
text 110,652,"Speed: "
text 914-text width("Fuel: "),652,"Fuel: "

`**End loop**
Sync
Loop


`**Functions**
Function Start_World_Physics()
Dim Gravity(0) as float
Dim Drag(0) as Float
Dim Windx(0) as float
Dim Windy(0) as float
Dim Windz(0) as float
Endfunction

Function Start_Object_Physics(o)
Dim Positionx(o) as Float : Dim Positiony(o) as Float : Dim Positionz(o) as Float
Dim OldPosX(o) as Float : Dim OldPosY(o) as Float : Dim OldPosZ(o) as Float
Dim Anglex(o) as Float : Dim Angley(o) as Float : Dim Anglez(o) as Float
Dim MoveAngX(o) as Float : Dim MoveAngY(o) as Float : Dim MoveAngZ(o) as Float
Dim LongForceX(o) as Float : Dim LongForceY(o) as Float : Dim LongForceZ(o) as Float
Dim TracForceX(o) as Float : Dim TracForceY(o) as Float : Dim TracForceZ(o) as Float
Dim DragForceX(o) as Float : Dim DragForceY(o) as Float : Dim DragForceZ(o) as Float
Dim AngForceX(o) as Float : Dim AngForceY(o) as Float : Dim AngForceZ(o) as Float
Dim AngVelX(o) as Float : Dim AngVelY(o) as Float : Dim AngVelZ(o) as Float
Dim AngAccX(o) as Float : Dim AngAccY(o) as Float : Dim AngAccZ(o) as Float
Dim VelX(o) as Float : Dim VelY(o) as Float : Dim VelZ(o) as Float
Dim AccX(o) as Float : Dim AccY(o) as Float : Dim AccZ(o) as Float
Dim DirX(o) as Float : Dim DirY(o) as Float : Dim DirZ(o) as Float
Dim OffsetX(o) as Float : Dim OffsetY(o) as Float : Dim OffsetZ(o) as Float
Dim MoveTime(o) as Float : Dim MoveDist(o) as Float : Dim Power(o) as Float
Dim Movespeed(o) as Float : Dim MoveConstant(o) as Float : Dim Acceleration(o) as Float
Dim TractionConstant(o) as Float : Dim AngularTractionConstant(o) as Float
Dim Mass(o) as Float : Dim Weight(o) as Float
Dim KineticEnergy(o) as Float: Dim WorkDone(o) as Float
Dim LatForce(o) as Float : Dim PushX(o) as Float : Dim PushY(o) as Float : Dim PushZ(o) as Float
Dim Force(o) as Float
Endfunction

Function Set_Standard_Gravity()
Gravity(0)=9.8
Endfunction

Function Set_Standard_Drag()
Drag(0)=0.43
Endfunction

Function Set_World_Standards()
Gravity(0)=9.8
Drag(0)=0.43
Endfunction

Function Set_Gravity(g#)
Gravity(0)=g#
Endfunction

Function Set_Drag(d#)
Drag(0)=d#
Endfunction

Function Create_Chassi(o,w#,h#,l#,mv#,t#,a#,m#)
make object box o,w#,h#,l#
Start_Object_Physics(o)
Set_Object_Constants(o,mv#,t#,a#,m#)
Endfunction

Function Create_Tyre(o,w#,r#,mv#,t#,a#,m#)
make object sphere o,100 : scale object o,w#,r#,r#
Start_Object_Physics(o)
Set_Object_Constants(o,mv#,t#,a#,m#)
Endfunction

Function Set_Object_Constants(o,mv#,t#,a#,m#)
Set_Move_Constant(o,mv#)
Set_Traction_Constant(o,t#)
Set_Angular_Traction_Constant(o,a#)
Set_Body_Mass(o,m#)
Endfunction

Function Set_Move_Constant(o,c#)
MoveConstant(o)=c#
Endfunction

Function Set_Traction_Constant(o,c#)
TractionConstant(o)=c#
Endfunction

Function Set_Angular_Traction_Constant(o,c#)
AngularTractionConstant(o)=c#
Endfunction

Function Set_Body_Mass(o,m#)
Mass(o)=m#
Endfunction

Function Weigh_Object(o)
Weight(o)=Mass(o)*Gravity(0)
Endfunction

Function Set_Linear_Force(o,x,y,z)
LongForceX(o)=x : LongForceY(o)=y : LongForceZ(o)=z
Endfunction

Function Set_Linear_Velocity(o,x,y,z)
VelX(o)=x : VelY(o)=y : VelZ(o)=z
Endfunction

Function Set_Angular_Force(o,x,y,z)
AngForceX(o)=x : AngForceY(o)=y : AngForceZ(o)=z
Endfunction

Function Set_Angular_Velocity(o,x,y,z)
AngVelX(o)=x : AngVelY(o)=y : AngVelZ(o)=z
Endfunction

Function Move_Object(o,spd)
VelX(o)=dirX(o)*spd
VelY(o)=dirY(o)*spd
VelZ(o)=dirZ(o)*spd
Endfunction

Function Push_Object(o,spd)
LongForceX(o)=dirX(o)*spd
LongForceY(o)=dirY(o)*spd
LongForceZ(o)=dirZ(o)*spd
Endfunction

Function Control_Object_Using_Forces(o,push,turn)
if upkey()=1 then Push_Object(o,push)
if downkey()=1 then Push_Object(o,-push)
if leftkey()=1 then set_angular_Force(o,0,-(MoveSpeed(o)/turn),0)
if rightkey()=1 then set_angular_Force(o,0,(MoveSpeed(o)/turn),0)
if upkey()=0 and downkey()=0 then Push_Object(o,0)
if leftkey()=0 and rightkey()=0 then set_angular_force(o,0,0,0)
Standard_Update(o)
Endfunction

Function Control_Object_Using_Velocities(o,spd,turn)
if upkey()=1 then Move_Object(o,spd)
if downkey()=1 then Move_Object(o,-spd)
if leftkey()=1 then set_angular_velocity(o,0,-turn,0)
if rightkey()=1 then set_angular_velocity(o,0,turn,0)
if upkey()=0 and downkey()=0 then Move_Object(o,0)
if leftkey()=0 and rightkey()=0 then set_angular_velocity(o,0,0,0)
Standard_Update(o)
Endfunction

Function Get_Movement_Angle(o)
MoveAngX(o)=AtanFull(VelY(o),VelZ(o))
MoveAngY(o)=AtanFull(VelX(o),VelZ(o))
MoveAngZ(o)=AtanFull(VelX(o),VelY(o))
Endfunction

Function Set_Next_Angle_X(a,b,smth#)
AngleX(a)=Curveangle(AngleX(b),AngleX(a),smth#)
Endfunction

Function Set_Next_Angle_Y(a,b,smth#)
AngleY(a)=Curveangle(AngleY(b),AngleY(a),smth#)
Endfunction

Function Set_Next_Angle_Z(a,b,smth#)
AngleZ(a)=Curveangle(AngleZ(b),AngleZ(a),smth#)
Endfunction

Function Position_At_Object(a,b)
PositionX(a)=PositionX(b)
PositionY(a)=PositionY(b)
PositionZ(a)=PositionZ(b)
Endfunction

Function Get_Old_Position(o)
OldPosX(o)=PositionX(o)-VelX(o)
OldPosY(o)=PositionY(o)-VelY(o)
OldPosZ(o)=PositionZ(o)-VelZ(o)
Endfunction

Function Get_Object_Direction(o)
DirX(o)=sin(AngleY(o))*cos(AngleX(o))
DirY(o)=-sin(AngleX(o))
DirZ(o)=cos(AngleY(o))*cos(AngleX(o))
Endfunction

Function Get_Traction_Force(o)
TracForceX(o)=DirX(o)*MoveConstant(o)
TracForceY(o)=DirY(o)*MoveConstant(o)
TracForceZ(o)=DirZ(o)*MoveConstant(o)
Endfunction

Function Get_Drag_Force(o)
DragForceX(o)=-Drag(0)*VelX(o)*MoveSpeed(o)
DragForceY(o)=-Drag(0)*VelY(o)*MoveSpeed(o)
DragForceZ(o)=-Drag(0)*VelZ(o)*MoveSpeed(o)
Endfunction

Function Get_Force(o)
Force(o)=Sqrt(LongForceX(o)^2+LongForceY(o)^2+LongForceZ(o)^2)
Endfunction

Function Get_Speed(o)
MoveSpeed(o)=sqrt(VelX(o)^2 + VelY(o)^2 + VelZ(o)^2)
Endfunction

Function Get_Acceleration(o)
Acceleration(o)=Sqrt(AccX(o)^2+AccY(o)^2+AccZ(o)^2)
Endfunction

Function Set_Acceleration(o)
AccX(o)=LongForceX(o)/Mass(o)
AccY(o)=LongForceY(o)/Mass(o)
AccZ(o)=LongForceZ(o)/Mass(o)
Endfunction

Function Set_Velocity(o)
VelX(o)=VelX(o)+AccX(o) : VelX(o)=VelX(o)*TractionConstant(o)
VelY(o)=VelY(o)+AccY(o) : VelY(o)=VelY(o)*TractionConstant(o)
VelZ(o)=VelZ(o)+AccZ(o) : VelZ(o)=VelZ(o)*TractionConstant(o)
Endfunction

Function Get_Angular_Acceleration(o)
AngAccX(o)=AngForceX(o)/Mass(o)
AngAccY(o)=AngForceY(o)/Mass(o)
AngAccZ(o)=AngForceZ(o)/Mass(o)
Endfunction

Function Get_Angular_Velocity(o)
AngVelX(o)=AngVelX(o)+AngAccX(o) : AngVelX(o)=AngVelX(o)*AngularTractionConstant(o)
AngVelY(o)=AngVelY(o)+AngAccY(o) : AngVelY(o)=AngVelY(o)*AngularTractionConstant(o)
AngVelZ(o)=AngVelZ(o)+AngAccZ(o) : AngVelZ(o)=AngVelZ(o)*AngularTractionConstant(o)
Endfunction

Function Update_Coordinates_Straight(o)
PositionX(o)=PositionX(o)+(DirX(o)*Movespeed(o))
PositionY(o)=PositionY(o)+(DirY(o)*Movespeed(o))
PositionZ(o)=PositionZ(o)+(DirZ(o)*Movespeed(o))
Endfunction

Function Update_Coordinates_Slide(o)
PositionX(o)=PositionX(o)+VelX(o)
PositionY(o)=PositionY(o)+VelY(o)
PositionZ(o)=PositionZ(o)+VelZ(o)
Endfunction

Function Update_Angles(o)
AngleX(o)=Wrapvalue(AngleX(o)+AngVelX(o))
AngleY(o)=Wrapvalue(AngleY(o)+AngVelY(o))
AngleZ(o)=Wrapvalue(AngleZ(o)+AngVelZ(o))
Endfunction

Function Get_Lateral_Forces(o)
LatForce(o)=abs(AngVelY(o))*(MoveSpeed(o)/Mass(o))
If AngVelY(o)>0.0
   PushX(o)=(sin(Wrapvalue(AngleY(o)-90))*LatForce(o))
   PushZ(o)=(cos(Wrapvalue(AngleY(o)-90))*LatForce(o))
endif
If AngVelY(o)<0.0
   PushX(o)=(sin(Wrapvalue(AngleY(o)+90))*LatForce(o))
   PushZ(o)=(cos(Wrapvalue(AngleY(o)+90))*LatForce(o))
endif
Endfunction

Function Calculate_Offsets(o,sax#,say#,out#)
OffsetX(o)=(Sin(Wrapvalue(AngleY(o)+say#))*Cos(Wrapvalue(AngleX(o)+sax#)))*Out#
OffsetY(o)=-Sin(Wrapvalue(AngleX(o)+sax#))*out#
OffsetZ(o)=(Cos(Wrapvalue(AngleY(o)+say#))*Cos(Wrapvalue(AngleX(o)+sax#)))*out#
Endfunction

Function Calculate_Offsets_AB(a,b,sax#,say#,out#)
OffsetX(a)=(Sin(Wrapvalue(AngleY(b)+say#))*Cos(Wrapvalue(AngleX(b)+sax#)))*out#
OffsetY(a)=-Sin(Wrapvalue(AngleX(b)+sax#))*out#
OffsetZ(a)=(Cos(Wrapvalue(AngleY(b)+say#))*Cos(Wrapvalue(AngleX(b)+sax#)))*out#
Endfunction

Function Place_On_Matrix(o,m,h#)
PositionY(o)=Get ground height(m,positionX(o),PositionZ(o))+h#
Endfunction

Function Update_Object_Position(o)
Position Object o,PositionX(o),PositionY(o),PositionZ(o)
Endfunction

Function Update_Position_With_Push_Forces(o,a,b,c)
Position object o,Positionx(o)+(PushX(o)*a),Positiony(o)+(PushY(o)*b),Positionz(o)+(PushZ(o)*c)
Endfunction

Function Update_Position_With_Offsets(o,a,b,c)
Position Object o,Positionx(o)+(offsetx(o)*a),Positiony(o)+(offsety(o)*b),PositionZ(o)+(offsetZ(o)*c)
Endfunction

Function Fully_Update_object_Position(o)
Position Object o,PositionX(o)+PushX(o)+OffsetX(o),PositionY(o)+PushY(o)+OffsetY(o),PositionZ(o)+PushZ(o)+OffsetZ(o)
Endfunction

Function Update_object_Rotation(o)
Rotate object o,AngleX(o),AngleY(o),AngleZ(o)
Endfunction

Function Standard_Update(o)
Get_Old_Position(o)
Get_Object_Direction(o)
Get_Traction_Force(o)
Get_Drag_Force(o)
Get_Force(o)
Get_Speed(o)
Get_Acceleration(o)
Set_Acceleration(o)
Set_Velocity(o)
Get_Angular_Acceleration(o)
Get_Angular_Velocity(o)
Update_Coordinates_Straight(o)
Update_Angles(o)
Get_Lateral_Forces(o)
Fully_Update_Object_Position(o)
Update_Object_Rotation(o)
Endfunction

Function Update_Speed(o)
if MoveSpeed(o)<0.9 then MoveSpeed(o)=MoveSpeed(o)*0.0
Endfunction

Function Get_Movement_Time(o)
If MoveSpeed(o)<>0 then MoveTime(o)=Movetime(o)+1.0 else MoveTime(o)=0.0
Endfunction

Function Get_Movement_Distance(o)
MoveDist(o)=MoveSpeed(o)*MoveTime(o)
Endfunction

Function Get_Kinetic_Energy(o)
KineticEnergy(o)=(Mass(o)/2)*(Movespeed(o)^2)
Endfunction

Function Get_Work_Done(o)
WorkDone(o)=Force(o)*MoveDist(o)
Endfunction

Function Get_Excerted_Power(o)
If MoveTime(o)>0.0 then Power(o)=KineticEnergy(o)/MoveTime(o) else Power(o)=0.0
Endfunction

Function create_meter(var#,speedox,speedoy,speedow,speedl,col1,col2,m1$,m2$,m3$,m4$)
ink col1,0
circle speedox,speedoy,speedl
line speedox,speedoy,speedox+sin(wrapvalue(int(360-var#)+180))*(speedow/2),speedoy+cos(wrapvalue(int(360-var#)+180))*(speedow/2)
line speedox,speedoy,speedox-sin(wrapvalue(int(360-var#)+180))*(speedow/2),speedoy-cos(wrapvalue(int(360-var#)+180))*(speedow/2)
line speedox+sin(wrapvalue(int(360-var#)+180))*(speedow/2),speedoy+cos(wrapvalue(int(360-var#)+180))*(speedow/2),speedox+sin(wrapvalue(int(360-var#)+270))*speedl,speedoy+cos(wrapvalue(int(360-var#)+270))*speedl
line speedox-sin(wrapvalue(int(360-var#)+180))*(speedow/2),speedoy-cos(wrapvalue(int(360-var#)+180))*(speedow/2),speedox+sin(wrapvalue(int(360-var#)+270))*speedl,speedoy+cos(wrapvalue(int(360-var#)+270))*speedl
ink col2,0
center text speedox-(speedl*1.5),speedoy,m1$ : center text speedox,speedoy-(speedl*1.5),m2$
center text speedox+(speedl*1.5),speedoy,m3$ : center text speedox,speedoy+(speedl*1.5),m4$
endfunction

Zotoaster - generally LOTS of updates to the entry and also good racer physics.

Ric - VERY close. Good circuit. This could have won if it had more development time

Riidii - DAMN good physics there. Shame it couldn't have evolved into a fully fledge game rather than a tech demo though