plot2d discrete
- Subject: plot2d discrete
- From: Mahery Raharinjatovo
- Date: Wed, 27 Aug 2008 13:37:39 +0400
Hello,
I need help
I made the program below
I would want to plot the curve define by points of coordinates (t,v)
chute():=block([g,k,r,rho,rhoh,vm,V,v],
g:9.8,
k:9.33/100,
rho:2400,
rhoh:950,
r:1.5/1000,
V:(4/3)*3.14*r*r*r,
vm:(g*rho*V/k)*(1-rhoh/rho),
dt:1.5/10000,
a:-k/(rho*V),
b:g*(1-(rhoh/rho)),
v:0,t:0,
while v<0.999*vm do (print(v,"?",t),dv:(a*v+b)*dt,v:v+dv,t:t+dt)
)$
(%i20) chute();
0 ? 0
8.8812499999999992E-4 ? 1.4999999999999999E-4
0.0014097334461253 ? 2.9999999999999997E-4
0.0017160815064425 ? 4.4999999999999999E-4
0.0018960040800842 ? 5.9999999999999995E-4
0.0020016751669921 ? 7.4999999999999991E-4
0.00206373729951 ? 8.9999999999999987E-4
0.0021001872727748 ? 0.00105
0.0021215948605761 ? 0.0012
0.0021341678414158 ? 0.00135
0.0021415521319143 ? 0.0015
0.002145889030766 ? 0.00165
0.0021484361531543 ? 0.0018