Advanced Search

salim-barzani

260 Reputation

3 Badges

0 years, 220 days
Contact salim-barzani

MaplePrimes Activity


These are questions asked by salim-barzani

How factoring two term after substituti...

salim-barzani 260 Maple 2024
pde differential-equations coefficients
December 27 2024
1 6

I already get the same results, but there's something about the factoring process that I encountered for the first time in this ODE. In the paper, it says that G′ satisfies a certain condition, but I’m not sure exactly what that means. Did the author use it for substitution, or did they change (m+F) into another variable and then solve? I’m not exactly sure what approach was taken. Does anyone have any idea or insight into this?

> 

restart
> 

with(PDEtools)
> 

with(LinearAlgebra)
> 

NULL
> 

with(SolveTools)
> 

undeclare(prime)

`There is no more prime differentiation variable; all derivatives will be displayed as indexed functions`

 (1)
> 

declare(Omega(x, t)); declare(U(xi)); declare(u(x, y, z, t)); declare(Q(xi))

Omega(x, t)*`will now be displayed as`*Omega

  

U(xi)*`will now be displayed as`*U

  

u(x, y, z, t)*`will now be displayed as`*u

  

Q(xi)*`will now be displayed as`*Q

 (2)
> 

tr := {t = tau, x = (-ZETA*c[3]-tau*c[4]-`Υ`*c[2]+xi)/c[1], y = `Υ`, z = ZETA, u(x, y, z, t) = U(xi)}

{t = tau, x = (-Zeta*c[3]-tau*c[4]-`Υ`*c[2]+xi)/c[1], y = `Υ`, z = Zeta, u(x, y, z, t) = U(xi)}

 (3)
> 

pde1 := diff(u(x, y, z, t), `$`(x, 3), z)-4*(diff(u(x, y, z, t), x, t))+4*(diff(u(x, y, z, t), x))*(diff(u(x, y, z, t), x, z))+2*(diff(u(x, y, z, t), `$`(x, 2)))*(diff(u(x, y, z, t), z))+3*(diff(u(x, y, z, t), `$`(y, 2))) = 0

diff(diff(diff(diff(u(x, y, z, t), x), x), x), z)-4*(diff(diff(u(x, y, z, t), t), x))+4*(diff(u(x, y, z, t), x))*(diff(diff(u(x, y, z, t), x), z))+2*(diff(diff(u(x, y, z, t), x), x))*(diff(u(x, y, z, t), z))+3*(diff(diff(u(x, y, z, t), y), y)) = 0

 (4)
> 

``
> 

L1 := PDEtools:-dchange(tr, pde1, [xi, `Υ`, ZETA, tau, U])

(diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*c[1]^3*c[3]-4*(diff(diff(U(xi), xi), xi))*c[4]*c[1]+6*(diff(U(xi), xi))*c[1]^2*(diff(diff(U(xi), xi), xi))*c[3]+3*(diff(diff(U(xi), xi), xi))*c[2]^2 = 0

 (5)
> 

map(int, (diff(diff(diff(diff(U(xi), xi), xi), xi), xi))*c[1]^3*c[3]-4*(diff(diff(U(xi), xi), xi))*c[4]*c[1]+6*(diff(U(xi), xi))*c[1]^2*(diff(diff(U(xi), xi), xi))*c[3]+3*(diff(diff(U(xi), xi), xi))*c[2]^2 = 0, xi)

c[1]^3*c[3]*(diff(diff(diff(U(xi), xi), xi), xi))+3*c[2]^2*(diff(U(xi), xi))-4*c[4]*c[1]*(diff(U(xi), xi))+3*c[1]^2*c[3]*(diff(U(xi), xi))^2 = 0

 (6)
> 

ode := %

c[1]^3*c[3]*(diff(diff(diff(U(xi), xi), xi), xi))+3*c[2]^2*(diff(U(xi), xi))-4*c[4]*c[1]*(diff(U(xi), xi))+3*c[1]^2*c[3]*(diff(U(xi), xi))^2 = 0

 (7)
> 

F := sum(a[i]*(m+1/(diff(G(xi), xi)))^i, i = -1 .. 1)

a[-1]/(m+1/(diff(G(xi), xi)))+a[0]+a[1]*(m+1/(diff(G(xi), xi)))

 (8)
> 

D1 := diff(F, xi)

a[-1]*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-a[1]*(diff(diff(G(xi), xi), xi))/(diff(G(xi), xi))^2

 (9)
> 

S := diff(G(xi), `$`(xi, 2)) = -(2*m*mu+lambda)*(diff(G(xi), xi))-mu

diff(diff(G(xi), xi), xi) = -(2*m*mu+lambda)*(diff(G(xi), xi))-mu

 (10)
> 

E1 := subs(S, D1)

a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2

 (11)
> 

D2 := diff(E1, xi)

2*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^4)-2*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^3)-a[-1]*(2*m*mu+lambda)*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)+2*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)*(diff(diff(G(xi), xi), xi))/(diff(G(xi), xi))^3+a[1]*(2*m*mu+lambda)*(diff(diff(G(xi), xi), xi))/(diff(G(xi), xi))^2

 (12)
