Advanced Search

delvin

35 Reputation

5 Badges

2 years, 347 days
Contact delvin

MaplePrimes Activity


These are replies submitted by delvin

@mmcdara I want to solve equation 3...

delvin 35
November 23 2023
0 0

@mmcdara 

I want to solve equation 3-4. For this, we assume that we have the equation 2-9 where the values of m = n = 3. Then the answers to 3-5 and 3-6 are obtained. But it does not happen in practice???

Why didn't you consider the first two terms?

@dharr   Thanks, Now the las...

delvin 35
November 18 2023
0 0

@dharr 

 Thanks,

Now the last answers of the program (photo) are not available. What should I do?

@mmcdara   I really appreciat...

delvin 35
October 09 2023
0 0

@mmcdara 

 I really appreciate your taking the time to explain things.

@mmcdara   Does it mean that ...

delvin 35
October 09 2023
0 0

@mmcdara 

 Does it mean that sometimes luck works and it is wrong? There is no way to get values like the previous file that answered?

@acer  I appreciate your taking th...

delvin 35
September 29 2023
0 0

@acer 

I appreciate your taking the time.

@acer I appreciate it....

delvin 35
September 29 2023
0 0

@acer 

I appreciate it.

@mmcdara Much obliged.Yes, I don...

delvin 35
September 29 2023
0 0

@mmcdara 

Much obliged.

Yes, I don't know why this happened.

I can't do the second part either, that is, when the number becomes a fraction, the variable is changed.

@acer  I’m so grateful. How...

delvin 35
September 28 2023
0 0

@acer 

I’m so grateful.

How do I do this?

03.mw

@mmcdara Thank you for helping me.B...

delvin 35
September 28 2023
0 0

@mmcdara 

Thank you for helping me.

But unfortunately there is an error.

``

``

restart

alias(U = U(xi))

``

"Balance:=proc(EQ)"

Error, unterminated procedure

"Balance:=proc(EQ)"

  

_local(eq); 1; eq := eval(EQ, Diff = diff)

EQ

 (1)

_local(HODD); 1; HODD := degree(eval(eval(eq, U = exp(__k*xi)), xi = 0), __k)

0

 (2)

_local(HNLD); 1; HNLD := degree(eval(eq, U = __C), __C)

0

 (3)

M = solve(HNLD*M = HODD+M)

M = 0

 (4)

"end proc:"

Error, unable to parse

"end proc:"

  

eq := a*U+b*U^2+diff(U, `$`(xi, 2))+diff(U, xi)

``

Balance(eq)

Balance(a*U+b*U^2+diff(diff(U, xi), xi)+diff(U, xi))

 (5)

eq := a*U+b*U^2+Diff(U, `$`(xi, 2))+Diff(U, xi)

``

Balance(eq)

Balance(a*U+b*U^2+Diff(U, xi, xi)+Diff(U, xi))

 (6)

eq := a*U+b/U^2+diff(U, `$`(xi, 5))+diff(U, `$`(xi, 3))+c*U^3; Balance(eq)

Balance(a*U+b/U^2+diff(diff(diff(diff(diff(U, xi), xi), xi), xi), xi)+diff(diff(diff(U, xi), xi), xi)+c*U^3)

 (7)

``

Download Balance.mw

How can we solve EQ equations?...

delvin 35
September 26 2023
0 0

Hello
How can we solve EQ equations?

``

``

``

NULL

``

> 

restart
> 

with(student)

U := a[0]+sum(-a[i]*tanh(xi[n])^i, i = 1 .. 1)+sum(-b[i]*tanh(xi[n])^(-i), i = 1 .. 1)

a[0]-a[1]*tanh(xi[n])-b[1]/tanh(xi[n])

 (1)
> 

u(xi[n+1]) := a[0]-a[1]*(tanh(xi[n])+tanh(d))/(1+tanh(xi[n])*tanh(d))-b[1]*(1+tanh(xi[n])*tanh(d))/(tanh(xi[n])+tanh(d))

