Maple Questions and Posts

These are Posts and Questions associated with the product, Maple

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

``

``

restart:with(plots):

mb:=765 : mp:=587 :Ib:=76.3*10^3 :Ip:=7.3*10^3 :l:=0.92 :d:=10:F:=0.042:omega:=0.56 :

eq1:=(mb+mp)*diff(x(t),t$2)+mp*l*cos(alpha(t))*diff(alpha(t),t$2)+mp*[l*[diff(alpha(t),t)]^2*sin(alpha(t))]-F*sin(omega*t)=0;

1352*(diff(diff(x(t), t), t))+540.04*cos(alpha(t))*(diff(diff(alpha(t), t), t))+[540.04*[diff(alpha(t), t)]^2*sin(alpha(t))]-0.42e-1*sin(.56*t) = 0

(1)

eq2:=-mp*l*sin(alpha(t))*diff(alpha(t),t$2)+mp*[l*[diff(alpha(t),t)]^2*cos(alpha(t))]-9.81*(mp+mb)-F*sin(omega*t)=0;

-540.04*sin(alpha(t))*(diff(diff(alpha(t), t), t))+[540.04*[diff(alpha(t), t)]^2*cos(alpha(t))]-13263.12-0.42e-1*sin(.56*t) = 0

(2)

eq3:=mp*[d+l*cos(alpha(t))]*diff(x(t),t$2)+[Ip+mp*l^2+mp*d*l*cos(alpha(t))]*diff(alpha(t),t$2)-mp*sin(alpha(t))*[-l*d*alpha(t)^2]+mp*[l*9.81*sin(alpha(t))]=0;

587*[10+.92*cos(alpha(t))]*(diff(diff(x(t), t), t))+[7796.8368+5400.40*cos(alpha(t))]*(diff(diff(alpha(t), t), t))-587*sin(alpha(t))*[-9.20*alpha(t)^2]+[5297.7924*sin(alpha(t))] = 0

(3)

eq4:=mp*l*cos(alpha(t))*diff(x(t),t$2)+(Ip+mp*l^2)*diff(alpha(t),t$2)-mp*9.81*l*sin(alpha(t))=0;

540.04*cos(alpha(t))*(diff(diff(x(t), t), t))+7796.8368*(diff(diff(alpha(t), t), t))-5297.7924*sin(alpha(t)) = 0

(4)

CI:= x(0)=0,alpha(0)=0,D(x)(0)=0,D(alpha)(0)=0;

x(0) = 0, alpha(0) = 0, (D(x))(0) = 0, (D(alpha))(0) = 0

(5)

sys := eq1, eq2, eq3, eq4:

``

solution:=dsolve([sys,CI],numeric);

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

 

 

NULL


thanks for your help

Download tangage.mw

Is it possible to numerically calculate  the integral

int((-12*y^2+1)*ln(abs(Zeta(x+I*y)))/(4*y^2+1)^3, [y = 0 .. infinity, x = 1/2 .. infinity])

in Maple?

The code

int((-12*y^2+1)*ln(abs(Zeta(x+I*y)))/(4*y^2+1)^3, [y = 0 .. infinity, x = 1/2 .. infinity],numeric,epsilon=0.1)

has been executed on my comp  without any output since this morning.

 

 

 

Hello everybody, i need to graphic a couple of functions just like this one:

 http://temasmatematicos.uniandes.edu.co/Casquetes_cilindricos/Pags/Anim_1.htm

i have been watching this:

 http://www.maplesoft.com/teachingconcepts/detail.aspx?cid=12 VISUALIZATION --> Animation 2

i've tried with (plots) (plottools) animate, etc. but i can't figure out how to do it. 

It would be very helpful if someone explain me how to do this.

Thank you all!

I have setup two groups in MapleCloud, to be used for distribution of materials (a readonly group) and as a drop for students to deposit their exercise sheets (an opaque group). It is my understanding that e.g. students have to request joining a group and I get to approve them.

How does one join a group managed by someone else? If I try from another account I cannot see the groups I just created, so I cannot request joining them. I seem also to not be able to just approve someone (by their email e.g.) without that person having requested to join first.

???

Thanks

USPAS2014

Is there a way of ploting y=x^x for when x<0 and y is a noncomplex solution?

 

I have some questions about Mobius Project:

  1. What is the license of the worksheets posted on Mobius Project's website?
    I can not find any information about the license of the materials submitted to the Mobius Project.
    Is my worksheets are protected by intellectual property laws and copyrighted by me?
    Or it is freely available for any use for everyone?
  2. What are the terms of use of Mobius Project and its service? I can not find it anywhere.
  3. What is the copyright for http://mobius.maplesoft.com/?
    I do not see any copyright in the footer of the site or elsewhere.
    So, no copyrights mean the materials/resources on this website can be freely copied/used by anyone?

with pointplot3d and 14,000 points when I enter symbol=point I get an empty plot.

Only when I set symbolsize=1 (a point) do I get points appearing in the graph.  Bug?

Dear Users

I have a problem for solving a system of linear equations that arise from collocation method for getting approximate solution of a coupled PDE and ODE in Food engineering problems.

When it reach to the fsolve command it takes long time!!!

I used maple 13.

If kindly is possible, please help me in this special case.

With kind regards,

Emran Tohidi.

 