> 

E2 := subs(S, D2)

2*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^4)-2*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^3)-a[-1]*(2*m*mu+lambda)*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)+2*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2/(diff(G(xi), xi))^3+a[1]*(2*m*mu+lambda)*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2

 (13)
> 

D3 := diff(E2, xi)

6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^4*(diff(G(xi), xi))^6)-12*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^5)-6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)*(2*m*mu+lambda)*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^4)+6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^4)+6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)*(2*m*mu+lambda)*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^3)+a[-1]*(2*m*mu+lambda)^2*(diff(diff(G(xi), xi), xi))/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-6*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(diff(diff(G(xi), xi), xi))/(diff(G(xi), xi))^4-6*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)*(2*m*mu+lambda)*(diff(diff(G(xi), xi), xi))/(diff(G(xi), xi))^3-a[1]*(2*m*mu+lambda)^2*(diff(diff(G(xi), xi), xi))/(diff(G(xi), xi))^2

 (14)
> 

E3 := subs(S, D3)

6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/((m+1/(diff(G(xi), xi)))^4*(diff(G(xi), xi))^6)-12*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^5)-6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(2*m*mu+lambda)/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^4)+6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^4)+6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(2*m*mu+lambda)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^3)+a[-1]*(2*m*mu+lambda)^2*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-6*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/(diff(G(xi), xi))^4-6*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(2*m*mu+lambda)/(diff(G(xi), xi))^3-a[1]*(2*m*mu+lambda)^2*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2

 (15)
> 

NULL
> 

NULL
> 

K := U(xi) = F

U(xi) = a[-1]/(m+1/(diff(G(xi), xi)))+a[0]+a[1]*(m+1/(diff(G(xi), xi)))

 (16)
> 

K1 := diff(U(xi), xi) = E1

diff(U(xi), xi) = a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2

 (17)
> 

K2 := diff(U(xi), `$`(xi, 2)) = E2

diff(diff(U(xi), xi), xi) = 2*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^4)-2*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^3)-a[-1]*(2*m*mu+lambda)*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)+2*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2/(diff(G(xi), xi))^3+a[1]*(2*m*mu+lambda)*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2

 (18)
> 

K3 := diff(U(xi), `$`(xi, 3)) = E3

diff(diff(diff(U(xi), xi), xi), xi) = 6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/((m+1/(diff(G(xi), xi)))^4*(diff(G(xi), xi))^6)-12*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^5)-6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(2*m*mu+lambda)/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^4)+6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^4)+6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(2*m*mu+lambda)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^3)+a[-1]*(2*m*mu+lambda)^2*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-6*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/(diff(G(xi), xi))^4-6*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(2*m*mu+lambda)/(diff(G(xi), xi))^3-a[1]*(2*m*mu+lambda)^2*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2

 (19)
> 

NULL
> 

L := eval(ode, {K, K1, K2, K3})

