My teacher tried helping me build this code for a falling object with resistance. (We did it on paper since we didn't have a computer near at the time.) Given only the acceleration a = g+kv2. With y(0) = 1000, and v(0) = 0. Need to find the time, and velocity when object hits the ground.
I'm getting empty sets as my outputs....suggestions?
h = 1000;
t0 = 0; tfinal = 15;
tspan = [t0,tfinal];
options = odeset('Events',@landing);
[t,y,te,ye,ie] = ode45(@drag,tspan,[h;0],options);
function dydt = drag(t,y)
k = 1/10;
g = 9.8;
dydt = [-y(2); g + k*-(y(2)).^2];
check = y(1);
stop = 1;