# Untitled

Oct 30th, 2020
602
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
1. for nig = 1:NoOfShifts
2.
3.     if nig == NoOfShifts
4.         DragTimeDistance(1,nig) = vpaintegral( fMig(GBGeometricRatios(nig)) /sigFx(GBGeometricRatios(nig)),...
5.             symVx , fNe2V(Ne(1) , GBGeometricRatios(nig))/3.6 , 0.99*MaxSpeed/3.6 , 'MaxFunctionCalls',Inf);
6.
7.         DragTimeDistance(2,nig) = vpaintegral( symVx * fMig(GBGeometricRatios(nig)) /sigFx(GBGeometricRatios(nig)),...
8.             symVx , fNe2V(Ne(1) , GBGeometricRatios(nig))/3.6 , 0.99*MaxSpeed/3.6 , 'MaxFunctionCalls',Inf);
9.         break
10.     end
11.
12.     solvx = vpasolve(fFx(GBGeometricRatios(nig))==fFx(GBGeometricRatios(nig+1)),symVx);
13.     UpShiftSpeeds(nig) = double(solvx(end))*3.6;
14.
15.     solvx = vpasolve(fFx(GBProgressiveRatios(nig))==fFx(GBProgressiveRatios(nig+1)),symVx);
16.     UpShiftSpeedsProg(nig) = double(solvx(end))*3.6;
17.
18.     DragTimeDistance(1,nig) = vpaintegral( fMig(GBGeometricRatios(nig)) /sigFx(GBGeometricRatios(nig)),...
19.         symVx , fNe2V(Ne(1) , GBGeometricRatios(nig))/3.6 , UpShiftSpeeds(nig)/3.6 , 'MaxFunctionCalls',Inf);
20.
21.     DragTimeDistance(2,nig) = vpaintegral( symVx * fMig(GBGeometricRatios(nig)) /sigFx(GBGeometricRatios(nig)),...
22.         symVx , fNe2V(Ne(1) , GBGeometricRatios(nig))/3.6 , UpShiftSpeeds(nig)/3.6 , 'MaxFunctionCalls',Inf);
23.
24. end
RAW Paste Data