Advanced Search

Saha

70 Reputation

4 Badges

3 years, 75 days
Contact Saha

MaplePrimes Activity


These are questions asked by Saha

How to solve this solution in R-K Method...

Saha 70 Maple
ode homework dsolve
June 19 2022
1 4


 

Please Help! To Clear this Error and Plot this series in R-K Method
> 

restart
> 

with(PDETools, declare):
> 

with(LinearAlgebra):
> 

with(PolynomialIdeals):
> 

with(plots):
> 

declare(f(x));

f(x)*`will now be displayed as`*f

 (1)
> 

eq1 := (diff(f(x), `$`(x, 3)))/((1-phi)^2.5*(1-phi+phi*rho[s]/rho[f]))+(1/2)(`cosα`)*eta*(diff(f(x), `$`(x, 2)))+(1/2)(`sinα`)*f(x)*(diff(f(x), `$`(x, 2)))-M(diff(f(x), `$`(x, 1))-1)/(1-phi+phi*rho[s]/rho[f]) = 0;

(diff(diff(diff(f(x), x), x), x))/((1-phi)^2.5*(1-phi+phi*rho[s]/rho[f]))+(1/2)*eta*(diff(diff(f(x), x), x))+(1/2)*f(x)*(diff(diff(f(x), x), x))-M(diff(f(x), x)-1)/(1-phi+phi*rho[s]/rho[f]) = 0

 (2)
> 

eq2 := K[nf]*(diff(Theta(x), `$`(x, 2)))/K[f]+(1/2)*Pr*(1-phi+phi*rho[s]*c[s]/(rho[f]*c[f]))(eta*`cosα`+f*`sinα`)*(diff(Theta(x), `$`(x, 1))) = 0;

K[nf]*(diff(diff(Theta(x), x), x))/K[f]+(1/2)*Pr*(1-phi(eta*`cosα`+f*`sinα`)+phi(eta*`cosα`+f*`sinα`)*rho[s](eta*`cosα`+f*`sinα`)*c[s](eta*`cosα`+f*`sinα`)/(rho[f](eta*`cosα`+f*`sinα`)*c[f](eta*`cosα`+f*`sinα`)))*(diff(Theta(x), x)) = 0

 (3)
> 

bcs := eta = 1, rho[s] = 5200, rho[f] = 997.1, c[s] = 670, c[f] = 4179, K[s] = 6, K[f] = .613, K[nf] = K[f]*(K[s]+2*K[f]-2*phi*(K[f]-K[s]))/(K[s]+2*K[f]+phi*(K[f]-K[s])):

f(0) = 0, (D(f))(0) = 0, Theta(0) = 1, Theta(10) = 0, (D(f))(10) = 1

 (4)
> 

a1 := [cos(beta) = 1, sin(beta) = 0, M = 1, Pr = 6.2, phi = .1];

[cos(beta) = 1, sin(beta) = 0, M = 1, Pr = 6.2, phi = .1]

  

[cos(beta) = 0, sin(beta) = 1, M = 1, Pr = 6.2, phi = .1]

  

[cos(beta) = 1/2, sin(beta) = (1/2)*3^(1/2), M = 1, Pr = 6.2, phi = .1]

  

[cos(beta) = 1, sin(beta) = 0, M = 1, Pr = 6.2, phi = .1]

  

[cos(beta) = 0, sin(beta) = 1, M = 1, Pr = 6.2, phi = 0.]

  

[cos(beta) = 1/2, sin(beta) = (1/2)*3^(1/2), M = 1, Pr = 6.2, phi = 0.]

 (5)
> 

b1 := subs(a1, eq1);

1.301348831*(diff(diff(diff(f(x), x), x), x))/(.9+.1*rho[s]/rho[f])+(1/2)*eta*(diff(diff(f(x), x), x))+(1/2)*f(x)*(diff(diff(f(x), x), x))-(1)(diff(f(x), x)-1)/(.9+.1*rho[s]/rho[f]) = 0

  

1.301348831*(diff(diff(diff(f(x), x), x), x))/(.9+.1*rho[s]/rho[f])+(1/2)*eta*(diff(diff(f(x), x), x))+(1/2)*f(x)*(diff(diff(f(x), x), x))-(1)(diff(f(x), x)-1)/(.9+.1*rho[s]/rho[f]) = 0

  

1.301348831*(diff(diff(diff(f(x), x), x), x))/(.9+.1*rho[s]/rho[f])+(1/2)*eta*(diff(diff(f(x), x), x))+(1/2)*f(x)*(diff(diff(f(x), x), x))-(1)(diff(f(x), x)-1)/(.9+.1*rho[s]/rho[f]) = 0

  

