Advanced Search

permanoon123

Mr. Erfan Permanoon

165 Reputation

5 Badges

5 years, 188 days
Contact Erfan Permanoon

MaplePrimes Activity


These are questions asked by permanoon123

how can i plot of uu10000...

permanoon123 165 Maple 2022
complex complexplot plotting
March 20 2023
3 3
> 

restart
> 

with(plottools)
> 

with(plots)
> 

with(CurveFitting)
> 

Digits := 10
> 

with(GaussInt)
> 

w := GInearest(0+I)

I

 (1)
> 

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
> 

v := .7

.7

 (2)
> 

t := 10000

10000

 (3)
> 

i = sqrt(-1)

i = I

 (4)
> 

"k(n) := evalf((2 *Pi*n)/(L))"

proc (n) options operator, arrow, function_assign; evalf(2*Pi*n/L) end proc

 (5)
> 

f(n) = (int(f(t)*exp(-w*k(n)*x), x = 0 .. L))/L
> 

NULL
> 

"C(x, t) := (∑) exp(-v* t*k(n)- D1 *t*(k(n)^())^(2)- epsilon *t*(k(n))^(4)) *f(n)* exp(w*k(n)* x)"

proc (x, t) options operator, arrow, function_assign; sum(exp(-v*t*k(n)-D1*t*k(n)^2-varepsilon*t*k(n)^4)*f(n)*exp(w*k(n)*x), n = 1 .. 10) end proc

 (6)
> 

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

[0.8582270020e-37+0.7071768085e-46*I, 0.8542115620e-37-0.8289193755e-38*I, 0.8422030662e-37-0.1650065052e-37*I, 0.8223146210e-37-0.2455737186e-37*I, 0.7947335229e-37-0.3238382907e-37*I, 0.7597194496e-37-0.3990668010e-37*I, 0.7176019535e-37-0.4705546418e-37*I, 0.6687772885e-37-0.5376326775e-37*I, 0.6137045955e-37-0.5996735584e-37*I, 0.5529014934e-37-0.6560976155e-37*I, 0.4869391210e-37-0.7063782859e-37*I, 0.4164366719e-37-0.7500470239e-37*I, 0.3420554942e-37-0.7866976371e-37*I, 0.2644927960e-37-0.8159900244e-37*I, 0.1844750367e-37-0.8376532722e-37*I, 0.1027510610e-37-0.8514880864e-37*I, 0.2008503850e-38-0.8573685436e-37*I, -0.6275071112e-38-0.8552431450e-37*I, -0.1449829443e-37-0.8451351706e-37*I, -0.2258447115e-37-0.8271423331e-37*I, -0.3045824843e-37-0.8014357420e-37*I, -0.3804631338e-37-0.7682581907e-37*I, -0.4527807059e-37-0.7279217843e-37*I, -0.5208629116e-37-0.6808049422e-37*I, -0.5840773014e-37-0.6273487939e-37*I, -0.6418370476e-37-0.5680530177e-37*I, -0.6936062957e-37-0.5034711554e-37*I, -0.7389050434e-37-0.4342054446e-37*I, -0.7773134955e-37-0.3609012248e-37*I, -0.8084758661e-37-0.2842409660e-37*I, -0.8321035975e-37-0.2049379645e-37*I, -0.8479779638e-37-0.1237297776e-37*I, -0.8559520434e-37-0.4137144557e-38*I, -0.8559520435e-37+0.4137144402e-38*I, -0.8479779639e-37+0.1237297769e-37*I, -0.8321035979e-37+0.2049379630e-37*I, -0.8084758664e-37+0.2842409653e-37*I, -0.7773134957e-37+0.3609012242e-37*I, -0.7389050442e-37+0.4342054433e-37*I, -0.6936062961e-37+0.5034711548e-37*I, -0.6418370481e-37+0.5680530171e-37*I, -0.5840773026e-37+0.6273487928e-37*I, -0.5208629121e-37+0.6808049418e-37*I, -0.4527807065e-37+0.7279217840e-37*I, -0.3804631352e-37+0.7682581900e-37*I, -0.3045824850e-37+0.8014357418e-37*I, -0.2258447130e-37+0.8271423327e-37*I, -0.1449829450e-37+0.8451351705e-37*I, -0.6275071183e-38+0.8552431450e-37*I, 0.2008503694e-38+0.8573685437e-37*I, 0.1027510603e-37+0.8514880865e-37*I, 0.1844750360e-37+0.8376532724e-37*I, 0.2644927944e-37+0.8159900249e-37*I, 0.3420554936e-37+0.7866976374e-37*I, 0.4164366713e-37+0.7500470242e-37*I, 0.4869391201e-37+0.7063782866e-37*I, 0.5529014927e-37+0.6560976160e-37*I, 0.6137045944e-37+0.5996735595e-37*I, 0.6687772877e-37+0.5376326784e-37*I, 0.7176019529e-37+0.4705546426e-37*I, 0.7597194490e-37+0.3990668022e-37*I, 0.7947335225e-37+0.3238382916e-37*I, 0.8223146207e-37+0.2455737195e-37*I, 0.8422030659e-37+0.1650065065e-37*I, 0.8542115619e-37+0.8289193856e-38*I, 0.8582270020e-37]

 (7)
 

