From 482822af466cac7b97270f80d2859e7648fdfab7 Mon Sep 17 00:00:00 2001 From: Ayush Rathore Date: Fri, 17 Nov 2017 21:17:02 -0500 Subject: [PATCH] Problem 5 --- README.md | 11 ++++++++++- figure5.png | Bin 0 -> 22746 bytes sphere_drag.m | 26 +++++++++++++++++++++----- 3 files changed, 31 insertions(+), 6 deletions(-) create mode 100644 figure5.png diff --git a/README.md b/README.md index 250460b..17c04bd 100644 --- a/README.md +++ b/README.md @@ -3,6 +3,7 @@ ME 3255 Homework 5 ##Problem 2 + Least Squares function ```matlab function [a,fx,r2] = least_squares(Z,y); @@ -159,6 +160,14 @@ Output: - L = 5.0108 m #Problem 5 ###Part A +```matlab +cd_out_linear = sphere_drag(300,'linear') = 0.1750 +cd_out_spline = sphere_drag(300,'spline') = 0.1809 +cd_out_pchip = sphere_drag(300,'pchip') = 0.1807 +``` + +### Part B +![Drag Force vs. Velocity (3 Interpolation Methods)](./figure5.png) #Problem 6 ###Analytical @@ -190,7 +199,7 @@ tildex = (b-a)/2*x + (b+a)/2; f = (1/6).*(tildex).^3+(1/2).*(tildex).^2+(tildex); value2 = sum(tildec.*f) ``` -###3-point +##3-point ```matlab %Three Point x = [-0.77459666,0,0.77459666]; diff --git a/figure5.png b/figure5.png new file mode 100644 index 0000000000000000000000000000000000000000..61658e6484b1250519430d3b145afe62ee285a0d GIT binary patch literal 22746 zcmafb1z45a*6t!CR6syN+C_<^w6wUC21N zlD=cP_xaB`_utRGJkJ*PT3^g>&N;?-$2;Ef1>IAYAvjNU9)Um*$jPEq5eQ5k1Og)s z4--B~eigV0zwqs3A2}cpgiYxGFrwKBso_IhM>(ZCxIZrtUt+)DAU=BnA2B)J(R6%Z z`^3@3@abcO(yfPzx1>239-2L7P%?Bh)3AMN>A=AtXJgDIra4Xt$0O0l-+ubo(9zu1 zhT);P^<%^hAx>pmXHH(88$t{>Ziw*niQKsHYA3)8fnY$$p>92VmasbRF7@#7spKY+ zDJ8C5nxX$Kt`C$^&)WL&xM!mabjaA6EWHN4jsGgEavt&O)vze-nm%qr9^%oL%=}ljsXt7p@A@8y(FEm@)xy)pEOe^hsOH}Uw(e)YZ0?eYImih z$R3b~u^P=ZaltPHF(y7dD5ftq0ukvW1(PHv#ejkEVPc3tAf!kx$-+SK@t6>Z2yAQ; z1VR?`&*2#Vakw2>HKKbZ=TY(euetX2815uK(=fZu&%Rh2bw}&)i0ul-m|2v z%z1x(+{VTx@yM*2eR0xZ@W`mBmF4AxL;Ku=!*SD` zm9_q8y;^rO=y?mB`}Zt*-0qurOp6x#KODMRUo!?RhH%#z$kf5eibuqre1G@zO91J8 zC#T&F?=yFOef?U1$uHDn4jY*eD&-d@&-+p<~btye&d1Yl}>*dRrF?w$1 ztr1t@iw`p4r<-4hMz(||o*9poS+8bhswAsTsC{FBQjw$hShgJ*u$ShlE*y6v|`WfzP(|C%fb(l*rt?)g$1?C=YfGz zl_UxHsQoSHScy655fa7kzPCC&JPc>u%T-s27Yq&TOp-i{Ca)65B~>Jqu-(*;N{8B-n!v?A70J@t#>N*H$j+Z5-i5uiyF9Q_yH7{l5hoCG zau`M{8o0RCLSN1E`$?>sAr4Ah>^ z1am^&yLa#NsrbWt_v~iA66~#3yX}4|UBVE5@}1Io9(Ir{D!=};YIn%j@!;3vSVfko zwxZ(VT=nb*gNp4=(V%P0%rA6d*oop-9;m9)xhl@1Fz>MO9G^(T8{RaI5<`ZcZ$ z(J?W_+J=JhBVXb89Cb@;>)Zx|;ac~7{nK9*&hyjFAwi4VU6SQmxg#SZWs^4~CG}aN zN{MC5G}p`Qv2vBNWiK$Bue~kX3OLt|Px%-ZNy%^4-4sOH*ViX`cD(J2g&Wr~9~g)S z=71EjpT;Iu)6h8EE3aop?p^!-{d-BZw3ZfC*Y@kzufaNmh+SP>@7%ewHe5Kfl3)L1 zO&LBvIofAU9Vs^LdM{qR_UqZRXJUtoAF7?#4wtj@pL8a~zI{u6`7*zGFH?2j;m%UA zQFAb@sPpPbF)J_cul}qDjCWsEoShyiC@8=_>h0})NkBcnwkB@;^_5z-`tDL+CS$pD z>ywl#)(y|kp=RgNKY}jnKW1lVS>HdK{vOMhshmjSwF3JPwKL^MRJ8H2Z|NlweW&WD zJ1UYVJ0I_V7fsuu-0@z)(*lGDM6=YO zu&^gYbuTVZoX3Udc%*V#<+S3wHribFsO9B&^2^ENc_9?P)o?*HyGD^-Qbfe1L5DBS zw{PFp$k(GKA(2MK3tAUdV;=1WfCY94O`1>~_a&V+edF^w{w)}M0uNl+_m<;PV^~^R z+NIZbKO8K62#<`6)SM)7FB&sEN6Gg@(Zs~W%{g@ALUtP5l z`|CFoUT^QSQp=&D9-6@? zH$HD4Glq*b5)x7izN)?V{cBC(M(fw}het==hfA3p&j{%ytoDPfu~n&s3S1dVofWi? zzF;n5Q&{$Sdmti>cw-^1uFuRcI_8zM8Wp~I`*sgP@Bzc~!^6X(B3?m%>|)gR#j|7U!T*!$QrP+S<^& zGiX%Pndts|^JkaD(Q01x`nUFW(rPnAVt^E8viHe~{_5P18^&Kfj@HV(yu6HJK6>I3 z(r6YxQS$PV7^CHT@;y(tdcAUi&R*f1ohTErtXr*wkdQ|9KvAcvD30IH-gJdVuK1d? zWWe85gN4IZcEqNQvhvYCrlgoby5v4%%%|?MV8+_zHUvaOUEjaI3JcRKw@Dnw z2Lr>$$G1wQK#qd<2Dhp$Hk>~9@_qA{FSl2bjyS1!ZX~ZQWyQJzUc6vtW?pg>@oTt5 zOuRLqp4Rdqz>(TO=eZ;S|a^H)&K6P`e8ovq7 zm4J{?kn&=qiGD8aUx%_(fBW_=<#xdNYWehDn|GJGw!zlpg>2!d=;-JhH*RQ`S($Gc zY=@K&GoC|qTh7hT2kZKm1=7H9#_B3Q+b4{_(E;!UF1$NOk-gK2+FjUZ z8$T+el;oyYRzy7ycOZC>6FJ+rVIV#iq6nEn?>fwXzhXFK#$Qy)a_+^87w-G(WX>xq zE01{}h2T~vmmm;RVn|$7O-)MT!rk}V+uMA+yex|Eq)_0;08d>-Mzvk0r|-&G?T7M}_AG&l?u{c|~@E}9}k$>a)ZWG)v!z{M4F@xWu zMzQY61?{fb=g%%63@-R$uNhx4&SIOl42dsO7=CFtt6r!WS&93U&Npvj#$+A+FRhjZ zV)4l(S!G8?z8)R>M%9z-CCdZ-jnT86mjm-MVcaB{HIp0>ef&ijiHuB-_~(t*OYm zwqB?>Fql_Ut1x;_^XQIQaBNudyH=y+;n)FndaKipSa@kR8Kf9}H}9}) ze16GGF6bdZ{5pN~+b0*eZrz*?2g7v+5tjX)Zvt8xNmAj&KRHrIe={g3-MEnX@>kT` zI3tH<{GY-XSeZy2UBK`!kaO@*g@SQ2Fiid4nvsx@u)3i-4}U@210->U-|Uj8QK2Td zKpkD(wS@)nRlWA-Ipx<@BthmL+JI2e+-!T2yB!!X!Tca_lzF}UTOG~zQS_r(BoWhZu zB9@%;eeZCdx3*GI1*Ikv6)3y@ff>0^Ose&rwUaSo8_$th)y#Ngh+iY2fPi((1a94; zjAb^IN1s>YIJJ;>_G51UpjncCcu)4plljy;A%NtIy-wV?x%0u~=^%M%4koj=w?|zj zB{k-!YI${K8*(0N2uc8xSfbd(#5{%z3?V?ZK+Z$Kp~++3n?^t_SW#64c^)?s`FzH0 z9$}6hNh%QHtG_0HJcye$9%xoMYU}!ro|ugCB^jmCxtF8;A<;72R2MF0n9MQ#ob`G6 z6iH<%B7>S)?mWl;rqKH3i>QN-dabPF3UrK(@fz=i!=+xxa^8zEZqkb8|E7#{9Ij+1a)z z)DB>haGIM>;dvIil9EJSxI(f$kM`C8<@Kh?!AIQO-1|E_myZ;LUpFNMGmzb~C}$|8 zK9?Q$tL>oyf4UwIU9G0;PHJJ(KzyJk3jY>a_dSV0N3MbtE35IF)gHVtOO{W0uP=r= ze`ympsdaC4Q&6(EUMCzZ?=LxFHo_kj3_OuxX9j4iaXlK1I@PtbAECba1(3gfC*?s%d9gN}OTV)+v{AiL=RO_8s^`8|dhP-_#OS)Muhgj71=ie+>5F(wc0LTH@!b_N z*7Kwx*9O=4I~_aAewk1xbHs{ec2G#XP{}8?GKrio7Q#rtS#siU^mPMlj`16GfLF7`grcT z>C~kt7A37xi-FCbU%oUp!o4JaPhn|j3}wD|vw)l!jS3)x1YGZSa_E%3J235gv7&tF$tRO{8>j>q~cR2{G}TnwAw`;E+7$<4$M&E*?{6-^y-fWu;bbS#!F? zMVv3`3NilUSCyClHkhmZ<;xd7KE4khJ`Ck*o4QCpefm^d`nhz9 zLA9vMdV;XMQ49qsO-^WMrVzKFvv8cR|C{iv>MMRiglf0Vt z+VnS_VHuwK4I^E7={bCaO~7TWpoN3qzkfrXkf;nMh(sb&KYZ|e{@khy&>~V^zCN)1xqzPipfk${s!mY&sy+x zWqzH`Xn&`EA^4VA-FiLD3y;ZXjImkjQ%(+zh~vkC?1$I4fBmvjnR2cB_|Xh<5FEc0 zfUqnqPk>e7)RSNIj(pb6_T!WOLSF&q7e6q|zHh*`2wKYZWxswF(LC<%dRHu2%Lj`6EE?qg=p-ucu z$I6*ol}ppnpFYT`$%w`&<~TQ(MSmDu4w25Ll}*^BeB;F{Xa7y}g>=JqO@U>hrjD*R zS;Fb*DeQ=y{e4!XzWyZS%BD}AI6JZ3YMiCAZhdD?q@>&MM0?>4mrX-L%QX!oZ;U{l z?a}r(donrOclk@&HI1Q3?;lfx^BE1hVi@SUUH%3kG7S(A?a;?NLaXuFM)Cy=ks+0%iZ8czj1hE@U37tYl?nb##8f!~x3}o0wSZyrxy-k8N#7^r*`fyGXs_ z+3XlLtu1QN|Iuk>&~vi^m!<_C zFo53c*AU9D!HnH0PB%U@I56s;%=CoRANB2z@F`ozPq|H4w~}2AjjMO+$9Yv6#bP*T zL`+6T=6QCq50O&8$fPZhn9*y$dJ-MUQ&XppHtJbZVdK1)@N|PnNkBj_>2;KwnaQeM zVGBdB)=A0aqELt=U_t~KJFpi7Z46rxzONLe+rvu`GCxii6|d#+^yjlI`n__dNAbH? zuk!9#L&*LWVG855w7cscqK_PRRJ!F@l;Zm9y}i*?Hp=nz_<-K)_s{q4>xPn{uh2hY z&o&LtPEST=O#uxB=;hlUo&OR*v@vbl^Z1f%cSwODAvlw46-joE5Iv1zq2wihWYj)w zmd5vVa-^rdM^5)%FIUPyr;e{w zLgiJq*@bfprdJBpBz*=e+`{)SPWC}`&!G*cH3{kcgbY6kpHD1b4@}80d?uMsGc?+5 zr`?EshhUC1QRC^NbKTjgyWp_}&@{lu0IG$|3`p^TfdSOcQfAWWl8SfL0K4O@!sh2x zBs+H#WRe>h-_}YqeJ55%-9jbhdfIKa%OkZMPLfb2TzwQREzW5NgyXW^394B-6?Q3i zKiqU#*HKj^Ar7S$S`T7oTHqprl_q%U=A>^la4Th>;e=o6?`d)C+u z8z;3KKe%~do6e03kWbN~{2*)_FYUaMA4`E43D5fax+H8OC#O=?tjmO_TjAawKy$qy zJHW)m^xv9NhrACPdmT{e|iUONNs{MJt4Y!2w zyTJ1oGlqlBbFhLQM~?q+xrj50@q-1v&wl&jK#5=vnTqMRbfsPKYE)cT%Qu6!YCdDW zQkZ6B)b=gp4<>iNM89bWvF6Jh#=J30ZULSdwQXm7@;J z9u~||Jl%vMKfyIXdR6O8z@!b&<*Z1cj-M#kclyg4%)8@>tno_M==&gF=Wcv}!DpxY z_3mp_vSsjCYsG7dZgZJ0KgUoF1(Z-vF}$)Y=zqgskGCTY-ftOez>FRhKLI`g{9%r6wVpeGQR)#q|Q?)d+-$h3h;0l!l#Y-z3OZS|lW040e-3W{f)? zAC4w$SUS|X(BJS4)TwPLeDRv+L*gQ2D)wi1Xv7P?oSr@jm)^r}r06l&+|8ReJBEZT zPD8+>nL3;UgctqxElE#Hb+1B*CZZFouo@}y_s53&&aY7G)_IJq zQ^ad{LJ0vNZ(KygXJ^2}e&FJ7mTsmlnhQuEGK~matlWpfHIyp}wkg5f1#Vuka>T+J z^_`lF(bt1;0Q7W{>wIirVBquTn9Hk&cY*p6cH6nuVh_*3eeGHWkjw6R?zAvhHqjAs z=WzOLiO!HYat(S4jo8u*&z7??p%lTY6vq$L51QH4Z5kAcyXi0wge7`Ix@@Zc8?=K- z@tJle5J~+0JWz^TW?QM5t6tbX#c}oORhTcM7)&?x5|pn5u%k-EzT0kBp5?Z+F{EN- zlZt0{t)F+2O3lPFkN zXMt|g^u~+4krb7CeIC-loXBDthuwPaHg=O`8&P*u8xiYi&A(P>x8$_(Q_Lc8#6PvQGc_A;V&;5dA-Jh`J$}x)pYYZIR60 zr`t)4jEq6oUUXU$+@t zxWVU>cXdnI-7LPl{5J$FG9(nr>E!Ety|8Su^sbZGNWT-+QLvxC8eM! zf3`Kv#mObbm&0ei*?t9b(zb7H;SV$MB`vN=%KdH;kqQX<5%3QXUPzMj>eZjVxO9}f zCV_UdKuR<;G~oC(SP7QfpG2?_a^Lik(p)|65H$)6%wI3eSCw+<52f=1V(;02 zbkkL1!LD68v!K@aUhDuMpq_B!w#WVZJ@Xy$=pqt;84A5;DWxk5g7(kaFsur9*&tk~ zR1z>Ex798f&I-xi7IM27{D;O=_G$PVc7OV0(U+mrcu`pq2vif3S->7iuV34PmoRP# zTO08cx^zpi(P$wj5*v}o@RWh=^7;1*FQB+$emnv6%`d9zp}AW@35fk4JhN+?rtK$nI})4Xi(5j4GkMDWd3-W_>}rdr6KdhBbv<~yH3M~1lpgv)(vbnWYX+QIIeZ1ObV?q~}CoV4T_3QhZniVb^9kJ1WaiBwF zRLNvkh1)_Kfj_gU!};GB2dnvnk+{CmMKhx585zp4yv&l4l8lVT;Py~&qNCG$%#M$b ztz>*94`!p_vyfM>zJn!%htrF@l_+Pj|Dj2{f%!xo5h4U3w9PF{JZ5XC!12*$dCFi< zQx1d*K>S6=s)U3D>#?%J!ts9*q}HySQ4tXl6%`f05&{pyD=O-~Hd+cWkqIfgR%SFf zoSJwN(C37lwhU?Z9wJJP|77ct_#%QGe;E>Q&T z1=Z@phDnr378`Ln4`Lh-6aVwpWxx>x7fY}fg|SI|W}++Hw0ZQgK66SlU)J#zG`euQ9=ZVaNKb^d0|~L zictc*@6_d67D9$^5-;XDBWOVaV)s@Q%hAg3LDjoReu%Hvk! z`OB!cR6AEON^bun6(B^E=I`&1OZSX{l{GU|;e|gx7vs21AH5IXWrq)eb>!8TyR-4b zZ!JX`F{r7l6Vi%299Fv?4b^Aog#;0&W#9-8C{nLGgpw8SAFn+kGBBHDoZ@X|H$tL1!FzN=+$P4w1NU2q`Z<+ zBJX2)02_Djeg&r&bd8-|;X}{wgM(b{a%JIzHrAvBAzM8YlQgsDb7r~`qltO^HpX6R zK4o%+FPrRt?5&N#z0e@yP6#;|ZfNM5vZA7x*Kx&gI#z{14mEz10HQ=8SbzyhC1j(s zy?Ox3i1?&~3HJ>rxgdlHV{1F-Sig6&_D$jx*hAK;Iae`u#6TO+%>~_u-B69m zM-Z?=f&!(PY28aIK)0|HJ0FWa@~dFRh(H^w!;;dwGm7Nf=+Bw?M+_Tneag^hgg6kH z{rJr5iLeGEcBDm5Fcw&P<@5e^tJ^OPC;xmd9Mc`pw0kFt_Rn>d1D3$7K>;U ziuq=mHGKxrgIS-?mPK!NC-gTLepy%*75Bdh)Fwf&p$%F0exJLnzppbQ7yWjB7jcXz{X=!_TIfbGDP3UUWPjKKbMj{H?M8pVrN}jWi1)w~UXKnA1kp{`Hk79W zJKu{_)6iJ)P!4nXG>s>8=HYx;MQF|criyGfa$U(S2N_y5sF6m#e3|(gN(}{WTT51r z|FY^bRy&%^zNKgHCm}|H1wm1NZDpnXU$B4x=QgA}ii+Rga%yYU)OE8(Lrf0j#1S*opzqg=<3Zj$N*4F(=ZwRh4P!Mbhqa737 zLc+9n@7|HItJ6q&B@DR36iXa$enE>RO1DAb(6!x|hpvWNr}WOx%oIB=^`aqiN5@)} zip19b{$!;?W_`UR@xkGta*~8NA75XgKKAEAe-=?Xb_ffN4Ze;bnwmV_D=(8_QBJrP z)docdEiwHMJqb4*?hXr3g@VMQj8&(~aS3F2q|Ue(Ded3{fGXg2&g1^U5wXT~2pG;9L&qBbM#nh4QZBalwX)#DRYxXu3GbLtazV51zp3kZX_;lFVx+Szf#*?i;w-wnhyF*5RWL`*>KI1ks(q zlSm5)P|JT`4WlYfaS1@0S?8AT!v_zFp^lA)BY}Z3AU6daEa4fR&-yjlwG5$Mit5XV z-npqi! zTAhOS$i+oPH0Y%vULckl$ZY01=YzUg@zN1+F`(#SI?yk%qyjb2;&n;N<)O*1C)8pr zmI$-2Tb{BkQCSDXf~B-l+!V6|71v4bdSUS|Oi2X?2EG%reSeRhH@?GiFsH|i|J@q3 z56>eksdyS4dpVGl1_uKxzA{$M-C_^op^*6aF&I3}5gTIs&BW~5ZPcF`yy3e8su^d< zexMcw17K#plRp#_6QiB0UaCcV`SNT;IARe4$f<|&{t`?+Z-rSRfXLw_zd#K9j_2vY zf);2Q2SjFJC4;(6hJ2V`BjW5886Z9k4?|+_4lV|?NnCkxx9P$t#2m+mX^q2bEFg^aW`4j>|d zjAhrz`%!2tvrT#mAP%g14>Cj0pFIP?D9t7 z{fy6U`XxpW)XAGSX?;rnB}AN8-^r%~>ln}b_!}rg3z!%j zFNAJqqtWseT%w~>&=A{tcl?4fO|*P^-%_9~h-%-9-4que-}9fEngY{~mcN5~TA8rm zT0G(rPk+Aq-1Gj>jjL#=bd#+vU+dAMn;<_i)W~``KQ}k>9&|LynKQ0l8xyqYaY8d#0MY&=}rq^I>T)m}Xo=dIwJdF|*I1WELK2d{=3NMM)ui>7C zMrMy$7`1Ot%5F6I%UdLZ@5&GgR@g?_Wuqr zOO2SAj(YcXXU<_Tbka^^TzJgd$qKLlye9t1Eor0LVQ7F|p%L(G3o6 zZWR?3Q&xnme#B^Q*(*GN$8q*{QR5(IMMpJY1ktfa!czny0$Te8eiA4raxNW8YU;5R z-v5n(JI2Q114&3p!|5fJQ7+?7Ll?MfE|HRMH<7BK?%w6nt>Q-_8`0rb4jr^U{|OsY zb09ab1G)NO8P!8rRjA*?U9~yhf-gDQuRl|}cA4bQzCrYJdA})s8i?hD1S)5*3GH^<9#ba z5l{wx|E~Bj9N|xe&YD++A7|jeXY;nSS3?n$3Z-Ox0to(v-%@}vs=Q9Bct*lUtWIG< zKtn#WNrbtW0UT*6NgcTWlaIc=9Ro=1VG*{SacgU9iDeO(mEnUrzJyE&cS5?yF%c1) zFfwr5>K}S&KokxIC`!=y0_<5DR>4@r=nlX(Xld|$sHV0ndA2Rd8*BC12LozsAmR7{ z@h$?QwtJ|0# zeJ&e01XFW^ zZwaJmGcDn^<5f=G$+sYHgKCgBRL($&06-s8!D$8uQjRhjv2S@rMZa3;y-miTkk;AO z2Hb6zuFLq`Un5=qD#1kODT@%K(YcN0NwsZjnODOoZ(y2Z8_4AK6ovZQm=$4Fv~3l8 z5>fHD77k%A8VZLa%I+ITMbUm%BD5)T2hCy(m#!Pg-r>??dAEJvKs0;aMuo6~cn3WC zO{6@OXJlnB5q>i%0*EE@*P{>bOWk>a&cZI^;;LF$d@|<7>t-xGP`mId7*qL?c$D2rilfDDUT!sBNfG1Uye&Sca?Cl2F_>B{ z@a?XmOFZ}l&WKljled@s-@JJPxz@{o04Sf|m66GZyldT`Q@e~t!b2EPDQK>s`b+P9 zTqvC)tE{XHE5bqyU>>datgcx_lwSb%g-v)N4H!*M&4MX)?~s31Soo7dU&TJ;Cp4nY zis`Dv_d&x5Sqc=8;i^zT2Jc6!=H-*ew2fE7f61^%6G0t@&tz-6F*-UrJ^l8~Viy#- zrTU5r3J%Z_?ZLn@o(ed@a4RajK?6spFJd8%38+uMu98=wvk ze%l&Y5IU$70x%crq97v!joF9+&qH)<#s{~%B zdo5dHUG_rlpY3%l`v{?m0$m#(fKXxo0XGT7Sr4cigUB-sbY82&g-{-GRC|pq z^(oi?zB57)WTXF~chLS=^4Qc=@>eS(L@f{oFcTmAk(~BTP9_CHp$`Ja``Fl#$B%LB z@DxpP4P<8jFP3Kf=9Mc~!U6k2Psi_J>{f&^gI;sh+wBoXjQu^mr}N4{^@w3!5lhu5P;{qU%C zB97OfU8a8n!>*k}D90c84XeOoI+7Oqv2xqV0z*IL%=WgnkuvM}@^T?!D1CLt39t(a zqPw%eZqU^je|^`T9vLC)(SOKKA$fk^nc-3T#>@7T z{7Y!rkH*H95e2TB#F8Y*9h)Mo#nOt51GO97gg(H02Ko^ku^{PyC%OdfXHpz_N-JLp zn_6uVNFY=K??6HT89XJhJi;&|ma?#}f6lu7@hUp%aN_%*fbm2zzOqCh&?htxutq4E z2EdwKgsXW_4k4s~a4BpR&@vPe1upVW16z;Xy3!04ltY0cWk6PL+T>tHGeP0h=h^R! z`_AZFfd563%xg4DiT=5&ESKUSa6mwR2I4`2JdD26-8}P)Ou#8IzCdSKzq5Vsv!Up8Ogj?5QB&8QF|^)BMizmgT3T5DEdeK{r3Gxx zlwb!%9@vi4!-d9ZP9-%pRaI3L^2a;uH>6D4JWB|Mf*{?Kic@m_M>YoJ^m_$gci(tuGWEJvrZs2Rvt|AJ7#!&C%n zR>G(;G)D*)|FgM&OY+ViZ4*$aYC>T%9OMcirf>h#iY>RWXD`Vn@VIwFK>B#jmrmg? zq&jSpdGFe@_Y}^MPiI?RXaGgbNh@PGhi+9#=`yj;4t_B%@M=E5y8ail;k*D17-FZt zXVBAO94vasy#E2%4f?tZ|`MKW{a@Sq~lL1Y!whHfm zl{#Nw5(Gq~JnxdxHnB5nKU z_)D@+^qsdW#!< zQF7%a-m-NXt+l_-dCs${oxvoyu+J)Kz;=F77ls`P&R}Thg*8oWVEOB07kg-jMWO!H z=<(9uzr=nNS{gwn)Ko_Bmrj$O;0kXOgw4uyY&%{MjY2*8K@~Rt8<0|pKBjcuF`9!C z&+r`3&(lK@V2_Y^RCu4ALhI4`aN&p2z6ofr1QE^oYG6WvUOCc@Lf6Q3RW%W;MjKyw zXA~Fka-F-|19fgsy8^s8a3!k{hwEI3<7Lp|hE=-izRqxo-vqJT7Y|I=UVo zbNNqq3na6=b~L7inya0K;=XME9NkNvo%~=4;}rDa!}FBoO%JiN8wy{%`+N*dtU}TB zf0En0yfuwcNLQiBZ)syg1yxs92Wc5}-8dStVN^AByq=V&u%H_9foT!{uf`IfpURH3 z|Eskm_Q5w?2AtT)kLEg9#D=a+i_9p(KQ$PmliT120S}O}C^Do%U5s5r%Q_3Xn6|dI zfGuN41ycBS03k3WP5??0n%BoQ1@nBnEpTgJLY@8pnSV}<+JZ-IYHmi8NFXYyKj^Sk z&(XYLKi$N|A>dP7z^49)o`iZx%B<8+?@O~qX|OJiDX4>%$WNXjp6iEk^MaK2p*F!; z4r}UaH>**5CUH*)5n-6OS;26+-@4LaIBLt(0)fD=ES-77l@mGZesDHy-ye}I1U4%h z2ywMLlwV*tBrr^m9~qV00<{gNxOvG>kxA=XvLB zC~$tU2TtmpUKS4v=EX#&6hkq@CS zNd5N(L8`a$uFS`JXx0QwJjvi2&SDVS1sRKWhru&{6kKpm8CYMu^Zl4W{1TiK41)lw zdhsq22@Z|uCCnqzFj@EZ0Y3DCe9WRY4_1CRLR4hF!ssu)`7!k(0i*7K$wRc#ucm z6N)W<9=_^@bur<~(!JA54UG~Y!9)8L`!s53hdTwp0}P&Zk~3iNJT@XJiWNRPt}zgg z8xC2p)?5vNBjx_NRJ>8M-3dXny2t>bdzaEoBfi=XZ6FZX|5BJk3OTzwq>nD?KpL%c z@1B8+GN?e#j=aw{^1V;8DbHaqo_8LXbmtUMcK$GRV&wPE?Q%>?x8HHX*AjASS$VTQp?g1Z!;cIT$sNyn;l9=@cf(~|4Rh^ zA3n}n|I;|vH3VDKKE1wzXkcyVtoVdylYd;otFP*!^*2@%K2WgF+~H*UOKr`9uC;E9 zwTyF6w%C)0G_TW)ks{o#|J)rGIf?0)&mM<(4lCNl(A-k>AFfL^rwv0mzB?2>!Tfvl zD%LjB!MF!|UF6Z1>366?=RrBunAc+l1)Yq*IgCs9lYAXI$XEk8@7#hj%;M8ChS6PJ z53I7-496P7pG+(ujT%IxgMI&1w0V%B*!CwO8Ln~VLnkEW{8Rvl3NPSwJNY+E?%vVA zUQ+n}!GZ=;?X8@F_4L{Yezo8SjP7qw1l96)OlPJXGwEx84=efF)ozi+Mn3AxEF1hU zhh{aL+_wbH2A~;%7HlQdzvLjHWYsi1hVB|3OBxHg0{wEfjkEIaG{KZ#AS?&18;` z5(%7)z1YSpcnHzWtj^*_%u-ORhtlLdb)Bdv``q{Ozm-HOh+{5W$|4LDD35hB&msv< zpG2pD3J{C>Term!pQc6d3$CQFhdy8~McA63ay?21RIN+3poc5CW2jdXuT-|L#jB#I zA8nAG1J#ApLE>q%RJkxFE3?~RLG52pUYCvnSZjML|J>&1+tEqQ&36j9O4Q`3At(yb z-;bE{BtTG=zWZY%e`$_t%{pfM5x9tP#LC^G{>wv4`>;amj4Z95C004(Cww{0tI-@CZN=f4&Q-x=%Fj(4A zTRF_z1_qja`cmir4m|^Y($w0pAWmYiK_V=_UWiT$P_}k|Jr6;l%Ksp)?z5RiMTSre z8IdD9YhIJ;14SKtd*bxn0F3k=loH^>vwgz39SL)dxDB zhmX`|Z|q%WL{7F#=H|M1J|Oc`7Y(?^^Gts*= zxy@5|7E4$Gi@0=?hK!`uV?g>Xq9|+S$Z{1#Y#L;_oZc!N>HXX4924eH*B|2G1k*S~&ZzJ?X zo~ATF7;xUWFYnq&X$nsh`0=xu5nmTMbm!qbkh(nr%2Q*d1x)I1fo9OxMlR<3rv>O$ zpr*jb((oxUOtR5R6w`*J9%UT&<6dtjWCWt>is3X4SY}QK7ZWls{{pe@eRH7m%J-@{ z-a?OQg3Qb-Uf%uBHd{Q_%7HZMJ%+Duqo6@12oeFY0Uws#3rld40-A48qEjfr5E`TryJTUjq2?x_Y+S)N zrPN~2!&ooHU&xDRiJ*m?2`J@pwQbF5F-+>br2~uXxvw4dLmmQcZ%|5V*J68LJSZ!^ zh>?BX@s@-W(^GVhWkTSpg&0S$$S<}W`j39Ob{xK7&o9a{J7zdeWRS3Q?#0^oK6@zX z$8URr< zglW(GHZHPkCLqotU_6`sV;|ItQ4AP6npSb9+X9}|2dEC@r_zNA-6Orj-4uXNCZsnp zJY0QK%FHR_gPQT%<3DGZ85-VJBCtvPZLvx~r@S#=B_s^>!x{Wz{jVD{jvJ zGjH0CB1xJq3YMmV%EF#szG;lvh)y*i)53H}udS>!vgyLYYS8B@E*X`RTToXW$ztP1 z{x#w5vH43uXxxD`Vrv$ROB!jok0=|)zR3~0dRU^>BvPg0SpxmCm3Q0ObVqeDTEp(S z&xm2hSxDE4_^v_$Bap%v!;O1yl;ubt+8eCBn~H)5K2nAA=kBTN_Y3nrPk(3`O{Sh@ z;exos7xVZ%&HbSCjNB-z;s7`O8Eq9v3w_A(9x*H4>lajyCmxI;4RbWf>YMl9PahUk zpjBT2*ycW_P*;l!uPvfUwy%hj;UF#*HmG7O^iNxw1Fz6cc1f^1kI+dO-h3nsFbm#= z^w-Oc7=olAv4@^ZD`OvzZBh`{gbSjJ6=`XNz){F;1xc4ntYPUl)Hub&N!Q?X4cPy(A|~<*HA*Z`3y4k{kar&-`Xo zoGBG@+FHH*z9lPD+t>ZpR(jI-6L--fm|xJex3y(h9a$AnZ?tBQen}CKex~nTpNy86 z0~|7<$UKLB>jYGUBr@87RcPXRT?-NrcaUi}=3(1r>EbgxDG1q4h*^)8pihG`E4)#Q zH&#(SX#5Mg!nG`Svr;Wm5|Rv+q|;;QXX;#aOVDdU$?zoU z>^NXeewl7`$=tWUyPF5P72tJA*#lbl@81B8G?4a>R;5RdzR;g;La|bJRkh^189%(U zWE&dZiRt!LtJ6Ca%%Xca(R~MSJ!}5|)5DpEL%r{Dd_;_zYzJj;&^BB$gmNk-G$>87 z#z+oXM@YobW|~CZ$8zkd>q_bN%C)hv)hIe7Dcn z=ly!W75k^aYcIj1LmyMEcBU2oK9$?!--bYqRgs_=pnAgXK_4?*{!wF4UXTCn(o$rJ zwIjyHWxdsysv+_`lDOJ>j zc@W%Tddvi0Czg=5v;5^Z+s=2!jNbv8bi(=r+zoSPs1 zuSw2ycT3CkwcgQL-dxTLpX2{IFgnq{9YE_3O5v#7@L4(azdX%}%NRH0#gdYe5!YMP ziH>LWd%jKa54TN;O_Kt2(Afb8v#^0GUPGUEigUPu69Hc%0kv-5kt+ZLEte}AO;s+V z-J+Tk_12Mme0=V^Uk7)Qot>SLfi@4Qo{yJ89bV}R=sYHv-Zt*Wp5O=+^5}AX9WZd2OHx*p8HEY9`B+9SJ*Mm zfoqv}aAOK=@?l|(9EP75z|%{|6|;fWR#R_+AA6;4d1nIlVOZ(X|X=Po-7 z#f#q3);2pg2iR?9Ztmvkx3LE<$2!EhUe^4Btf1H0$!U7&=5X8bMPHv_X@&AHs0jrH zVxz0(W9es9yo@zPW6g@}Z8Nq$8hN~Sm!5<(IY>`}c%Ae?yP7OSAbhEO-J^={3swJG z{`hc(0Ww7h&S$x(287j}NCPBz2q_{0kf3e42gUO2O5f!}vW&1s+_?hV)XUp@R;{DG zy&eLnk>Ox{44cX1A5acK;~VIJIT8KG>wOI%C}BU~-nryn`uWqRyjPDibzlI`&-X6$ zCV0qqCZ65!Te&k@i~dzvy9#^!(IlOmn5Zwu7)Y9!n4B4kgYFpp@GsbZ0lb*pU^S1f z1w%^m(<-tb8yceaG4s*4tA)l+o|Gb-IL8&yc406QyL|7RS=trD(3H_zcNvt;8)HvL zNBN5lX2cl8Akmfl4YHuDpdY$sHV#k(^Nw$qhmH7oQ(`8@;2Orm=80$d3SnUjb?F%R z+JbRhp#qu{{i+6F(>Q1P@NFm#RRDj@s`W_8Kq~|B2#H#qx@UkfLEeRV7}0Wwegdw} zVYs=uu~;nRY7{$_HhPA@5=mQ*`~nRxhqwCn8gD#mTr=>aC+Z&)kh^8{af4;p#dXKLAJ!1?dth4|n z)Kpc0*4yf!k-ZmNASx=#-Tj(n9$77XlYBzm-Mau?4I+{Bur?jR2^t*mK41>5d)CWb z!epb?*zRWd`4wHit~TqaEsp@FV2H42ejdJUNTpUFa|1OP`}9G`u%(ricx@IfmqAHI z0$FJ7ujpQb*+DJGGO7#vATxWz7*S9I1BX*ORuZ_FnbmuxM6%m%mbl2t$|7C_ucppH zcB0y@=u-%bKq}d>N=cO+jh^LL7JFE(I0s!_vy+k@y$!v3^(q4MlvPyVtc#Zdy&LsS z;D@xiiw&8uf`Nb8xpQaH?r(SxZTSE+yhmzV0d(~Rh95jAU3&*4y@5gZjdtwm3W_zR zi@}gGlY(WI$z(z*6EOp~gA^*n*`c zTn-)jq;zPdhp~si;3R?DjO;NmOwIbs6*reb|7WWgmBR=RH>^5_b8XKa;n_K?0Km+@ zRrYtsrKdZ&x+>8MYis)<)c5yiD%aH1?2eonG4coGphKYyK<9PS0hqdf|Nea4o)2MQ zXi%hcC2f)B7>#1SSxG4ymx<)#nCb&(?GhbOCJBAR|AQc7jwX zFAt%sSyNJV-0x*X^ZWa2VPyl1$w_?V;TsT;ce(>_|8L%SQE>vIo)l%iXf126#l@|? zy=fhzhhj;h!?563=O;fMMmUG&h7FVay7UEWZlf0!(hyM8Vc^}9lQCdhPK>rko;sCi z10EIYx!L0L^eX@g>kE=oV7jJIDDv9QrdO8n++2pp%9UMCc|?&F{(n4{r0D8Gkl6J# zt7CK%t_Btsc+e>}J0dI$j+%HkdnOa1i=-5*7=p)pwY)ChDMi9o{j0VG_7oLmWzqQ( zi{c)C1I3vk=Oqh(9Aw?W%&qH4mU&6v`Lt19Pn5U;qFB literal 0 HcmV?d00001 diff --git a/sphere_drag.m b/sphere_drag.m index e75818f..204b485 100644 --- a/sphere_drag.m +++ b/sphere_drag.m @@ -1,10 +1,26 @@ function [Cd_out] = sphere_drag(Re_in,spline_type) % interpolation for drag coeff based upon Reynolds number - % output is + % output is % Cd: drag coeff - % input is + % input is % Re_in: Reynolds number - % and + % and % spline_type: 'linear', 'cubic', or 'spline'/'pchip' - - \ No newline at end of file + + %function [Cd_out]=sphere_drag(Re_in,spline_type) + data = [2 0.52 + 5.8 0.52 + 16.8 0.52 + 27.2 0.5 + 29.9 0.49 + 33.9 0.44 + 36.3 0.18 + 40 0.074 + 46 0.067 + 60 0.08 + 100 0.12 + 200 0.16 + 400 0.19]; + + Cd_out = interp1(data(:,1),data(:,2),Re_in,spline_type); +end