MaplePrimes Questions

How do I plot the graph of 3 functions into one page (that is, as a series). Thank you.

Please this is urgent! I can't seem to be able to upload my job on here.

Thanks.

Hi Everyone, i have question about Piechart.. http://www.mmsonline.com/cdn/cms/Sandvik-pie-chart.jpg

is something like this possible in Maple?

 

Thank you

Hello guys,

 

I'm starting to learn MAPLE this week, so first of all I'm sorry for my lack of knowledge about the software. 

I got a worksheet with a ready code, and i'm trying to run it, but I'm receiving an error message that I don't know how to solve, by the fact that I don't know very well yet. But I need to run this URGENTLY, so my hope is that you guys could help me.

 

The line of code is the following: 

 

S1:-Results(solutionpoint);for cont2 from 1 to 2 do VNUM[jj,cont2]:=(rhs(S1:-Results(solutionpoint)[de+s+m+cont2]))end do;seq(-add(v[k]*x[jj,k],k=1..de)+add(u[r]*y[jj,r],r=1..s)<=0,j=1..n);seq(add(u[r]*y[jj,r],r=1..s)-add(v[k]*x[jj,k],k=1..de)<=0,j=1..n),

 

And the error message is the following: "Error, invalid sequence". I think that maybe is something really simple like a comma or something like that. Can you please help me? Thank you very much

Hi,

I would like to compute the elements of the permutation group, let us say S10 or S20.

Is there any method to compute all the elements.

And can we make a list of the tranposition and cycles.

Many thinks.

 

Dear Colleagues,

 

1) is there any way to plot the difference of two ode solutions? Smth like display(p1-p2)?
2) is there a way to add a regular plot, i.e. explicit function to the set of plots? I need to have regular plotsand ode plots all together. 

So if you are given a 3x3 matrix i was just wondering how to use plot3d and plot the 3 column vectors given by it?  Thanks.

Hello,

I have a quite long document (around 200 equations + comments).

I need to share my document with collegues.

However, they don't have maple et moreover i would like only to give them my work in a paper format and in a synthetic way. Indeed, they only need to see the main steps of my calculation.

Consequently, i need to make a synthetic presentation of the work i have done on a classic maple worksheet.

For that purpose, i believe that i have 2 ways to do.

1) Create a slideshow. However, i find the mode not very complete in Maple ? Indeed, for the moment, when i have a big section, i didn't manage to scroll on it so as to read the whole section.

2) create a document mode and emphasize only the main equations. 

Questions:

1) What is the best solution for you so as to make a presentation of a work made with classic worksheet?

2) Is it possible to transform a worksheet in a document mode? and if yes, how?

3) Have you some ideas to compact the calculation (or emphasize only the main steps) ?

Thank you for your help

im solving 6 ODE which is the equations are unsteady with boundary conditions.. the program can be run when A=0 but when A=0.2 or others value .. its cannot be run... A means for unsteadiness... before this i solve for steady equations.. this is first time i solve for unsteady using maple.. anyone know where i am wrong??? thanks for helping :)

 

restart; with(plots); n := 2; Ec := 2.0; Pr := .72; N := .2; M := .1; l := 1; Nr := 1; y := 1; blt := 2.5; B := .1; a1 := 1; rho := .5

Eq1 := diff(f(eta), eta, eta, eta)+f(eta)*(diff(f(eta), eta, eta))-(diff(f(eta), eta))^2+l*B*H(eta)*(F(eta)-(diff(f(eta), eta)))-M*(diff(f(eta), eta))-A*(diff(f(eta), eta)+.5*eta*(diff(f(eta), eta, eta))) = 0;

diff(diff(diff(f(eta), eta), eta), eta)+f(eta)*(diff(diff(f(eta), eta), eta))-(diff(f(eta), eta))^2+.1*H(eta)*(F(eta)-(diff(f(eta), eta)))-.1*(diff(f(eta), eta))-A*(diff(f(eta), eta)+.5*eta*(diff(diff(f(eta), eta), eta))) = 0

(1)

Eq2 := A*(F(eta)+.5*eta*(diff(F(eta), eta)))+G(eta)*(diff(F(eta), eta))+F(eta)^2+B*(F(eta)-(diff(f(eta), eta))) = 0;

A*(F(eta)+.5*eta*(diff(F(eta), eta)))+G(eta)*(diff(F(eta), eta))+F(eta)^2+.1*F(eta)-.1*(diff(f(eta), eta)) = 0

(2)

