ex6_8 - r3b=-10*8*(x(3)/ft)^2*(x(2)/ft);

Info iconThis preview shows page 1. Sign up to view the full content.

View Full Document Right Arrow Icon
% "ex6_8" % m-file to solve example 6.8 % % x(1)=fa % x(2)=fb % x(3)=fc % x(4)=fd % x(5)=fe % x(6)=ff % xdot(1)=d(fa)/dv, xdot(2)=d(fb)/dv, xdot(3)=d(fc)/dv, etc. % function xdot=ex6_8(v,x) f ft=x(1)+x(2)+x(3)+x(4)+x(5)+x(6); r1a=-5*8*(x(1)/ft)*(x(2)/ft)^2; r2a=-2*4*(x(1)/ft)*(x(2)/ft);
Background image of page 1
This is the end of the preview. Sign up to access the rest of the document.

Unformatted text preview: r3b=-10*8*(x(3)/ft)^2*(x(2)/ft); r4c=-5*3.175*(x(3)/ft)*(x(1)/ft)^(2/3); ca=2*x(1)/ft; c xdot(1,:)=r1a+r2a+2*r4c/3; xdot(2,:)=1.25*r1a+0.75*r2a+r3b; xdot(3)=-r1a+2*r3b+r4c; xdot(4)=-1.5*r1a-1.5*r2a-r4c; xdot(5)=-0.5*r2a-5*r4c/6; xdot(6)=-2*r3b;...
View Full Document

Ask a homework question - tutors are online