function x_dot = f1(t,x)
global b;
x_dot = -b*x;