Eq3 := .5*A*(G(eta)+.5*eta*(diff(G(eta), eta)))+G(eta)*(diff(G(eta), eta))+B*(f(eta)+G(eta)) = 0;

.5*A*(G(eta)+.5*eta*(diff(G(eta), eta)))+G(eta)*(diff(G(eta), eta))+.1*f(eta)+.1*G(eta) = 0

(3)

Eq4 := G(eta)*(diff(H(eta), eta))+H(eta)*(diff(G(eta), eta))+F(eta)*H(eta) = 0;

G(eta)*(diff(H(eta), eta))+H(eta)*(diff(G(eta), eta))+F(eta)*H(eta) = 0

(4)

Eq5 := (1+Nr)*(diff(theta(eta), eta, eta))+Pr*((diff(theta(eta), eta))*f(eta)-2*(diff(f(eta), eta))*theta(eta))+N*Pr*a1*(theta1(eta)-theta(eta))/rho+N*Pr*Ec*B*(F(eta)-(diff(f(eta), eta)))^2/rho+Pr*Ec*(diff(f(eta), eta))^2-.5*A*Pr*(4*theta(eta)+eta*(diff(theta(eta), eta))) = 0;

2*(diff(diff(theta(eta), eta), eta))+.72*(diff(theta(eta), eta))*f(eta)-1.44*(diff(f(eta), eta))*theta(eta)+.2880000000*theta1(eta)-.2880000000*theta(eta)+0.5760000000e-1*(F(eta)-(diff(f(eta), eta)))^2+1.440*(diff(f(eta), eta))^2-.360*A*(4*theta(eta)+eta*(diff(theta(eta), eta))) = 0

(5)

Eq6 := 2*F(eta)*theta1(eta)+G(eta)*(diff(theta1(eta), eta))+a1*y*(theta1(eta)-theta(eta))+.5*A*(4*theta1(eta)+eta*(diff(theta1(eta), eta))) = 0;

2*F(eta)*theta1(eta)+G(eta)*(diff(theta1(eta), eta))+theta1(eta)-theta(eta)+.5*A*(4*theta1(eta)+eta*(diff(theta1(eta), eta))) = 0

(6)

bcs1 := f(0) = 0, (D(f))(0) = 1, (D(f))(blt) = 0, F(blt) = 0, G(blt) = -f(blt), H(blt) = n, theta(0) = 1, theta(blt) = 0, theta1(blt) = 0;

f(0) = 0, (D(f))(0) = 1, (D(f))(2.5) = 0, F(2.5) = 0, G(2.5) = -f(2.5), H(2.5) = 2, theta(0) = 1, theta(2.5) = 0, theta1(2.5) = 0

(7)

L := [0., .2, .5];

[0., .2, .5]

(8)

for k to 3 do R := dsolve(eval({Eq1, Eq2, Eq3, Eq4, Eq5, Eq6, bcs1}, A = L[k]), [f(eta), F(eta), G(eta), H(eta), theta(eta), theta1(eta)], numeric, output = listprocedure); Y || k := rhs(R[3]); YP || k := rhs(R[5]); YR || k := rhs(R[6]); YQ || k := rhs(R[7]); YA || k := rhs(R[9]); YB || k := rhs(R[8]) end do

Error, (in dsolve/numeric/bvp) initial Newton iteration is not converging

 

P1 := plot([Y || (1 .. 3)], 0 .. 10, labels = [eta, (D(f))(eta)])

P2 := plot([YP || (1 .. 3)], 0 .. 10, labels = [eta, F(eta)])

plots:-display([P1, P2])

Error, (in plots:-display) expecting plot structures but received: [P1, P2]

 

``

 

Download unsteadyManjunatha.mw

Help!!

Why do the first two of the following 4 examples not work in Maple 15?

subs(m=21,`mod`(m, 4));
subs(m=21,m mod(4));
                               21
                               21
`mod`(21, 4);
21 mod(4);

                               1
                               1
Is there a (simple) workaround?

Thanks.

Hi:
i have two functions(f(x)&g(x)),i calculated values this functions with numerical method(dsolve,numeric),
now i will plot the f(x)-g(x) for 0<x<1,how do this work?

Can someone help me to solve system of equations please. I have a system of 8 complex valued equations, with 8 unknowns: _C1,_C2........_C8

Equation system looks like:

eq_system:={ -3.248046797 10 _C1 + 1.773373463 10 _C2 + (2.182313824 10 - 9.987524076 10 I) _C3 + 1.773373463 10 _C4 = -7.389056097 10 _C2- 7.389056097 10 _C4+ (4.161468365 10 + 9.092974265 10 I)_C3,

............}  its only 1st equation, others are similar.