a[0]-a[1]*(tanh(xi[n])+tanh(d))/(1+tanh(xi[n])*tanh(d))-b[1]*(1+tanh(xi[n])*tanh(d))/(tanh(xi[n])+tanh(d))

 (2)
> 

u(xi[n-1]) := a[0]-a[1]*(tanh(xi[n])-tanh(d))/(1-tanh(xi[n])*tanh(d))-b[1]*(1-tanh(xi[n])*tanh(d))/(tanh(xi[n])-tanh(d))

a[0]-a[1]*(tanh(xi[n])-tanh(d))/(1-tanh(xi[n])*tanh(d))-b[1]*(1-tanh(xi[n])*tanh(d))/(tanh(xi[n])-tanh(d))

 (3)

eq := diff(U, `$`(xi[n], 2))-(diff(U, xi[n]))-1-u(xi[n-1])+2*U-u(xi[n+1])

2*a[1]*tanh(xi[n])*(1-tanh(xi[n])^2)-2*b[1]*(1-tanh(xi[n])^2)^2/tanh(xi[n])^3-2*b[1]*(1-tanh(xi[n])^2)/tanh(xi[n])+a[1]*(1-tanh(xi[n])^2)-b[1]*(1-tanh(xi[n])^2)/tanh(xi[n])^2-1+a[1]*(tanh(xi[n])-tanh(d))/(1-tanh(xi[n])*tanh(d))+b[1]*(1-tanh(xi[n])*tanh(d))/(tanh(xi[n])-tanh(d))-2*a[1]*tanh(xi[n])-2*b[1]/tanh(xi[n])+a[1]*(tanh(xi[n])+tanh(d))/(1+tanh(xi[n])*tanh(d))+b[1]*(1+tanh(xi[n])*tanh(d))/(tanh(xi[n])+tanh(d))

 (4)
> 

U := a[0]-a[1]*tanh(xi[n])-b[1]/tanh(xi[n])

a[0]-a[1]*tanh(xi[n])-b[1]/tanh(xi[n])

 (5)

simplify(eq); fin1 := simplify(numer(%))

-2*tanh(xi[n])^10*tanh(d)^2*a[1]-tanh(xi[n])^9*tanh(d)^2*a[1]+(2*a[1]*tanh(d)^4+2*a[1])*tanh(xi[n])^8+(a[1]*tanh(d)^4+(a[1]+b[1]-1)*tanh(d)^2+a[1])*tanh(xi[n])^7+(-2*b[1]*tanh(d)^4+2*tanh(d)^2*b[1]-2*a[1])*tanh(xi[n])^6+((-a[1]-b[1]+1)*tanh(d)^4+(-a[1]-b[1])*tanh(d)^2-a[1]-b[1]+1)*tanh(xi[n])^5+(-2*a[1]*tanh(d)^4+2*a[1]*tanh(d)^2-2*b[1])*tanh(xi[n])^4+(b[1]*tanh(d)^4+(a[1]+b[1]-1)*tanh(d)^2+b[1])*tanh(xi[n])^3+(2*b[1]*tanh(d)^4+2*b[1])*tanh(xi[n])^2-tanh(xi[n])*tanh(d)^2*b[1]-2*tanh(d)^2*b[1]

 (6)

subs(tanh(xi[n]) = Psi, fin1); fin := simplify(%)

