javid basha jv

## 130 Reputation

6 years, 303 days

## Kindly do the needful for n=1.5, 2.5 ca...

Dear @Carl Love

I have tried for the n=1,2,3 like a round value; this is working. But when the values like 1.5 and 1.8, it is not executing.
Kindly do the needful to find a solution

Thank you.

## Is there any possible to compute the dec...

Dear maple users

Kindly do the needful to find a solution for the decimal power.

Thank you.

## For the decimal numerator, it is not wor...

In this equation I have the numerator in decimal("^n-1")

It works when power("^n") is 1,2...and so on.
But when power("^n") is decimal, like 0.8, 1.2, 1.4, it is not working.

Is there any way to tackle such a situation

(1+(diff(f(x, t), x))^(n-1))*(diff(f(x, t), x))/(x*R(z)^2)+(diff(f(x, t), x))^(n-1)*(n-1)*(diff(f(x, t), x, x))/R(z)^2

## Is this possible to rectify "HFloat(unde...

Dear @acer

I reworked the earlier question and found a new issue, so I made a new question.

I am getting the following issue. I though It makes an Imaginary part.

unable to store HFloat(undefined)+HFloat(undefined)*I when datatype=float[8]
Is this possible to get a solution?

1.mw

 > restart:
 > with(LinearAlgebra):
 > with(PDEtools):
 >
 > with(plots):
 > fcns := {f(x,t)};
 (1)
 > #GOLD
 > b1:=1.41:d:=0.5:xi:=0.1:ea:=0.5:c1:=1:eo:=0.2:c2:=1;tau:=0.5;m:=1.5:n:=0.8:ra:=1;alpha:=Pi/4;b2:=2:kappa:=Pi/4:
 (2)
 >
 >
 >
 > R := proc (z) options operator, arrow; piecewise(d <= z and z <= d+1/2, 1-2*xi*(z-d), d+1/2 <= z and z <= d+1, 1+(-1)*.5*xi*(1+cos(2*Pi*(z-d-.5))), 1) end proc;
 (3)
 > PDE1 :=-ra*tau*(diff(f(x,t),t))+b1*(1+ea*cos(c1*t))+b2*(cos(c2+kappa))+((m+(diff(f(x, t), x))^(n-1))/(x*R(z)*R(z))*diff(f(x, t), x)+(diff(f(x, t), x)/(R(z)*R(z)))*diff(m+(diff(f(x, t), x))^(n-1),x)+(1/(R(z)*R(z)))*(m+(diff(f(x, t), x))^(n-1))*diff(f(x, t), x\$2));
 (4)
 >
 > IBC := {D[1](f)(0,t)=0,f(1,t)=0,f(x,0)=0}:
 > z:=0.98:
 >
 > sol1:=pdsolve(evalf({PDE1}),IBC ,numeric, time = t,spacestep = 0.025, timestep=0.001): sol:=((evalf(Re(sol1))), output=listprocedure);
 (5)
 > p1 := (sol):-plot(f(x, t), t = 1.5,numpoints = 100);
 >
 >

## Is there any possible...

Hi maple users

Is there any possibility of executing this code at sa=4?.
Kindly do the needful.

## I found a solution...

I have found the results for sa=1 case, But it is not executing at sa=2 or 3 cases.

Is there any way to get the solutions

