MaplePrimes Questions

Let

and f=

The elements of W are none zero. I want a procedure that return "true" if f is none zero w.r.t. W and return

"false" otherwise.

hello 

 

if you have a function lets say: 2x+1/4-x3

 

 

Now if you have to plot the graf, how should you know what the x and y shoud be? I mean you do like: plot(f(x),x=..

 

what should it be? i find it hard :(

Hello, 

 

I'm trying built a graph in cylindrical coordinates. I have 3 matrix X(20,30); Y(20,30), and temperatura T(20,30). The matrices X and Y are already in cylindrical coordinates. In matlab is just "pcolor (X,Y,T)". Anyone can help me?

X and Y are the coordinates of cylinder discretized (z was not included)

I have this plot in matlab, but i need it in maple.

Thanks

Hello, forgive me if I used bad english, I am not a native speaker.

 

Anyhow: I have to decide the set of coordinates looking at the grafs intersections.  (I hope you understand that)

 

This is my function: f(x)=1/4*x3-x2-x+4

 

When you plot the function you see that the intersections is (-2,0) (2,0) and (4,0)

 

BUT I have to use a command to find these three intersections, plzz help me!

 

THANK YOU indeed.

 

 

[point of] intersection

Are there any commands in maple that will help me find a suitable function that approximates the numerical solution of:



  restart;
  PDE := diff(v(x, t), t) = diff(v(x, t), x, x);
  JACOBIINTEGRAL := int(JacobiTheta3(0, exp(-Pi^2*s))*v(1, t-s)^4, s = 0 .. t);
  IBC:= D[1](v)(0,t)=0,
        D[1](v)(1,t)=-0.000065*v(1, t)^4,
        v(x,0)=1;
#
# For x=0..1, t=0..1, the solution varies only very slowly
# so I have increased the timestep/spacestep, just to speed
# up results generation for diagnostic purposes
#
  pds := pdsolve( PDE, [IBC], numeric, time = t, range = 0 .. 1,
                  spacestep = 0.1e-1, timestep = 0.1e-1,
                  errorest=true
                )

diff(v(x, t), t) = diff(diff(v(x, t), x), x)

 

int(JacobiTheta3(0, exp(-Pi^2*s))*v(1, t-s)^4, s = 0 .. t)

 

(D[1](v))(0, t) = 0, (D[1](v))(1, t) = -0.65e-4*v(1, t)^4, v(x, 0) = 1

 

_m649569600

(1)

#
# Plot the solution over the ranges x=0..1,
# time=0..1. Not a lot happens!
#
  pds:-plot(x=1, t=0..1);

 

#
# Plot the estimated error over the ranges x=0..1,
# time=0..1
#
  pds:-plot( err(v(x,t)), x=1,t=0..1);

 

#
# Get some numerical solution values
#
  pVal:=pds:-value(v(x,t), output=procedurelist):
  for k from 0 by 0.1 to 1 do
      pVal(1, k)[2], pVal(1, k)[3];
  od;

 

t = 0., v(x, t) = Float(undefined)

 

t = .1, v(x, t) = .999977377613528229

 

t = .2, v(x, t) = .999967577518313666

 

t = .3, v(x, t) = .999959874331053822

 

t = .4, v(x, t) = .999952927885405241

 

t = .5, v(x, t) = .999946262964885979

 

t = .6, v(x, t) = .999939702966688881

 

t = .7, v(x, t) = .999933182128311282

 

t = .8, v(x, t) = .999926675964661227

 

t = .9, v(x, t) = .999920175361791563

 

t = 1.0, v(x, t) = .999913676928735229

(2)

 

 

 

Download PDEprob2_(2).mw

 

I am refering to the first graph, is there a way in maple to find an explicit suitable approximating function?

I.e, I want the function to have the same first graph obviously, it seems like addition of exponent and a line function, I tried plotting exp(-t)-0.3*t, it doesn't look like it approximates it very well. Any suggestion on how to implement this task in maple?

Thanks.

 

I'm have used a program to find the roots of a function 

 

f:=x*cos(x)-sin(x)*sin(x/1000);
/ 1 \
x cos(x) - sin(x) sin|---- x|
\1000 /

x_max:=50; x_min:=-50; step:=2; i_max:=(x_max-x_min)/step;
50
-50
2
50

j:=1:
for i from 0 to i_max by 1 do
x0:=x_min+i*step;
x[j]:=fsolve(f=0,x=x0);
j:=j+1;
end:

 

and my output was of the form of multiple "potential" roots and a bunch of which are the same. So I tried to get rid of the ones which were the same before actually finding the ones which ARE roots. To do that I done....

 

 

j := 1; for j to 50 do if x[j]-x[j+1] = 0 then ignore(x[j]) else print(x[j]) end if end do:

 

and it got rid of the ones which are of the above form but some roots are the same and seperated by more than 1 ... i.e x[ j ]= x[j + 2] or some other number. 

 

Basically I am trying to generalise the above for loop for all "numbers" instead of 1 but when I try some things the for loop doesnt like it. 

Any help would be good!

 

 

Say I have an equation of ax^2 + by^2 = 0,

I would like to plot a graph of y^2 against x^2...

How do i do that?

 

Hello i want to solve the differential equation but i have these problem i don't understand  why !?

restart;
with(LinearAlgebra):
with(student):
P:=180000:
Vt:=P/m;
m:=1.4787880*10^5-t*606.06:
g:=9.81:
T:=176:
pi:=evalf(Pi):

Euler_x := -(diff(lambda3(t), t)) = 0;
Euler_y := -(diff(lambda4(t), t)) = 0;
Euler_Vx := -lambda3(t)-(diff(lambda1(t), t)) = 0;
Euler_Vy := -lambda4(t)-(diff(lambda2(t), t)) = 0;
Euler_theta := lambda1(t)*sin(theta)-lambda2(t)*cos(theta) = 0;
transversalnost:=y(t)-lambda1(t)*Vt*cos(theta)-lambda2(t)*Vt*sin(theta)+lambda2(t)*g-lambda3(t)*Vx(t)-lambda4(t)*Vy(t)=0;

transversalnost:=eval(transversalnost,t=176);

sys:={Euler_x, Euler_Vx, Euler_y, Euler_Vy, theta=arctan(lambda2(t)/lambda1(t)), diff(Vx(t),t)=Vt*cos(theta), diff(Vy(t),t)=Vt*sin(theta)-g, diff(y(t),t)=Vy(t), diff(x(t),t)=Vx(t),x(0)=0, y(0)=0.302, lambda1(0)=0, lambda2(0)=0, lambda3(0)=0, lambda4(0)=0,transversalnost, x(176)=120};

ans:=dsolve(sys, numeric, output=listprocedure):

Error, (in dsolve/numeric/bvp/convertsys) unable to convert to an explicit first-order system

thanks for your help

 

Download 1234.mw

Hi everyone,

 

I'm trying to solve the following eqauation but Maple gives me the answer (( RootOf(mexp(-_Z*(m-1))*d^2-theta+_Z*theta-theta*c*t__kj) ))

 

The equation is:

solve(mexp(-(m-1)*t__ij)*d^2-theta+theta*t__ij-theta*(sum(t__kj, k = 1 .. c))-m*eta*(diff((1-1/exp(t))^m, t)) = 0, t__ij);

 

Could you please help me??

 

What is the meaning rootOF? Is there any explicit solution to that equation??

 

Thank you for your help

Please how do I solve the above 2nd order differential equation (F and K are constants)?

I have the following boundary conditions:

Many thanks!

Hi everyone,

I'm trying to solve the following equation and unfortunately I get this error ( Error, (in Engine:-Dispatch) invalid subscript selector ) each time I try. Could you please help me??

I will be thankful for your great help.

solve(m*exp(-(m-1)*tij)*d2-θ+θ*tij-θ*(Σ(tkj, k = 1 .. c))-m*η*(1-exp(-tij))(m-1) = 0, tij)

Hello !!

I want to replace the word "assuming" with another word. I tried to write:

> as:=assuming;

> alias(as=assuming);

> `&as`:=proc(a,b) a assuming b; end;

The first two methods are not allowed. The last one doesn't work like "assuming".

Is there any solution to this problem ?

 

How does one 3D plot the simple complex exponential e^I2pift or cos(2pift)+I*sin(2pift)  where f is freqency and t is time.  It shoud display a spiral aroung the time axis.

This may be a trivial question, but does this factor fully with the newer versions of Maple, say at 900 digits?

 

Digits:=900;

rho_poly := -2201506283520*rho^32+(-17612050268160+104204630753280*I)*rho^31+(2237195146493952+737798139150336*I)*rho^30+(14065203494780928-29153528496783360*I)*rho^29+(-260893325886750720-161432056834818048*I)*rho^28+(-1240991775275876352+1727517243589263360*I)*rho^27+(8952004373272068096+6696323263091441664*I)*rho^26+(25553042370906292224-37948239682297921536*I)*rho^25+(-135024511500569280512-65293199430849134592*I)*rho^24+(-79740262928225402880+401487130320847241216*I)*rho^23+(956745211126674882560-164797793704574713856*I)*rho^22+(-1213375867282228772864-1655554058430246551552*I)*rho^21+(-1483956336776821211136+3604946201834409820160*I)*rho^20+(6525094787202650144768-1597915397190007586816*I)*rho^19+(-8575469412912592879616-6168391294117580865536*I)*rho^18+(2408139380338842796032+15004449784317106323456*I)*rho^17+(10583091471310114717696-17047513330720373194752*I)*rho^16+(-22619716982813548707840+8898637295768494915584*I)*rho^15+(26538067620972845277184+5129530051326543351808*I)*rho^14+(-21415800164460070789120-17268159356969925234688*I)*rho^13+(11916012071577094946816+22601135173030541677568*I)*rho^12+(-3551246770922037813248-21229478915196610975744*I)*rho^11+(-977434486760953073664+16249214903618313346048*I)*rho^10+(1977414870691507931136-10721551032564274826240*I)*rho^9+(-1197394212949208115968+6172794574205050632192*I)*rho^8+(280273257275327368320-2996290081120136529792*I)*rho^7+(108849195761508531648+1152454823926345101504*I)*rho^6+(-119736267114490955904-327757949185254534784*I)*rho^5+(49149411853848597568+63563541902968683712*I)*rho^4+(-11524495997215059744-7307364351434838944*I)*rho^3+(1585189353379709888+299568910286253408*I)*rho^2+(-116032795768295808+25487628220230528*I)*rho+3299863116538269-2454681763039104*I;;

If I do even expansion for the function 2*x-2, x from [0,1], to the real line,how do I make Maple to produce

the Fourier cosine expansion in the form 2*x-2=a_0/2+Sum(a_n cos n*Pi*x)

 

?

 

First 1210 1211 1212 1213 1214 1215 1216 Last Page 1212 of 2429