c[1]^3*c[3]*(6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/((m+1/(diff(G(xi), xi)))^4*(diff(G(xi), xi))^6)-12*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^5)-6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(2*m*mu+lambda)/((m+1/(diff(G(xi), xi)))^3*(diff(G(xi), xi))^4)+6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^4)+6*a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(2*m*mu+lambda)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^3)+a[-1]*(2*m*mu+lambda)^2*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-6*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^3/(diff(G(xi), xi))^4-6*a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)^2*(2*m*mu+lambda)/(diff(G(xi), xi))^3-a[1]*(2*m*mu+lambda)^2*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2)+3*c[2]^2*(a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2)-4*c[4]*c[1]*(a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2)+3*c[1]^2*c[3]*(a[-1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/((m+1/(diff(G(xi), xi)))^2*(diff(G(xi), xi))^2)-a[1]*(-(2*m*mu+lambda)*(diff(G(xi), xi))-mu)/(diff(G(xi), xi))^2)^2 = 0

 (20)
> 

"collect(L,(m+1/(diff(G(xi),xi))))^( )"

Error, (in collect) cannot collect m+1/diff(G(xi),xi)

  
 

NULL

Download factoring.mw

Why in two platform for ODE equation i h...

salim-barzani 260 Maple 2024
ode dsolve differential-equations
December 24 2024
1 16

P1.mw

what is pronlem why i can't get solution...

salim-barzani 260 Maple 2024
syntax 2dmath
December 23 2024
1 9

i already use this method for a lot of equation but this time something not normal hapening what is problem?

``

> 

restart
> 

with(PDEtools)
> 

with(LinearAlgebra)
> 

with(Physics)
> 

with(SolveTools)
> 

``
> 

eq0 := -4*alpha*k^2*m^2*n^2*A[0]^2+4*beta*k*m*n^2*A[0]^3-4*gamma*k*m*n^2*A[0]^3+4*delta^2*m*n^2*A[0]^2-4*n^2*sigma*A[0]^4-4*m*n^2*w*A[0]^2 = 0
> 

eq1 := -8*alpha*k^2*m^2*n^2*A[0]*A[1]+12*beta*k*m*n^2*A[0]^2*A[1]-12*gamma*k*m*n^2*A[0]^2*A[1]+8*delta^2*m*n^2*A[0]*A[1]-16*n^2*sigma*A[0]^3*A[1]+2*a*alpha*m*n*A[0]*A[1]-8*m*n^2*w*A[0]*A[1] = 0
> 

eq2 := -4*alpha*k^2*m^2*n^2*A[1]^2+12*beta*k*m*n^2*A[0]*A[1]^2-12*gamma*k*m*n^2*A[0]*A[1]^2+4*delta^2*m*n^2*A[1]^2-24*n^2*sigma*A[0]^2*A[1]^2+a*alpha*m^2*A[1]^2+3*alpha*b*m*n*A[0]*A[1]-4*m*n^2*w*A[1]^2 = 0
> 

eq3 := 4*beta*k*m*n^2*A[1]^3-4*gamma*k*m*n^2*A[1]^3-16*n^2*sigma*A[0]*A[1]^3+alpha*b*m^2*A[1]^2+alpha*b*m*n*A[1]^2+4*alpha*c*m*n*A[0]*A[1] = 0
> 

eq4 := -4*n^2*sigma*A[1]^4+alpha*c*m^2*A[1]^2+2*alpha*c*m*n*A[1]^2 = 0
> 

C := solve({eq0, eq1, eq2, eq3, eq4}, {a, b, c, `__ `*A[0]})

Warning, solving for expressions other than names or functions is not recommended.

  

(1)
 

NULL

Download problem.mw

How remove error when changing delta par...

salim-barzani 260 Maple 2024
finance spam numeric
December 17 2024
0 14

This is my first time working with plotting data from a matrix. However, with the help of a friends on MaplePrimes, I learned how to plot the data in both Maple and MATLAB. Despite this, I am having trouble with visualization. When I change the delta value, my function experiences vibrations or noise, which is clearly visible in the plot. But when I change delta, I encounter errors with my matrix data. How can I fix this problem? and there is any way for get better visualization by Explore ? also How show this vibration or noise in 2D?

> 

restart;
> 

randomize():
> 

local gamma;

gamma

 (1)
> 

currentdir(kernelopts(':-homedir'))
> 

NULL
> 

T3 := (B[1]*(tanh(2*n^2*(delta^2-w)*k*t/((k*n-1)*(k*n+1))+x)-1))^(1/(2*n))*exp(I*(-k*x+w*t+delta*W(t)-delta^2*t))

(B[1]*(tanh(2*n^2*(delta^2-w)*k*t/((k*n-1)*(k*n+1))+x)-1))^((1/2)/n)*exp(I*(-k*x+w*t+delta*W(t)-delta^2*t))

 (2)
> 

NULL
> 

params := {B[1]=1,n=2,delta=1,w=1,k=3 };

{delta = 1, k = 3, n = 2, w = 1, B[1] = 1}

 (3)
> 

NULL

insert numerical values

> 

solnum :=subs(params, T3);

(tanh(x)-1)^(1/4)*exp(I*(-3*x+W(t)))

 (4)
> 

CodeGeneration['Matlab']('(tanh(x)-1)^(1/4)*exp(I*(-3*x+W(t)))')

Warning, the function names {W} are not recognized in the target language

  

cg = ((tanh(x) - 0.1e1) ^ (0.1e1 / 0.4e1)) * exp(i * (-0.3e1 * x + W(t)));

  
> 

N := 100:

use Finance in:
  Wiener := WienerProcess():
  P := PathPlot(Wiener(t), t = 0..10, timesteps = N, replications = 1):
end use:
> 

W__points := plottools:-getdata(P)[1, -1]:
t_grid := convert(W__points[..,1], list):
x_grid := [seq(-2..2, 4/N)]:
> 

T, X := map(mul, [selectremove(has, [op(expand(solnum))], t)])[]:

ST := unapply(eval(T, W(t)=w), w)~(W__points[.., 2]):
SX := evalf(unapply(X, x)~(x_grid)):

STX := Matrix(N$2, (it, ix) -> ST[it]*SX[ix]);

_rtable[36893490640185799852]

 (5)
> 

opts := axis[1]=[tickmarks=[seq(k=nprintf("%1.1f", t_grid[k]), k=1..N, 40)]],
        axis[2]=[tickmarks=[seq(k=nprintf("%1.1f", x_grid[k]), k=1..N, 40)]],
        style=surface:

DocumentTools:-Tabulate(
  [
    plots:-matrixplot(Re~(STX), opts),
    plots:-matrixplot(Im~(STX), opts),
plots:-matrixplot(abs~(STX), opts)
  ]
  , width=60
)

"Tabulate"

 (6)
> 

MatlabFile := cat(currentdir(), "/ST2.txt"); ExportMatrix(MatlabFile, STX, target = MATLAB, format = rectangular, mode = ascii, format = entries)

421796

 (7)

NULL

Download data-analysis.mw

Transform data in maple to Matlab or mak...

salim-barzani 260 Maple 2024
matlab excel export
December 16 2024
1 7

I have a matrix for data analysis that I want to plot. Ideally, I would like to use Maple, but I’m struggling to create a well-designed plot suitable for submission to journals. Because of this, I’m considering transferring the data to Excel or constructing a 3D graph using MATLAB.

My question is: how can I transfer this data to Excel? The data is currently saved as a Notepad file, but I’m unsure how to convert it into an Excel format. I will upload a figure to show the data structure.

also in last runig program give me error which is (Error, (in ExportMatrix) permission denied

Thank you in advance for any help!

> 

restart;
> 

randomize():
> 

local gamma;

gamma

 (1)
 
> 

T3 := (B[1]*(tanh(2*n^2*(delta^2-w)*k*t/((k*n-1)*(k*n+1))+x)-1))^(1/(2*n))*exp(I*(-k*x+w*t+delta*W(t)-delta^2*t))

(B[1]*(tanh(2*n^2*(delta^2-w)*k*t/((k*n-1)*(k*n+1))+x)-1))^((1/2)/n)*exp(I*(-k*x+w*t+delta*W(t)-delta^2*t))

 (2)
> 

``
> 

params := {B[1]=1,n=2,delta=1,w=1,k=3 };

{delta = 1, k = 3, n = 2, w = 1, B[1] = 1}

 (3)
> 

``

insert numerical values

> 

solnum :=subs(params, T3);

(tanh(x)-1)^(1/4)*exp(I*(-3*x+W(t)))

 (4)
> 

CodeGeneration['Matlab']('(tanh(x)-1)^(1/4)*exp(I*(-3*x+W(t)))')

Warning, the function names {W} are not recognized in the target language

  

cg = ((tanh(x) - 0.1e1) ^ (0.1e1 / 0.4e1)) * exp(i * (-0.3e1 * x + W(t)));

  
> 

N := 100:

use Finance in:
  Wiener := WienerProcess():
  P := PathPlot(Wiener(t), t = 0..10, timesteps = N, replications = 1):
end use:
> 

W__points := plottools:-getdata(P)[1, -1]:
t_grid := convert(W__points[..,1], list):
x_grid := [seq(-2..2, 4/N)]:
> 

T, X := map(mul, [selectremove(has, [op(expand(solnum))], t)])[]:

ST := unapply(eval(T, W(t)=w), w)~(W__points[.., 2]):
SX := evalf(unapply(X, x)~(x_grid)):

STX := Matrix(N$2, (it, ix) -> ST[it]*SX[ix]);

_rtable[36893489786521178348]

 (5)
> 

opts := axis[1]=[tickmarks=[seq(k=nprintf("%1.1f", t_grid[k]), k=1..N, 40)]],
        axis[2]=[tickmarks=[seq(k=nprintf("%1.1f", x_grid[k]), k=1..N, 40)]],
        style=surface:

DocumentTools:-Tabulate(
  [
    plots:-matrixplot(Re~(STX), opts),
    plots:-matrixplot(Im~(STX), opts),
plots:-matrixplot(abs~(STX), opts)
  ]
  , width=60
)

"Tabulate"

 (6)
> 

MatlabFile := cat(currentdir(), "/ST2.txt"); ExportMatrix(MatlabFile, STX, target = MATLAB, format = rectangular, mode = ascii, format = entries)

Error, (in ExportMatrix) permission denied

  
 

 

Download data-analysis.mw

5 6 7 8 9 10 11 Last Page 7 of 20