Advanced Search

Question: An easier way to plot data

Question:An easier way to plot data

raggapuffin 89 Maple
plot graph
April 15 2010
0
I have some code for evaluating a velocity given different variables: > h := 50; e := 5/6; > d := .218*100; Omega := 6*Pi*(1/100); > a := -h/e; b := h/(1-e); > Q := piecewise(x <><><><><><><><><><><><><><><><><><><> k := piecewise(0 <><><><> F := 0; > tau[1] := 1/(30*d); tau[2] := (1/4)*tau[1]; > PDE := diff(u(x, t), t) = (k+F)*(diff(u(x, t), x))/Omega+d*(diff(u(x, t), x, x)); > with(plots); a := 1/70; > IBC := {u(0, t) = 0, u(5, t) = 0, u(x, 0) = exp(-(x-2-e)^2/(2*a^2))/((2*Pi)^(1/2)*a)}; > pds := pdsolve(PDE, IBC, numeric, time = t, range = 0 .. 5, timestep = 1/10000, spacestep = 1/5000); > R := pds:-value(output = listprocedure); uu := subs(R, u(x, t)); > S := evalf(Int(x*uu(x, tau[1]+tau[2]), x = 0 .. 5, method = _Dexp)); B := evalf(Int(uu(x, tau[1]+tau[2]), x = 0 .. 5, method = _Dexp)); > V[1] := 10*(S/B-2-e)/(tau[1]+tau[2]); but i'd like to be able to plot the velocity on a graph against a range of values of say, tau_1 or any parameter. At the moment I'm just restarting this code and changing the values of the parameters to obtain datapoints which im then plotting with listplot, which is very tedious, and takes ages. Does anyone have an easier way of doing this? Ryan
Please Wait...