I have a equation

((D@@2)(theta))(eta) = -(1/2)*(D(theta))(eta)*(-2*(D(phi))(eta)*beta*epsilon*lambda*D[B]+2*(D(phi))(eta)*beta*epsilon*mu*D[B]+f(eta)*sin(alpha)*beta*nu+2*(D(theta))(eta)*gamma*epsilon*D[t]-2*(D(theta))(eta)*beta*epsilon*D[t]+cos(alpha)*beta*eta*nu)/(beta*sigma)

and a parameters expression

Pr:=nu/sigma; N[b] := epsilon*D[B](mu-lambda)/sigma; N[t] := epsilon*D[t](gamma-beta)/(gamma*sigma); Le := nu/D[B]

How can I seperate common terms and substitute this parameters and got this following expression

((D@@2)(theta))(eta) = -(1/2)*Pr*(D(theta))(eta)*eta*cos(alpha)-(1/2)*Pr*(D(theta))(eta)*sin(alpha)*f(eta)-N[b]*(D(theta))(eta)*(D(phi))(eta)-N[t]*(D(theta))(eta)

How to conver a patial differetial equation to ordinary differential equation with or without dchange?
 

Download pde_to_ode.mw

i want to plot this equations for P and t for  t= 0 t0 600 seconds..can u provide me how?

time-dependent_aerodynamic_pressure.mw

Download time-dependent_aerodynamic_pressure.mw

Download convert-latex.mw

Hi everybody

I have a differential equation where derivatives (for example d/dr) have been displayed by prime notation.

I want to convert this equation latex format by keeping the prime notation not d/dr.

I would appreciate it if anyone could help me.

I have attached a sample code.

Thank you

Hadi

Hello everyone! I have found this document from Maple 18 (2014) that will be very usefull for my school work, but when i copy it into Maple 2022 i get the error message, "action, does not evaluate to a module" Is there any way for me to fix this? What i have been able to read is that i have to go in and update the code, but i have absolute no idea on how to do it? 

Kind Regards Samuel 

Regressionmodeller._Skal_laves_til_maple_22_og_23.mw

could any one help  to plot this.

Here is my code.

PWS.mw

I need to make the graph i have attached match the image below.

Is there any wrong with unknown parameter i found?How to solve this.Please Help.

Download bug.mw

Hi

I think there is a bug in the "randpoly" command. please see the attached file line 7 of my procedure "randomzero". Why x^2-y-4 is created while terms=1 is considered and the outputs must contain binomial?

How to get U1,U2,..I dont know how to use this inverse transform.Please help to find the series.

CM.mw

Equilibrium._solutions.mw

I need to obtain the equilibrium solutions when infection is present as well as analyze the local asymptotic stability when infection is absent and present in the system

Hi, please can someone help on how non-dimensionalize PDEs. 

I have tried the following, but is not working:

restart:
eqn := (diff(theta(x, z, t), x))^2*(K[1]-K[3])*cos(theta(x, z, t))*sin(theta(x, z, t))+(diff(theta(x, z, t), x))*((diff(theta(x, z, t), z))*(-K[1]*cos(2*theta(x, z, t))+K[3]*cos(2*theta(x, z, t)))-(1/2)*gamma[1]*(4*sin(theta(x, z, t))^2*u(x, z, t)+2*u(x, z, t)*cos(2*theta(x, z, t))))+(diff(theta(x, z, t), z))^2*(K[3]-K[1])*cos(theta(x, z, t))*sin(theta(x, z, t))-(1/2)*gamma[1]*(diff(theta(x, z, t), z))*(4*sin(theta(x, z, t))^2*v(x, z, t)+2*v(x, z, t)*cos(2*theta(x, z, t)))+(diff(theta(x, z, t), z, x))*(-2*K[1]+2*K[3])*cos(theta(x, z, t))*sin(theta(x, z, t))-(diff(u(x, z, t), z))*((1/2)*gamma[2]*cos(2*theta(x, z, t))+(1/2)*gamma[1]*(2*sin(theta(x, z, t))^2+cos(2*theta(x, z, t))))-(diff(v(x, z, t), x))*((1/2)*gamma[2]*cos(2*theta(x, z, t))+(1/2)*gamma[1]*(-2*sin(theta(x, z, t))^2-cos(2*theta(x, z, t))))-(1/2)*gamma[1]*(4*sin(theta(x, z, t))^2*(diff(theta(x, z, t), t))+2*(diff(theta(x, z, t), t))*cos(2*theta(x, z, t)))+((diff(u(x, z, t), x))*gamma[2]-(diff(v(x, z, t), z))*gamma[2])*cos(theta(x, z, t))*sin(theta(x, z, t))+f[2](theta(x, z, t))*(diff(theta(x, z, t), x, x))+f[1](theta(x, z, t))*(diff(theta(x, z, t), z, z));

varchange := {t = T*tau, u = xi*h^2*U/alpha[4], v = xi*h^2*V/alpha[4], x = X*h, z = Z*h, K[3] = K[1]*k[3], f[1] = K[1]*F[1], f[2] = K[1]*F[2], gamma[1] = mu*Gamma[1], gamma[2] = mu*Gamma[2]};

PDEtools:-dchange(varchange, eqn, [tau, U, V, X, Z, k[3], F[1], F[2], GAMMA[1], GAMMA[2]]);

i want to label x,y axis with command...if i use manually ..i always need to do again after running..so how can i label x,y axis in this plot...? i watch manual and i din't find for axis..stress-strain.mwstress-strain.mw

in this equations...how to plot (a,e)? thank you guys

plot.mw

Can you change f(eta) to upflow curve and theta(eta) to downflow curve.

In my Problem,Boundary Conditions are

theta(infinity) = 0, (D(f))(infinity) = 1 , (Take, eta =infinity)

Flows will be correct for what value is taken for infinity .

I take  eta = 5. and also tried changing ranges  but could't find it.Please Help to fix the curve.

my code is,

SM.mw

How to solve ordinary differemtial equation system with initial conditions and boundary conditions. Here, some initial conditions are unknown variables. So how to find these  values of parameters.

eq1 := diff(f(x), x, x, x)+(1/2)*cos(alpha)*x*(diff(f(x), x, x))+(1/2)*sin(alpha)*f(x)*(diff(f(x), x, x)) = 0;

eq2 := diff(g(x), x, x)+diff(g(x), x)+(diff(g(x), x))*(diff(h(x), x))+cos(alpha)*x*(diff(g(x), x))+sin(alpha)*f(x)*g(x) = 0;

eq3 := diff(g(x), x, x)+diff(h(x), x, x)+1/2*(cos(alpha)*x+sin(alpha)*f(x)) = 0

ics:=f(0)=0, f'(0)=1, f''(0)=a[1], g(0)=1, g'(0)=a[2], h(0)=1, h'(0)=a[3];

bcs:=f(x) , g(x), h(x) tends to 0 ad x tends to infinity

How to find series values.I got this error.Please Help.

Maple code for the problem is

TFBE.mw