It looks rather simple though I am not able to solve it with solve or fsolve commands. What I'm doing wrong?

solve(eq_system,{_C1,_C2,_C3,_C4,_C5,_C6,_C7,_C8});

Hello! I wrote a program in Maple but it doesn't work. It crashes with the following error: "Error, (in r_nach_1[3]) too many levels of recursion". I am a newbie in Maple and don't know how to solve this problem. A part of my program that causes the error is given below. Many thanks!

maple.mws


restart;
t1F := 20: A11 := 2.5*10^(-3): c1 := 4*10^2: ss1 := 0.5*10^(-2):
pI1 := 0.4*10^(-1): r1max := 0.6*10^(-2): procent:=0.2:
A12 := procent*A11: B1:=-1: alpha:=0.000001:

beta1:=0.5: r_nach_1[1]:=0: r_nach_1[2]:=0.2*10^(-2):


iter_psi_c1_1:=proc(t) 0 end proc: iter_psi_c2_1:=proc(t) 0 end proc:
s1[1]:=proc(t) 0 end proc:



for i from 1 to 2 do

_p11 := evalf(dsolve({
diff(p1(t), t) = A11*(p1(t)-'s1'[i](t))+B1*(r_nach_1[i+1](t)-r_nach_1[i](t))+A12*(p1(t)-'s1'[i](t))^2+ss1,
p1(0) = pI1}, numeric,method = dverk78, abserr = 1.*10^(-8), relerr = 1.*10^(-8),optimize,output = listprocedure, known=[s1[i],iter_psi_c1_1,iter_psi_c2_1])):

p1F := rhs(_p11(t1F)[2]):
s1[i+1] := subs(_p11, p1(t)):

q11 := evalf(dsolve({diff(q1(t), t) = c1*r_nach_1[i+1](t), q1(0) = 0},numeric,range=0..t1F,known=[iter_psi_c1_1,iter_psi_c2_1])):

q1F := rhs(q11(t1F)[2]):


F[i]:= beta1*q1F;

_psi_c1_1 :=
evalf(dsolve({
diff(psi_c1_1(t), t) = A11+2*A12*('s1'[i+1](t)-'s1'[i](t))*psi_c1_1(t), psi_c1_1(t1F) =  -1+beta1},numeric,method = dverk78, abserr = 1.*10^(-8), relerr = 1.*10^(-8),optimize,known=[s1[i],s1[i+1]],output=listprocedure)):

_psi_c2_1 :=evalf(dsolve({
diff(psi_c2_1(t), t) = 0, psi_c2_1(t1F) =  -beta1},numeric,output=listprocedure)):
iter_psi_c1_1 := op([2,2],_psi_c1_1):
iter_psi_c2_1 := op([2,2],_psi_c2_1):

r_nach_1[i+2] := t-> r_nach_1[i+1](t) +  alpha*(B1*(iter_psi_c1_1(t)) + (c1*iter_psi_c2_1(t))):

_p11 := evalf(dsolve({
diff(p1(t), t) = A11*(p1(t)-'s1'[i+1](t))+B1*(r_nach_1[i+2](t)-r_nach_1[i+1](t))+A12*(p1(t)-'s1'[i+1](t))^2+ss1,
p1(0) = pI1}, numeric,method = dverk78, abserr = 1.*10^(-8), relerr = 1.*10^(-8),optimize, output = listprocedure, known=[s1[i+1],iter_psi_c1_1,iter_psi_c2_1])):

p1F := rhs(_p11(t1F)[2]):

s1[i+2] := subs(_p11, p1(t)):

q11 := evalf(dsolve({diff(q1(t), t) = c1*r_nach_1[i+2](t), q1(0) = 0},numeric,range=0..t1F,known=[iter_psi_c1_1,iter_psi_c2_1])):

q1F := rhs(q11(t1F)[2]):

F[i+1]:= beta1*q1F;


od:
Error, (in r_nach_1[3]) too many levels of recursion

Is it possible to generate a list of values for a function in the form of an ordered pair, like

a cartesian plane, i.e. ( x, f(x) )? And how can I plot a graph with this output?

thanks in advance.

f:=piecewise(x>0,1+x,x=0,2,x<0,1+x^2)

in my work i must use some parts of this piecewise function.

But i don't know how can i call the part that i will use.

How can i do that ?

Thanks for help

 

First 1425 1426 1427 1428 1429 1430 1431 Last Page 1427 of 2427