Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

Hi, 

 

I wish to be able to calculate the roots of the function f(p) by using the roots of the function h(p) and applying the bisection method due to the fact that the roots of h(p) bracket the roots of f(p) as can be seen in the graph below. I have done this before for another example when h:=J0(p); and hence i could use The commands BesselJZeros(0,n)/BesselJZzeros(0,n+1) to find the roots. So my problem arises with the finding the roots of h(p) and how to insert them into my bisection loop(underlined below).

Any advice would amazing. Many thanks


restart;

Digits := 30:
with (plots):
with(RootFinding):

#Define given parameters

R:=1: #external radius of particles, cm

d:=10^(-3): #diffusivity cm^2 per second

alpha:= 1: #fractional void volume

c0:=10^(-6): #concentartion of soltion in void volume of solid initially, moles per liter

C0:=0: #concentration of main body of solution initially, moles per liter

k1:=0.5: #constant in adsorption isotherm (ka)

k2:=0.75: #constant in adsorption isotherm (kd)

k:=2.5: #equilbrium constant for adsorption kinetics

n0:=(k1/k2)*c0:#initial amount absrobed on solid, moles per liter

V:=0.1: #volume of external solution, liters

W:=0.1: #weight of absorbant, grams

rho:=2.0: #solid aparrant density, g/cc

delta:=(1/d)*((p+alpha*k2+k1)/(p+alpha*k2));

beta:=W*alpha*d/(rho*V);

   

1000*(p+1.25)/(p+.75)

 

0.500000000000000000000000000000e-3

 

 

 

 

 


f:=p->(BesselJ(0,R*sqrt(-delta*p))*k*p-(R*sqrt(-delta*p))*BesselJ(1,R*sqrt(-delta*p))*(d*p/R + 2*beta*k/(R^2)));

proc (p) options operator, arrow; BesselJ(0, R*sqrt(-delta*p))*k*p-R*sqrt(-delta*p)*BesselJ(1, R*sqrt(-delta*p))*(d*p/R+2*beta*k/R^2) end proc

(2)

h:=p->(BesselJ(0,R*sqrt(-delta*p)));
plot([f(p),h(p)],p=-0.3..0,axis=-5..5,legend=["f(p)","h(p)"]);


 




proc (p) options operator, arrow; BesselJ(0, R*sqrt(-delta*p)) end proc

 

 

 

(3)



