Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- %X = [x_r; y_r; x_d; y_d]
- f = @(t,X)[v_r*(norm([X(1)-X(3); X(2)-X(4)])>.01);...
- 0;...
- v_d*(X(1)-X(3))./norm([X(1)-X(3); X(2)-X(4)])*(norm([X(1)-X(3); X(2)-X(4)])>.01);...
- v_d*(X(2)-X(4))./norm([X(1)-X(3); X(2)-X(4)])*(norm([X(1)-X(3); X(2)-X(4)])>.01)];
- [tout, xout] = ode45(f,[0 T],x0);
- delta = sqrt((xout(:,1)-xout(:,3)).^2+(xout(:,2)-xout(:,4)).^2);
Add Comment
Please, Sign In to add comment