1.301348831*(diff(diff(diff(f(x), x), x), x))/(.9+.1*rho[s]/rho[f])+(1/2)*eta*(diff(diff(f(x), x), x))+(1/2)*f(x)*(diff(diff(f(x), x), x))-(1)(diff(f(x), x)-1)/(.9+.1*rho[s]/rho[f]) = 0

  

1.000000000*(diff(diff(diff(f(x), x), x), x))+(1/2)*eta*(diff(diff(f(x), x), x))+(1/2)*f(x)*(diff(diff(f(x), x), x))-1.*(1)(diff(f(x), x)-1) = 0

  

1.000000000*(diff(diff(diff(f(x), x), x), x))+(1/2)*eta*(diff(diff(f(x), x), x))+(1/2)*f(x)*(diff(diff(f(x), x), x))-1.*(1)(diff(f(x), x)-1) = 0

 (6)
> 

c1 := subs(a1, eq2);

K[nf]*(diff(diff(Theta(x), x), x))/K[f]+3.100000000*(1-(.1)(eta*`cosα`+f*`sinα`)+(.1)(eta*`cosα`+f*`sinα`)*rho[s](eta*`cosα`+f*`sinα`)*c[s](eta*`cosα`+f*`sinα`)/(rho[f](eta*`cosα`+f*`sinα`)*c[f](eta*`cosα`+f*`sinα`)))*(diff(Theta(x), x)) = 0

  

K[nf]*(diff(diff(Theta(x), x), x))/K[f]+3.100000000*(1-(.1)(eta*`cosα`+f*`sinα`)+(.1)(eta*`cosα`+f*`sinα`)*rho[s](eta*`cosα`+f*`sinα`)*c[s](eta*`cosα`+f*`sinα`)/(rho[f](eta*`cosα`+f*`sinα`)*c[f](eta*`cosα`+f*`sinα`)))*(diff(Theta(x), x)) = 0

  

K[nf]*(diff(diff(Theta(x), x), x))/K[f]+3.100000000*(1-(.1)(eta*`cosα`+f*`sinα`)+(.1)(eta*`cosα`+f*`sinα`)*rho[s](eta*`cosα`+f*`sinα`)*c[s](eta*`cosα`+f*`sinα`)/(rho[f](eta*`cosα`+f*`sinα`)*c[f](eta*`cosα`+f*`sinα`)))*(diff(Theta(x), x)) = 0

  

eq4

  

K[nf]*(diff(diff(Theta(x), x), x))/K[f]+3.100000000*(1-(0.)(eta*`cosα`+f*`sinα`)+(0.)(eta*`cosα`+f*`sinα`)*rho[s](eta*`cosα`+f*`sinα`)*c[s](eta*`cosα`+f*`sinα`)/(rho[f](eta*`cosα`+f*`sinα`)*c[f](eta*`cosα`+f*`sinα`)))*(diff(Theta(x), x)) = 0

  

K[nf]*(diff(diff(Theta(x), x), x))/K[f]+3.100000000*(1-(0.)(eta*`cosα`+f*`sinα`)+(0.)(eta*`cosα`+f*`sinα`)*rho[s](eta*`cosα`+f*`sinα`)*c[s](eta*`cosα`+f*`sinα`)/(rho[f](eta*`cosα`+f*`sinα`)*c[f](eta*`cosα`+f*`sinα`)))*(diff(Theta(x), x)) = 0

 (7)
> 

d1 := dsolve({b1, bcs}, {bcs, c1}, numeric); d1(0)

Error, (in dsolve/numeric/process_input) invalid argument: {eta = 1, K[f] = .613, K[nf] = K[f]*(K[s]+2*K[f]-2*phi*(K[f]-K[s]))/(K[s]+2*K[f]+phi*(K[f]-K[s])), K[s] = 6, c[f] = 4179, c[s] = 670, rho[f] = 997.1, rho[s] = 5200, K[nf]*(diff(diff(Theta(x), x), x))/K[f]+3.100000000*(.9+.1*rho[s](eta*`cosα`+f*`sinα`)*c[s](eta*`cosα`+f*`sinα`)/(rho[f](eta*`cosα`+f*`sinα`)*c[f](eta*`cosα`+f*`sinα`)))*(diff(Theta(x), x)) = 0}

  

d1(0)

 (8)
> 

``

Warning,Solving for expressions other th...

Saha 70 Maple
solve duplicate-question
May 22 2022
2 0


 

Please,Help the Maple code for find the values of F[3],F[4],F[5],........How to rectify this error?
> 

