Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

I've got the following four differential equations :

v_x:=diff(x(t),t);
v_y:=diff(y(t),t);
d2v_x:=-((C_d)*rho*Pi*(r^2)*(v_x)*sqrt((v_x)^2 +(v_y)^2))/(2*m);
d2v_y:=-((C_d)*rho*Pi*(r^2)*(v_y)*sqrt((v_x)^2 +(v_y)^2))/(2*m)-g;


and the following initial value conditions:

x(0)=0,y(0)=0,v_x(0)=v0/sqrt(2),v_y(0)=v0/sqrt(2) given v0=65 

I need to solve these using the numeric type and then draw overlaid plots

(i) setting C_d=0

(ii) leaving C_d as a variable

before plotting y(t) vs x(t). The hint for this last part is that the path can be seeing using [x(t),y(t)] instead of [t,y(t)]

I've tried to do it but seemed to have several syntax errors.

 

 

I've got the following diff.eq

y'(x)=sin(x*y(x)) given y(0)=1 

and need to solve it numerically which is why I've used:

dy4:=diff(y(x),x);
eqn4:=dy4=sin(x*y(x));
ic1:=y(0)=1;
ans3:=dsolve({eqn4,ic1},y(x),type=numeric);

This code doesn't return a value though and in fact, ans3 is being displayed as a procedure

"ans3:=proc(x_rkf45) ... end proc"

I don't quite understand why and what I need to do to get the required numerical solution

 

Hello

Could you plese help me to plot my function. i want to plot a function in which variables are have able formatting. I have attache my file.

Thank you.


restart

L[b] := 400:

L[c] := 400:

q := 5:

M[ab] := VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(2, E), I[c]), 1/L[c]), theta[b]):``

M[ba] := VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(2, E), I[c]), 1/L[c]), VectorCalculus:-`*`(2, theta[b])):

M[bc] := VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(2, E), I[b]), 1/L[b]), VectorCalculus:-`+`(VectorCalculus:-`*`(2, theta[b]), theta[c])), VectorCalculus:-`*`(VectorCalculus:-`*`(q, L[b]^2), 1/12)):

M[cb] := VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(2, E), I[b]), 1/L[b]), VectorCalculus:-`+`(theta[b], VectorCalculus:-`*`(2, theta[c]))), VectorCalculus:-`-`(VectorCalculus:-`*`(VectorCalculus:-`*`(q, L[b]^2), 1/12))):

M[cd] := VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(2, E), I[c]), 1/L[c]), VectorCalculus:-`*`(2, theta[c])):

M[dc] := VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(VectorCalculus:-`*`(2, E), I[c]), 1/L[c]), theta[c]):

eq1 := VectorCalculus:-`+`(VectorCalculus:-`*`(VectorCalculus:-`+`(M[ab], M[ba]), 1/L[c]), VectorCalculus:-`*`(VectorCalculus:-`+`(M[cd], M[dc]), 1/L[c])) = 0:

eq2 := VectorCalculus:-`+`(M[ba], M[bc]) = 0:

solve([eq1, eq2], {theta[b], theta[c]})

{theta[b] = -(40000000/3)/(E*(I[b]+2*I[c])), theta[c] = (40000000/3)/(E*(I[b]+2*I[c]))}

(1)

theta[b] := rhs({theta[b] = -(40000000/3)/(E*(I[b]+2*I[c])), theta[c] = (40000000/3)/(E*(I[b]+2*I[c]))}[1]):

theta[c] := rhs({theta[b] = -(40000000/3)/(E*(I[b]+2*I[c])), theta[c] = (40000000/3)/(E*(I[b]+2*I[c]))}[2]):

subs(theta[b] = theta[b], M[ba])

-(400000/3)*I[c]/(I[b]+2*I[c])

(2)

subs(theta[b] = theta[b], M[bc]):

``

EXplore(plot(M[ba], I[b] = 1 .. 10), parameters = [I[c] = 1 .. 20])

Error, (in plot) unexpected option: I[b] = 1 .. 10

 

``


Download shiboft2.mw

Matlab seems to be pretty strong at doing color plots with separate color bars, e.g.

Is this also possible in Maple and somehow in combination with `plots[surfdata](...,color=zhue,...)`?

In a 3d coordinate system I have a circular spacecurve with z-minimum -4 and z-maximum +4. In the same 3d coordinate system I have a 3d surface plot with z-minimum -0.5 and z-maximum +1.3 . When I choose the color option "Z(Hue)" in order to color-code the z-values on the 3d surface and make the topography more clear, I mostly get a totally green 3d surface. It seems that the color scaling is coupled with the spacecurve with z-values of +-4 . How can I uncouple the color scaling from the spacecurve and couple it with the z-range of the surface, while the color-limits shall be at +-1.3 ?

I need to find the local maxima and minima of f(x,y)=x(x+y)*e^(y-x). I have tried to look for an appropriate method that I could use to achieve this, but got stuck. I also don't quite understand the math behind tying to obtain the local maxima and  minima for a function of this type.

I'm trying to evaluate the multidimensional limit: 

(1+y)^(x-1)-1/(1-cos((x-1)^2+y^2)^(1/4)) as (x,y)->(1,0) using the limit command :

limit((1+y)^(x-1)/(1-cos((x-1)^2+y^2)^(1/4)),{x=1,y=0});

 

but don't seem to get any output. Also, I think the limit for this function doesn't exist or is indeterminate on R2. Where am I wrong?

I am considering options to access Maple features from within a Fortran program using the OpenMaple API. I do not find any examples illustrating this (or even precise statements if this is possible in principle). (Quite surprising to me considering the presence of Fortran in scientific computing).

Is the only way to go to write a C wrapper for the API? If yes, is there an example for such a wrapper which performs common conversions between Fortran and C data types.

Thanks!

Greetings, seeking an expert to animate a plot.

see worksheet.posterior_graphs_(encapsulted)_1D.mw

before they play each other, each have a law (a normal distribution) plot-output 6.

after DD defeats CC, and a numerical integration is performed the new laws are given by plot-output 18.

as you can see, the laws of DD and CC are closer together.

if the calc was repeated (DD defeats CC again), the laws would be closer again.

so what i require is an animation of the new laws from game 1 to (say) game 6 (DD defeats CC every time). seeing the red and blue distributions merging would be ideal.

as an aside I heard maples FFT could simplify the complicated integration. any suggestions?

cheers

I have been Maple  18 with no problem. Then, today, the  "=" symbol and  "+" symbol don't work.  Does anyone know what to do about this?

 

Hi,
There a lot of symbols that don't work in Maple.
The symbols that dont work are shown as an "A".
Almost half of the symbols in the pallets are shown as an "A".

Hi,

I am trying to discretize a kernel of the form $K(x,y,t,s)$. I want to evaluate a four dimensional integral of the form

\int\int\int\int K(x,y,t,s)*h_m(x)*h_n(y)*h_p(t)*h_q(s) dsdtdxdy, where limits of integration are from 0 to 1.

$h_m()$ are function of one variable.

please suggest how to evaluate this.

thanks

 

In the attached Maple worksheet I attempt to plot the solution of an initial value problem for a first order ODE.  DEplot fails with a cryptic message.  Strangely enough, if I give the "arrows=none" option to DEplot, it produces the correct plot!

I see this behavior in Maple 17 and 18.

Maple 11, however, works fine with or without the "arrows=none" option.

Is there an explanation for this or is it a bug?

DEplot-bug.mw

First 1293 1294 1295 1296 1297 1298 1299 Last Page 1295 of 2248