Question: I can't run plot of time series!

restart

with(plottools)

with(plots)

with(CurveFitting)

Digits := 10

NULL

"f(t):=7.0*(e)^((-(t-13180)^(2))/(2000000))+4.7*(e)^((-(t-16000)^(2))/(3200000)):"

p1 := plot(f(t), t = 0 .. 20000, color = green); plots[display]({p1})

 

NULL

D1 := 15

epsilon := 200000

L := 6500

n := 200

t := 1000

1000

(1)

lambda := simplify(evalf(n*Pi*sqrt((1/2)*D1+sqrt((1/4)*D1^2+epsilon*(n*Pi/L)^2))/L))

.6928578233

(2)

b := 2*(int(f(t)*sin(m*Pi*x/L), x = 0 .. L))/L

-0.6366197724e-1*(0.1409730543e-28*cos(3.141592654*m)-0.1409730543e-28)/m

(3)

C(x, t) = sum(b*exp^(-lambda^2*t)*sin(m*Pi*x/L), m = 1 .. 2)

C(x, 1000) = 0.1794924675e-29*sin(0.4833219466e-3*x)/exp^(4800519633/10000000)

(4)

uu1000 := [seq(evalf(C(L-i, t)), i = 0 .. 6500, 100)]

[C(6500, 1000), C(6400, 1000), C(6300, 1000), C(6200, 1000), C(6100, 1000), C(6000, 1000), C(5900, 1000), C(5800, 1000), C(5700, 1000), C(5600, 1000), C(5500, 1000), C(5400, 1000), C(5300, 1000), C(5200, 1000), C(5100, 1000), C(5000, 1000), C(4900, 1000), C(4800, 1000), C(4700, 1000), C(4600, 1000), C(4500, 1000), C(4400, 1000), C(4300, 1000), C(4200, 1000), C(4100, 1000), C(4000, 1000), C(3900, 1000), C(3800, 1000), C(3700, 1000), C(3600, 1000), C(3500, 1000), C(3400, 1000), C(3300, 1000), C(3200, 1000), C(3100, 1000), C(3000, 1000), C(2900, 1000), C(2800, 1000), C(2700, 1000), C(2600, 1000), C(2500, 1000), C(2400, 1000), C(2300, 1000), C(2200, 1000), C(2100, 1000), C(2000, 1000), C(1900, 1000), C(1800, 1000), C(1700, 1000), C(1600, 1000), C(1500, 1000), C(1400, 1000), C(1300, 1000), C(1200, 1000), C(1100, 1000), C(1000, 1000), C(900, 1000), C(800, 1000), C(700, 1000), C(600, 1000), C(500, 1000), C(400, 1000), C(300, 1000), C(200, 1000), C(100, 1000), C(0, 1000)]

(5)

``

xx := [seq(k, k = 0 .. 6500, 100)]

NULL

p2 := plot(xx, uu1000, color = cyan)

Error, (in plot) two lists or Vectors of numerical values expected

 

plots[display]({p2})

Error, (in plots:-display) expecting plot structures but received: {p2}

 

NULL

Download easy_way.mw

Please Wait...