restart
> 

F[0] := 0;

0

  

0

  

A

 (1)
> 

for k from 0 to 10 do F[k+3] := solve(F[k+3]+((1/2)*(sum(cos*beta*delta(k-m-1)*(m+1)*(m+2)*F[m+2], m = 0 .. k))-(1/2)*sin(beta)*(sum(F[k-m]*(m+1)*(m+2)*F[m+2], m = 0 .. k))+M(k+1)*F[k+1])/((k+1)(k+2))(k+3), F[k+3]) end do

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

  

-cos*beta*(table( [( -1 ) = 0, ( 0 ) = 1, ( 1 ) = 0, ( 2 ) = 0, ( 3 ) = 0, ( 4 ) = 0, ( 5 ) = 0, ( 6 ) = 0, ( 7 ) = 0, ( 9 ) = 0, ( 8 ) = 0, ( 10 ) = 0 ] ))('k-1')*A

  

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

  

"F[4]:="

  

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

  

Error, (in thaw) invalid input: thaw uses a 1st argument, e, which is missing

  
> 

``

"How to find the values of this series? ...

Saha 70 Maple
solve recurrence
May 02 2022
3 5

How to find the values of X(1),X(2),..&Y(1),Y(2)...Plese help .
> 

restart;
 

for k from 0 to 5 do
X(k+1):=solve(2*(k+1)*X(k+1)+(k+1)*Y(k+1)-X(k)-Y(k)+(1)/k!,X(k+1));
Y(k+1):=solve((k+1)*X(k+1)+(k+1)*Y(k+1)+2*X(k)+Y(k)+(1)/k!,Y(k+1)); od;
> 


 

2

  

1

  

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

  

Error, (in solve) a constant is invalid as a variable, -(1/2)*Y(1)+1

  
> 

``

Download DE_Using_DTM-Ex-3(1).mw

 

"Help!for this series."...

Saha 70 Maple
solve homework
April 23 2022
0 1


could you please help me ,the maple code for this given series.

> 

restart
> 

U[0](x) := x;

x

 (1)
> 

"U[k+1](x):=solve((k+1)*U[k+1](x)+(x*(-1)^((k-1)/(2)))/(k!)-x^(2)*((e)^(x))/(10){6/(k!)-sum((2^(k[1]))/(k[1!])(5*delta[k[]-k[1]](x)+(2^(k-k[1])*(-1)^((k-k[1])/(2)))/((k-k[1])!)+(2^(k-k[1]+1)*(-1)^((k-k[1]-1)/(2)))/((k-k[1])!)),k[1]=0..k)}-(cos^(2)(x)+sin^()(x))*((∂)^2)/((∂)^( )x^2) [U[k](x)]-(e)^(x)[sum(1/(k[1]!){1/(k-k[1])(sum(sum(1/(k[3]!)*U[k[2]-k[3]](x)*U[k-k[1]-k[2]-1](x)},k[1]=0..k),k[2]=0..k-k[1]-1),k[3]=0..k[2])),U[k+1](x)]; od;"

Error, unable to match delimiters

"U[k+1](x):=solve((k+1)*U[k+1](x)+(x*(-1)^((k-1)/2))/(k!)-x^2*((e)^x)/10{6/(k!)-sum((2^(k[1]))/(k[1!])(5*delta[k-k[1]](x)+(2^(k-k[1])*(-1)^((k-k[1])/2))/((k-k[1])!)+(2^(k-k[1]+1)*(-1)^((k-k[1]-1)/2))/((k-k[1])!)),k[1]=0..k)}-(cos^2(x)+sin(x))*((∂)^2)/(∂x^2) [U[k](x)]-(e)^x[sum(1/(k[1]!){1/(k-k[1])(sum(sum(1/(k[3]!)*U[k[2]-k[3]](x)*U[k-k[1]-k[2]-1](x)},k[1]=0..k),k[2]=0..k-k[1]-1),k[3]=0..k[2])),U[k+1](x)]; od;"

  
> 

``


 

Download Chapter_6-Example-6.5.4.mw

''How do I solve a differential equation...

Saha 70 Maple
April 21 2022
1 2

I  can't solve this series,please could you help me ,the maple code for this.

restart;
U[0](x):=(e)^(x);
                             exp(x)
n:=7:  for k from 0 to n do  U[k+1](x):=solve((k+1)U[k+1](x)-((∂)^2)/((∂)^( )x^2) U[k](x)-2 U[k](x));  end do;
Error, (in solve) cannot solve expressions with diff(diff(U[0](x), x), x) for x

First 6 7 8 9 Page 8 of 9