> Restart;
print(`output redirected...`); # input placeholder
> h := 50; hm := 0.1e-3; rhodp := 1500; Y := 0.5e-1; T0 := 20; rhoair := 1.2041; Dair := 0.2e-8; DD := 0.85e-9; C := 3240; L := 0.4e-1; X0 := 1.5; V := .2; delta := 0.2e-2; Yair := 0.5e-1; nu := .2; Tair := 60; Hnu := 2400; rho := 1359; tt := 3;
%;
> N := 5; Digits := 20;
> X := sum(sum(a[m, n]*z^m*t^n, m = 0 .. N), n = 0 .. N); X := unapply(X, z, t); Xt := diff(X(z, t), `$`(t, 1)); Xt := unapply(Xt, z, t); Xz := diff(X(z, t), `$`(z, 1)); Xz := unapply(Xz, z, t); Xzz := diff(X(z, t), `$`(z, 2)); Xzz := unapply(Xzz, z, t); T := sum(b[n]*t^n, n = 0 .. (N+1)^2-1); T := unapply(T, t); Tt := diff(T(t), `$`(t, 1)); Tt := unapply(Tt, t); aw := exp(.914)*X(z, t)^.5639-.5*exp(1.828)*X(z, t)^(2*.5639); aw := unapply(aw, z, t); TT := 8.3036+3816.44*(1+T(t)/(46.13)+T(t)^2/46.13^2)/(46.13); TT := unapply(TT, t); pwv := 133.3*(1+TT(t)+(1/2)*TT(t)^2); pwv := unapply(pwv, t); Yi := .622*pwv(t)*aw(z, t)*(1+pwv(t)*aw(z, t)/rho+(pwv(t)*aw(z, t)/rho)^2)/rho; Yi := unapply(Yi, z, t);
%;
> S1 := {seq(seq(Xt(delta*i/N, tt*j/N)-DD*Xzz(delta*i/N, tt*j/N) = 0, i = 1 .. N-1), j = 1 .. N)};
> S2 := {seq(DD*rhodp*Xz(delta, tt*j/N)+hm*rhoair*Yi(delta, tt*j/N) = 0, j = 0 .. N)};
> S3 := {seq(Xz(0, tt*j/N) = 0, j = 0 .. N)};
> S4 := {seq(X(delta*i/N, 0) = 0, i = 1 .. N-1)};
> S5 := {seq(seq(rho*delta*C*Tt(tt*j/N)-h*(Tair-T(tt*j/N))+hm*Hnu*rhoair*(Yair-Yi(delta*i/N, tt*j/N)) = 0, j = 1 .. N), i = 0 .. N)};
print(`output redirected...`); # input placeholder
> S6 := {seq(rho*delta*C*Tt(0)-h*(Tair-T0)+hm*Hnu*rhoair*(Yair-Yi(delta*i/N, 0)) = 0, i = 0 .. N)};
%;
> SS := `union`(`union`(`union`(`union`(`union`(S1, S2), S3), S4), S5), S6);
> sol := fsolve(SS);

Hello, I have two plots:

Is it possible to display differene between these plots?

I don't understand how to create better parameters. Why is the graph not flush with the image (it appears that the graph is a set of points and NOT lines)

 

From the manual: coeffs - extract all coefficients of a multivariate polynomial. Is there a way of doing the reverse (giving the coefficients, obtain the multivariate polynomial) ? For univariate polynomials I know that the answer is yes because  PolynomialTools[FromCoefficientList] - return a univariate polynomial from list of coefficients. But what about multivariate polynomials?

ACP.mw

 

hi all

trying to modify some kitonum code to get the smallest solution, d=3,515,820, but i have problem....

http://en.wikipedia.org/wiki/Archimedes'_cattle_problem

 

hello,

i have unistalled maple, now i have to give the license server name and the purchase code to reactivate it, how can i do?? where i can find them?

someone can help me?!

Hello,

 

  I have a question. Consider

 

fsolve(x^2+3*x+1=3, x);

 

  I want to save the two roots into two variables. What kind of commend shall I use? 

 

P.S. My further aim comes from solving an equation without analytical solution. Therefore I cannot plug in the solution formula. 

 

 

restart;
Eq1 := diff(T1(t), t) = (W*Cp*(To-T1(t))+UA*(Ts-T1(t)))/(M*Cp);
Eq2 := diff(T2(t), t) = (W*Cp*(T1(t)-T2(t))+UA*(Ts-T2(t)))/(M*Cp);
Eq3 := diff(T3(t), t) = (W*Cp*(T2(t)-T3(t))+UA*(Ts-T3(t)))/(M*Cp);
sys := Eq1, Eq2, Eq3;

Operational Veriables

W := 100;
UA := 10;
Cp := 2;
M := 1000;
To := 20;
Ts := 250;

Initial Conditions

sys1 := {Eq1, Eq2, Eq3};

nsys := nops(sys1);

ics := {T1(0) = 20, T2(0) = 20, T3(0) = 20};
{T1(0) = 20, T2(0) = 20, T3(0) = 20}
nics := nops(ics);
for i from 1 to nics do Sol ||i:=dsolve({sys1, ics[i]},{T1(t),T2(t),T3(t)},numeric)od;
Error, unable to match delimiters
Typesetting:-mambiguous(Typesetting:-mambiguous( for i from 1 to

nics do Sol verbarverbariAssigndsolvelpar(sys1comma ics(i))

commalcubT1(t)commaT2(t)commaT3(t)rcubcommanumericrparod,

Typesetting:-merror("unable to match delimiters")))

 

First 1368 1369 1370 1371 1372 1373 1374 Last Page 1370 of 2248