``

Download 0_one.mw

I can't run plot of time series!...

permanoon123 165 Maple 2022
syntax user-error
March 09 2023
1 0
> 

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

Could you help me to solve Error of Inte...

permanoon123 165 Maple 2022
integration syntax 2dmath
March 06 2023
1 2
> 

restart
> 

with(plottools)
> 

with(plots)
> 

with(CurveFitting)
> 

Digits := 100
> 

"g(t):=10*(e)^((-(t-4000)^(2))/(1300000))+6*(e)^((-(t-6900)^(2))/(1400000))"

proc (t) options operator, arrow, function_assign; 10*exp(-(1/1300000)*(t-4000)^2)+6*exp(-(1/1400000)*(t-6900)^2) end proc

 (1)
> 

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

 

````

> 

v := .7
> 

disp := 15
> 

PDE := diff(C(x, t), t) = -v*(diff(C(x, t), x))+disp*(diff(C(x, t), x, x))
> 

IBC := C(x, 0) = 0, C(0, t) = g(t), (D[1](C))(10000, t) = 0
> 

pds := pdsolve(PDE, [IBC], time = t, range = 0 .. 10000, timestep = 10, numeric, spacestep = 10)

_m2712135358688

 (2)
> 

k := pds:-plot(x = 6500, t = 0 .. 20000, numpoints = 600); plots[display]({k})

 

NULL

> 

"f(t):= unapply(Spline( getdata(k)[3][..,1],getdata(k)[3][..,2],'t',degree=2),t): "
> 

D1 := 15
> 

E := 20000
> 

L := 6500
> 

n := 200
> 

lambda = `√`(n*Pi/L, (1/2)*D1+`√`((1/4)*D1^2+E(`nπ`/L)^2))
> 

"p(t):=(∫)[0]^(L)f(t) (e(-lambda^2 t) )^dt :"

Error, Got internal error in Typesetting:-Parse : "invalid subscript selector"

"p(t):=(∫)[0]^Lf(t) e(-lambda^2 t) dt :"

  
> 

C(x, t) = sum((2*sin(`nπx`/L)*exp(1)/L*sin(`nπx`/L))*p(t), n = 1 .. 500)

Error, (in sum) summation variable previously assigned, second argument evaluates to 200 = 1 .. 500

  
> 

uu20000 := [seq(evalf(C(20000-i, 5000)), i = 0 .. 20000, 100)]
> 

``
> 

xx := [seq(i, i = 0 .. 20000, 100)]
> 

p1 := pds:-plot(t = 33000, numpoints = 150, color = red)
> 

plots[display]({p1})

Download 1.mw

could you please help me to fix the mist...

permanoon123 165 Maple 2022
syntax pde numeric
February 28 2023
1 0
> 

restart
> 

with(plottools)
> 

with(plots)
> 

with(CurveFitting)
> 

Digits := 10

NULL

> 

pde := diff(C(x, t), t) = -.7*(diff(C(x, t), x))+15*(diff(C(x, t), x, x))+70000*(diff(C(x, t), `$`(x, 4))); ic := C(x, 0) = 0; bc := [eval(diff(C(x, t), x), x = 0) = 0, eval(diff(C(x, t), x, x), x = 0) = 0, eval(diff(C(x, t), x, x), x = L) = 0, C(L, t) = 3.65*exp(-(1/2000000)*(t-11150)^2)]; L := 5000; T := 7000; sol := pdsolve([pde, ic, bc], C(x, t), numeric, time = t, range = 0 .. T, spacestep = (1/50)*L); plot(sol(x, 7000), x = 0 .. L, y = 0 .. 4, axes = boxed, gridlines = true, labels = ["x", "C(x,7000)"])

diff(C(x, t), t) = -.7*(diff(C(x, t), x))+15*(diff(diff(C(x, t), x), x))+70000*(diff(diff(diff(diff(C(x, t), x), x), x), x))

  

C(x, 0) = 0

  

[eval(diff(C(x, t), x), {x = 0}) = 0, eval(diff(diff(C(x, t), x), x), {x = 0}) = 0, diff(diff(C(L, t), L), L) = 0, C(L, t) = 3.65*exp(-(1/2000000)*(t-11150)^2)]

  

5000

  

7000

  

Error, invalid input: diff received 5000, which is not valid for its 2nd argument

  

Warning, expecting only range variable x in expression sol(x,7000) to be plotted but found name sol

  

 

NULL

Download QR_method.mw

How to curve fitting in 3dpoint? and ma...

permanoon123 165 Maple 2018
3d curvefitting regression
December 18 2021
2 6

H.M.mw

> 

restart
> 

with(plottools):
> 

with(plots):
> 

with(CurveFitting):
> 

with(Statistics):
> 

Digits := 10:
> 

L := point([0, 0, 0], color = blue, symbol = cross, symbolsize = 50), point([0, 0, 1], color = red, symbol = cross, symbolsize = 50), point([0, 1, 0], color = black, symbol = cross, symbolsize = 50):
> 

display(L, axes = boxed, view = [-1 .. 1, -1 .. 1, -1 .. 1], orientation = [125, 65])

 
> 

``
> 

``
> 

``

``

Download H.M.mw

1 2 3 4 5 Page 1 of 5