Kindly do the needful.

 > restart:   with(plots):   PDEtools[declare](f1(x),f2(x),f3(x), t1(x),t2(x),t3(x)):
 (1)
 > gr:=5: pa:=5: sa:=1: br:=0.1: A1:=1: A2:=2: A3:=1:   Eq11:= (1-p)*((diff(F1, x\$2)+gr*T1)+pa)+p*((diff(F1, x\$2)+gr*T1)+pa):   Eq12:= (1-p)*(diff(T1, x\$2))+p*((diff(T1, x\$2)+br*(diff(F1, x))*(diff(F1, x)))):   Eq21:= (1-p)*((diff(F2, x\$2)+gr*A1*A2*T2)-sa*sa*F2+pa*A1)+p*((diff(F2, x\$2)+gr*A1*A2*T2)-sa*sa*F2+pa*A1):   Eq22:= (1-p)*(diff(T2, x\$2))+p*((diff(T2, x\$2)+A1*A3*br*((diff(F1, x))*(diff(F1, x))+sa*sa*F2*F2))):   Eq31:= (1-p)*br*((diff(F3, x\$2)+gr*T3)+pa)+p*((diff(F3, x\$2)+gr*T3)+pa):   Eq32:= (1-p)*br*(diff(T3, x\$2))+p*((diff(T3, x\$2)+br*(diff(F3, x))*(diff(F3, x)))):
 > for i from 0 to N+1 do       equ1[i] := coeff(Eq11, p, i) = 0:       equ2[i] := coeff(Eq12, p, i) = 0:       equ3[i] := coeff(Eq21, p, i) = 0:       equ4[i] := coeff(Eq22, p, i) = 0:       equ5[i] := coeff(Eq31, p, i) = 0:       equ6[i] := coeff(Eq32, p, i) = 0:   end do:
 > con1[0]:= f1[0](-1) = 0, f1[0](0) = f2[0](0), D(f1[0])(0) = D(f2[0])(0),             f2[0](1) = f3[0](1), D(f2[0])(1) = D(f3[0])(1),f3[0](2)=0:   con2[0]:= t1[0](-1) = 0, t1[0](0) = t2[0](0), D(t1[0])(0) = D(t2[0])(0),             t2[0](1) = t3[0](1), D(t2[0])(1) = D(t3[0])(1), t3[0](2)=1:   for h to N do       con1[h]:= f1[h](-1) = 0,  f1[h](0) = f2[h](0), D(f1[h])(0) = D(f2[h])(0),                 f2[h](1) = f3[h](1), D(f2[h])(1) = D(f3[h])(1), f3[h](2)=0:       con2[h]:= t1[h](-1) = 0, t1[h](0) = t2[h](0), D(t1[h])(0) = D(t2[h])(0),                 t2[h](1) = t3[h](1), D(t2[h])(1) = D(t3[h])(1), t3[h](2)=0:   end do:
 > for i from 0 to N do       P:= dsolve( {con1[i], con2[i], equ1[i], equ2[i], equ3[i], equ4[i], equ5[i], equ6[i]},                   {f1[i](x), t1[i](x), f2[i](x), t2[i](x), f3[i](x), t3[i](x)}                 ):       f1[i](x):=rhs(P[1]):       f2[i](x):=rhs(P[2]):       f3[i](x):=rhs(P[3]):       t1[i](x):=rhs(P[4]):       t2[i](x):=rhs(P[5]):       t3[i](x):=rhs(P[6]):   end do;
 (2)
 >

## Solution...

Hi maple users,

Is there any way to get the solutions?
Kindly do the needful.

## Is there any alter way to tackle the err...

I hope everything goes well

Is there any alter way to tackle the error

I have noticed one of your replaies, I have tried that one also.   I found the error.

Kindly look at one of your replaies

https://www.mapleprimes.com/questions/224981-Invalid-Subscript-Selector

Kindly do the needful to find a solution.
Thank you.

## we need to calculate the integral consta...

Dear @vv

Here we need to calculate the stream function. So should be calculated the integral constant values that time we can found the stream function.

Is this possible to found an integral constant(C1) value.

Dear @vv

Have a good day.

## Thank you...

Dear @tomleslie

Thanks for your explanation. Now I got clarity. I will cross-check the equation.

Have a good day.

## Thank you for your response...

Dear @tomleslie

I have to calculate the values like

L(0.71)=0 than f(x,t)=?

L(0.71)=0.1 than f(x,t)=?

..

L(0.71)=1 than f(x,t)=?.

and the x and t are 0.71 and 1.12.

But here various values of Z only the answer is coming.

The actual values are,

L(0.71)=0 than f(x,t)=0.5859

L(0.71)=0.1 than f(x,t)=0.5829

...

L(0.71)=1 than f(x,t)=0.

How to get the actual values.

## Thank you...

Dear@Carl Love

Thanks for the help

Have a good day.

## Thank you...

Dear @Kitonum

Thanks for the help

Have a good day.

## Thank you so much...

Dear @tomleslie

Thank you so much.