clear

theta = 0.2 ;
thetaprim = 0.0;
grav = 9.81 ;
mass = 1.0

nstep = 20000 ;
deltatime = 0.001;
time=0.0;

for i = 1:nstep	

  thetabiss = - grav*sin(theta);
  thetaprim = thetaprim + deltatime*thetabiss;
  theta = theta + deltatime*thetaprim;
  time = time + deltatime;
  tid(i) = time ;
  vinkel(i) = theta ;

end

plot(tid,vinkel,'-')
pause(0)

