clc
clear

x=0;
dx=0.01;

while 10*(1-exp(-x/4)) < 9.8
end

xmax=x;

x_vals=0:dx:xmax;
y_vals=10*(1-exp(-x_vals/4));

plot(x_vals,y_vals,'LinsWidth',2)
xlabel('Time(seconds)')
ylabel('Force (Newtons)')
title ('Force vs Time')