function a=accel_midt1(t,b)
	g=9.8;
	tc=10;

	a=zeros(1,length(t));

	a=100*(exp(-(b*tanh(b*t).*t/10).^2))-g;
	a(t>tc)=-g;
end