points:=5:
rts:= Array(1..points):
for n from 1 by +1 to points do
pl:=evalf(#**first root of h**);
pu:=evalf(#**second root of h, i.e n+1 root**);
pe:= (pl+pu)/2;
while abs(f(pe))>10^(-6) do
pe:=(pu+pl)/2;
if f(pu)*f(pe) <0 then
pl:=pe;
elif f(pl)*f(pe)<0 then
pu:=pe;
end if;
od;
rts[n]:=pe;

od;
rts[n]:=p[n];

 

 

Download spherical_continue.mw


 

Download spherical_continue.mw

Hi,

 

I am trying to plot union of regions in Maple 15. Each region is formed by a set of inequalities. I can see from online such a feature is available in Maple 17.

The code is:

> with(plots);

> I1 := inequal({x < 1, y < 1}, x = 0 .. 2, y = 0 .. 2, optionsfeasible = [color = "Brown", thickness = 2]);

> I2 := inequal({y > 1, x < 1}, x = 0 .. 2, y = 0 .. 2, optionsfeasible = [color = "Yellow", thickness = 2]);

> display(I1);      ---> This is working 

>display(I2);      ---> This is also working

>display({I1,I2});  ---> Not working.

The above is also attached in the worksheet.

Is there any other way to plot the two region. I am new to Maple. Please suggest me the alternative in Maple 15.

 

Thanks!

 

Hello there,

i got a question regarding derivatives in Matlab.

 

I got a function for a example:

f:=f1(y)*f2(x)

 

from this function i need the partial derivatives with respect to x and y, i can easily get them with

diff(f,x)

diff(f,y)

 

now i want to compute the derivative of f with respect to time t, assuming that y and x both depend on t - how can i tell maple that y and x depend on t?

 

thanks

ben

I want to find the solution in a special form.
How can I do it?
Here is what I tried:

(Maple)

(Maple)


In the left hand side u_1 is not changed in  D(u_1).
I want to substitute and evalute (differentiate) it.

Thanks,  Sandor

 

 

 

 

assume input 0.567

how to get a list [5, 6, 7] or [5;6;7]?

(a) Design your own 3-stage explicit Runge-Kutta method with one-step error O(h4).

(b) Test your method by solving y= −y. Confirm that the global error in your numerical solution

is O(h3).

Write a Maple procedure that solves for y(1) in the initial value problem

                     y= f(y),     y(0) = 1,

 

 

using a numerical stencil based on the nth order Taylor series expansion of y. The procedure’s arguments should include an arbitrary function f, an integer n representing the accuracy of the Taylor series expansion, another integer N representing the number of steps between x = 0 and x = 1. Pick a test problem and compare your results with the output of dsolve/numeric.

Hello Maplers, i have encountered a little annoyance with Maple, that i would like to ask, whether it can be solved. 

 

It's when i try to define a function, like f(x)=2x and try to define it with a command f(x):=2x, a pop-up box comes up, asking me whether i'd like to use a 'function definition', or 'remember table assignment', and i would like to make Maple remember my choice that i want a function. 

 

I know i can write it like f:=x->2x, but i hate to look at that, to be frankly..

 

So, is there any way to solve this?

 

I am using "Maple 15" and "Windows 7" and I want to share the tip for using "Classic worksheet" and "GR Tensor". Windows 7 is 64-bit but you can use Classic worksheet and GR Tensor, which work basically in 32-bit, by simply downlording "32-bit Maple 15" !! I haven't tried other more current versions of Maple but Maple 15 works well: you can choose 32-bit Maple in the download options.

Enjoy the Classic worksheet and GR Tensor. Good luck !  

Please tell me how to do about the following problem to me.

 

g:=(b*y)^k*k*y;

simplify(%);

Then, what I obtained was (b*y)^k*k*y, not (b^k)*k*y^(k+1).

expand of the command brought the same answer not (b^k)*k*y^(k+1).

Please tell me what was wrong to my calculation.

 

Thank you in advance.

 

Taro.

 

How to solve this equation

cos(x)^2+2*x*cos(x)-2*sin(x)+4=0

with

solve(%,x) and evalf(%) ?

I get no solution. Is there any option to force Maple to solve this equation by solve and evalf?


Student[Calculus1]:-NewtonsMetnod(%,x=2) solves the problem.
If I use Series for sin(x) and cos(x) Maple can solve it too.

Thank you

 

Dear all,

A stupid question.

I tried to create a density plot of the data but I only saw the method to plot graphs through equations. (Actually, I'm not very familiar with related statistical commands in Maple)

This is the histograph of my data:

And I want to get the density plot like this:

It seems that we can get the graph through "InteractiveDataAnalysis" but cannot print the graph into the worksheet or assign the graph to the variable.

How to create a density plot through data? Thanks.

AOA... I want to export system of linear equation into matrix form...

Help.mw

restart:

Eq1:=S*diff(f(x,t),x,t)+diff(f(x,t),x)^2-f(x,t)*diff(f(x,t),x$2)=diff(f(x,t),x$3);

BCs := {D[1](f)(0,t)=cos(t), f(0,t)=0,D[1](f)(L,t)=0};

ICs := {f(x,0)=0};

S:=10:L:=5:
smod3:= pdsolve(Eq1,ICs union BCs,numeric,range=0..L);

smod3:-plot(t=0,  color=red):

it seems to me that the problem is due to the mixed bcs. Any way around?

Cheers!

Question_.pdf

 

this question about midpoint method I need help with the part c

 

thanks

First 1366 1367 1368 1369 1370 1371 1372 Last Page 1368 of 2223