function z=suuntakf(t,x)
m=1;b=1; g=9.81; 
z=g-(b*x.^2)/m;
