javid basha jv

## 130 Reputation

6 years, 336 days

## How to change a computed data value ...

Maple 18

Dear maple users,

A fine day wishes to all.

How to replace a value in already computed data value.

In my problem, I have evaluated a set of values (with help of #ma1 := evalf(seq((A4(t)), t = 0.0..1, 0.1));)

ma1 := 0., .1941703021, .3203871063, .4089371834, .4712881303, .5145114133, .5435036431, .5617715009, .5718586242, .5756277760, .5744585726

I need to change 0 value as a 1.

I have used subs(ma[1]=1,ma1), but its coming wrong.

## How to find out f'(x,t) and f''(x,t) val...

Maple 18

Dear maple users,

In this code, how to find out the f'(x,t) and f''(x,t) values.
How to export the computed values in the excel file.JVB.mw

 > restart:
 > with(PDEtools):
 > with(plots):
 > fcns := {f(x,t)};
 (1)
 > ra:=2:b1:=1.41:na:=0.7:we:=0.5:eta[1]:=4*0.1:d:=0.5/1:xi:=0.1:m:=na:ea:=0.5:pr:=21: gr:=0.1:
 > R:=0.9323556933;
 (2)
 > PDE1 :=ra*(diff(f(x,t),t))=+b1*(1+ea*cos(t))+(1/(R^2))*((diff(f(x,t),x,x))+(1/x)*diff(f(x,t),x));
 (3)
 > IBC := {D[1](f)(0,t)=0,f(1,t)=0,f(x,0)=0};
 (4)
 > sol :=  pdsolve({PDE1}, IBC, numeric,spacestep = 0.025, timestep=0.0001) ;
 (5)
 > sol:-plot[display](f(x, t), t = 1.2, linestyle = "solid", title = "Velocity Profile", labels = ["r", "f"]);
 >

## How to plot the figure at time ...

Dear maple users,

Greetings.

In this code, the problem has executed. But I unable to get the graph.
Kindly do the needful to plot the figure at time t=0.45.
Also, please explain how to import the computed values into an excel file.

 > restart:
 > with(PDEtools):
 > with(plots):
 > ra:=1:b1:=4:na:=0.7:we:=0.5:eta[1]:=4*0.1:beta:=0.5/1:xi:=tan(0):m:=na:ea:=0.5:
 > R:=z->piecewise(beta<=z,(1+xi*z)*(1-(64/10)*eta[1]*((11/32)*(z-beta)-(47/48)*(z-beta)^2+(z-beta)^3-(1/3)*(z-beta)^4)),z<=beta+(3/2),(1+xi*z)*(1-(64/10)*eta[1]*((11/32)*(z-beta)-(47/48)*(z-beta)^2+(z-beta)^3-(1/3)*(z-beta)^4)),(1+xi*z));
 (1)
 >
 > PDE :=ra*(diff(f(x,t),t))=+b1*(1+ea*cos(t))+(1/(x*R(z)^2))*diff((x*(m+(1-m)*(1+(we*we*((diff(f(x,t),x)))^((na-1)/2)))^(-1))*(diff(f(x,t),x))),x);
 (2)
 > IBC := {D[1](f)(0,t)=0,f(1,t)=0,f(x,0)=0};
 (3)
 > z:=0.77:
 >
 > sol := pdsolve(PDE, IBC, numeric, spacestep = 0.1e-1);
 (4)
 >

 >

## How to solve this PDE equation...

Maple 18

Dear maple users

Greetings.

I hope you are all fine.

In this code, I am solving the PDEs via pdsolve with numeric.

There is some mistake in the boundary condition and pdsolve.

Kindly help me that to get the solution for this PDE.

In this problem h(z) is piecewise

Bc:

code:JVB.mw

Note: z=0.5:

## How to plot the excel data...

Maple 18

Dear maple users,

Greeting for the day.

I have some data in an excel file. How to plot the excel data in contours with labels.

Also, how to show the x-axis range =-1.5..1.5 and y-axis range = -0.5..1 with the contour sequence -2..2.

XLdat.xlsx