(-2*(Psi^6*a[1]+(1/2)*Psi^5*a[1]+(-(1/2)*a[1]-(1/2)*b[1]+1/2)*Psi^3+(1/2)*Psi*b[1]+b[1])*(Psi+1)^2*(Psi-1)^2*cosh(d)^4+2*(Psi+1)*(Psi-1)*(Psi^8*a[1]+(1/2)*Psi^7*a[1]-Psi^6*a[1]+(-a[1]-(1/2)*b[1]+1/2)*Psi^5+(-a[1]+b[1])*Psi^4+((1/2)*a[1]+b[1]-1/2)*Psi^3+Psi^2*b[1]-(1/2)*Psi*b[1]-b[1])*cosh(d)^2+2*Psi^2*(Psi^6*a[1]+(1/2)*Psi^5*a[1]-Psi^4*b[1]+(-(1/2)*a[1]-(1/2)*b[1]+1/2)*Psi^3-Psi^2*a[1]+(1/2)*Psi*b[1]+b[1]))/cosh(d)^4

 (7)

for i from 0 to degree(fin, Psi) do EQ[i] := simplify(coeff(fin, Psi, i)) end do

-2*sinh(d)^2*a[1]/cosh(d)^2

 (8)

example("??")

``

Download 01.mw

@Carl Love  It is solved for anoth...

delvin 35
August 18 2023
0 0

@Carl Love 

It is solved for another equation.

How is u defined in Maple?...

delvin 35
July 31 2023
0 0

How is u defined in Maple?

@mmcdaraI REALLY APPRECIATE IT. Yes...

delvin 35
July 30 2023
0 0

@mmcdara

I REALLY APPRECIATE IT. 

Yes, unfortunately, I have many faults.
What is add? Why didn't you put sum?

I want a code that solves a differential equation. in such a way that it shows the relationship of 10 in the same way despite the presence of psi. And in the next series, the psi function should be defined and then it should be in such a way that the values of r and s are given for example 1, -2, i... be) and then put in the main equation to simplify the equation and get the unknown values.

``

``

> 

restart
> 

with(student)
> 

NULL

u(xi[n]) := a[0]+sum(a[i]*psi(xi[n])^i, i = 1 .. 1)+sum(b[i]*psi(xi[n])^(-i), i = 1 .. 1)+sum(c[i]*(diff(psi, xi[n]))^i/psi(xi[n])^i, i = 1 .. 1)

a[0]+a[1]*psi(xi[n])+b[1]/psi(xi[n])

 (1)
> 

NULL
> 

NULL
> 

u(xi[n+1]) := u(xi[n]+d)

u(xi[n]+d)

 (2)
> 

u(xi[n-1]) := u(xi[n]-d)

u(xi[n]-d)

 (3)
> 

NULL

eq := c*(diff(U, xi[n]))*(U(xi[n])+u(xi[n-1]))*(U(xi[n])+u(xi[n+1]))-(2*(u(xi[n-1])-u(xi[n+1])))*(U(xi[n])^2)(1-U(xi[n])^2)

-2*(u(xi[n]-d)-u(xi[n]+d))*(U(xi[n]))(1-U(xi[n])^2)^2

 (4)

``

Download abs123.mw

Thanks...

delvin 35
July 28 2023
0 0

@C_R 

How do I do this? Can you modify my instructions?

@C_R  Difficulty in drawing a 3D d...

delvin 35
July 28 2023
0 0

@C_R 

Difficulty in drawing a 3D diagram
Hello, I want to draw a 3D diagram that has i inside, but it gives me an error. Can anyone help me?

clc
clear
close all
% plot grid
t = linspace(-5, 5, 100);
n = linspace(-5, 5, 100);
% compute plot grid values
[T, N] = meshgrid(t, n);

% compute Z values
Z = -2.*sin(2)*(cos(T).*sin(T).*i - cos(T).^2 + 1./2).*cos(2)./((4.*cos(2).^2 - 2).*sin(T).*cos(T));

% plot
surf(T, N, Z);
title(' Diagram3D With Surfc','fontsi',20,'interpreter','latex')
xlabel('X','fontsi',20,'interpreter','latex')
ylabel('Y','fontsi',20,'interpreter','latex')
zlabel('Z','fontsi',20,'interpreter','latex')
colorbar;
colormap hsv

 

1 2 3 4 5 6 7 Page 2 of 7