MaplePrimes Questions

Search Questions:

Latest Questions Latest Questions Feed

Hy i need your help experts,

I have à small program for some variables calculations and im intrested to plot loop result,

Here is the program:

 



I just want to assign the results of my loop to a matrix to be easily ploted.

Thak you for help.

 

Dear all;

Please how can I plot the error between the two function.

 

Hi:

can I use Laplace Transform Method for solve Nonlinear ODE in maple?

This has been buggin me all day, but if i define a variable with an expresion say "a:=3 x+5" and i then want to assign this expression to a function like so "f:=x->a". if i then call f(7), maple return with "3 x+5" so just returning the value of a without substituting x with 7.

Is there any way in which i can define a function this way?

Hello everybody,

Since I am geting older, I need to fix the font size form 12 to 14.  I have been tryong to find where to do that permanent change.  But without success.  Could someone help to to resolve this problem.  Thank you in advance for your help.

 

--------------------------------------
Mario Lemelin
Maple 18 Ubuntu 13.10 - 64 bits
Maple 18 Win 7 - 64 bits messagerie : mario.lemelin@cgocable.ca téléphone :  (819) 376-0987

I have a function. It's calculating coefficients of a polynomial.

So when it's done, I have the variables a0, a1, a2,a3 ... aN

where N is an input value and could be any positive integer.

 

Now, I just want to show them on the screen. So if I have, say , a0=1.3 and a1=5.76, I want this to appear:

1.3x + 5.76x2

 

That's all. I just want a loop that does this:

for term from 1 to N do:

    print( a0 ); print (x^N); print ("+");

od:

so all the terms come out on the same line, as in my example above. I've tried so many ways, but each time it ignores the numerical values of a0 and just writes "a0" on the screen, or breaks it up into a set, or puts || marks everywhere, or does it on sequential lines, or something or something or something.

I have Maple 16 for Windows 64-bit. How to get a printable Maple 16 Quick Reference Card? Also, how to save the Maple 16 Quick Reference Card as a separate file? (When I try to copy the Card I get en error message "selection too large").

 

How can I use separation of variable to solve heat equation ut=uxx on a rod of length pi/4 with u(0,t)=0, u(pi/4)=0 and u(x,0)=x(pi/4-x) 

then solve heat equation by laplace transform

Dear Maple experts,

I would like to generate population data that is the best possible approximation of a multivariate normal distribution with a specified covariance matrix and vector of means. I do not want to draw a sample from a multivariate distribution, but I want the population values itself which are approximately multivariate normal distributed. The size of the datamatrix should be limited, otherwise I could draw a huge sample from a multivariate normal distribution. For instance, I would like to generate a 200 by 6 data matrix that is the best (or at least good enough) approximation of a MVN distribution. For a bivariate normal distribution one could calculate the probalities of a grid by integrating the density, but for six variables that seems undoable. 

Before trying the invent the wheel again, I think I will ask this question to experts, because it is unlikely that there is no already existing algorithm that does the job pretty well.

Thanks in advance,

Harry Garst

i have a list of data, i can plot them with point plot,but i can not plot it in line style not point, could anyone help ? tnx for your help. 

restart:m1:=3:m2:=3:m3:=3:

P:=unapply(Vector(9, {(1) = 0, (2) = 0, (3) = 0, (4) = -50, (5) = -100/3, (6) = -25, (7) = 0, (8) = 0, (9) = 0}),t):

K:=unapply(Matrix(9, 9, {(1, 1) = -200000000, (1, 2) = -600000000, (1, 3) = -1000000000, (1, 4) = -100000000*a[4](t)-600000000*a[5](t)-1000000000*a[6](t), (1, 5) = -2000000000/3*a[5](t)-2100000000*a[6](t), (1, 6) = -1620000000*a[6](t), (1, 7) = 0, (1, 8) = 0, (1, 9) = 0, (2, 1) = -200000000, (2, 2) = -1600000000/3, (2, 3) = -900000000, (2, 4) = -100000000*a[4](t)-1600000000/3*a[5](t)-900000000*a[6](t), (2, 5) = -600000000*a[5](t)-1920000000*a[6](t), (2, 6) = -1500000000*a[6](t), (2, 7) = 0, (2, 8) = 0, (2, 9) = 0, (3, 1) = -200000000, (3, 2) = -500000000, (3, 3) = -840000000, (3, 4) = -100000000*a[4](t)-500000000*a[5](t)-840000000*a[6](t), (3, 5) = -560000000*a[5](t)-1800000000*a[6](t), (3, 6) = -9900000000/7*a[6](t), (3, 7) = 0, (3, 8) = 0, (3, 9) = 0, (4, 1) = -200000000*a[4](t)-600000000*a[5](t)-1000000000*a[6](t), (4, 2) = -600000000*a[4](t)-4000000000/3*a[5](t)-2100000000*a[6](t), (4, 3) = -1000000000*a[4](t)-2100000000*a[5](t)-3240000000*a[6](t), (4, 4) = -100000000*a[4](t)^2+(-900000000*a[5](t)-1500000000*a[6](t))*a[4](t)-6300000000*a[5](t)*a[6](t)-4860000000*a[6](t)^2-175000000/3-2000000000*a[5](t)^2, (4, 5) = -1400000000*a[5](t)^2-6480000000*a[5](t)*a[6](t)-9900000000*a[6](t)^2-175000000, (4, 6) = -35100000000/7*(a[6](t)^2)-875000000/3, (4, 7) = -87500000, (4, 8) = -875000000/9, (4, 9) = -306250000/3, (5, 1) = -200000000*a[4](t)-1600000000/3*a[5](t)-900000000*a[6](t), (5, 2) = -1600000000/3*a[4](t)-1200000000*a[5](t)-1920000000*a[6](t), (5, 3) = -900000000*a[4](t)-1920000000*a[5](t)-3000000000*a[6](t), (5, 4) = -100000000*a[4](t)^2+(-800000000*a[5](t)-1350000000*a[6](t))*a[4](t)-5760000000*a[5](t)*a[6](t)-4500000000*a[6](t)^2-175000000/3-1800000000*a[5](t)^2, (5, 5) = -1280000000*a[5](t)^2-6000000000*a[5](t)*a[6](t)-64800000000/7*(a[6](t)^2)-1400000000/9, (5, 6) = -4725000000*a[6](t)^2-262500000, (5, 7) = -700000000/9, (5, 8) = -87500000, (5, 9) = -280000000/3, (6, 1) = -200000000*a[4](t)-500000000*a[5](t)-840000000*a[6](t), (6, 2) = -500000000*a[4](t)-1120000000*a[5](t)-1800000000*a[6](t), (6, 3) = -840000000*a[4](t)-1800000000*a[5](t)-19800000000/7*a[6](t), (6, 4) = -100000000*a[4](t)^2+(-750000000*a[5](t)-1260000000*a[6](t))*a[4](t)-5400000000*a[5](t)*a[6](t)-29700000000/7*(a[6](t)^2)-175000000/3-1680000000*a[5](t)^2, (6, 5) = -1200000000*a[5](t)^2-39600000000/7*a[5](t)*a[6](t)-8775000000*a[6](t)^2-437500000/3, (6, 6) = -4500000000*a[6](t)^2-245000000, (6, 7) = -218750000/3, (6, 8) = -245000000/3, (6, 9) = -87500000, (7, 1) = 0, (7, 2) = 0, (7, 3) = 0, (7, 4) = 87500000/3, (7, 5) = 350000000/9, (7, 6) = 43750000, (7, 7) = 173500000/9, (7, 8) = 42250000/3, (7, 9) = 32500000/3, (8, 1) = 0, (8, 2) = 0, (8, 3) = 0, (8, 4) = 175000000/9, (8, 5) = 87500000/3, (8, 6) = 35000000, (8, 7) = 43250000/3, (8, 8) = 101000000/9, (8, 9) = 80750000/9, (9, 1) = 0, (9, 2) = 0, (9, 3) = 0, (9, 4) = 43750000/3, (9, 5) = 70000000/3, (9, 6) = 87500000/3, (9, 7) = 11500000, (9, 8) = 83750000/9, (9, 9) = 22900000/3}),t):

M:=Matrix(9, 9, {(1, 1) = 13/5, (1, 2) = 39/20, (1, 3) = 39/25, (1, 4) = 0, (1, 5) = 0, (1, 6) = 0, (1, 7) = 0, (1, 8) = 0, (1, 9) = 0, (2, 1) = 39/20, (2, 2) = 39/25, (2, 3) = 13/10, (2, 4) = 0, (2, 5) = 0, (2, 6) = 0, (2, 7) = 0, (2, 8) = 0, (2, 9) = 0, (3, 1) = 39/25, (3, 2) = 13/10, (3, 3) = 39/35, (3, 4) = 0, (3, 5) = 0, (3, 6) = 0, (3, 7) = 0, (3, 8) = 0, (3, 9) = 0, (4, 1) = 0, (4, 2) = 0, (4, 3) = 0, (4, 4) = 260, (4, 5) = 195, (4, 6) = 156, (4, 7) = 0, (4, 8) = 0, (4, 9) = 0, (5, 1) = 0, (5, 2) = 0, (5, 3) = 0, (5, 4) = 195, (5, 5) = 156, (5, 6) = 130, (5, 7) = 0, (5, 8) = 0, (5, 9) = 0, (6, 1) = 0, (6, 2) = 0, (6, 3) = 0, (6, 4) = 156, (6, 5) = 130, (6, 6) = 780/7, (6, 7) = 0, (6, 8) = 0, (6, 9) = 0, (7, 1) = 0, (7, 2) = 0, (7, 3) = 0, (7, 4) = 0, (7, 5) = 0, (7, 6) = 0, (7, 7) = 13/6000, (7, 8) = 13/8000, (7, 9) = 13/10000, (8, 1) = 0, (8, 2) = 0, (8, 3) = 0, (8, 4) = 0, (8, 5) = 0, (8, 6) = 0, (8, 7) = 13/8000, (8, 8) = 13/10000, (8, 9) = 13/12000, (9, 1) = 0, (9, 2) = 0, (9, 3) = 0, (9, 4) = 0, (9, 5) = 0, (9, 6) = 0, (9, 7) = 13/10000, (9, 8) = 13/12000, (9, 9) = 13/14000}):

 

 

 

 

 

 

 

 

 

 

V:=unapply(Vector(9,i->a[i](t)),t):

Vdot:=unapply(Vector(9,i->dd1[i](t)),t):

V2dot:=unapply(Vector(9,i->dd2[i](t)),t):

 

 

 

 

for i to 9 do
a[i](0):=0
od:

for i to 9 do
dd1[i](0):=0:
od:

V2dot(0):=M^(-1).(P(0)-K(0).V(0)):

 

 

#K(t)(8..10,8..10);

 

 

#Newmark step by step SOLVER #

 

 

h:=1/100;

1/100

(1)

t:=0;

0

(2)

hh:=0.5;

.5

(3)

for i from 0 to 0.5 by h do


KC:=K(t)+(4/h^2)*M:
PC:=(P(t+h)-P(t))+M.((4/h)*Vdot(t)+2*V2dot(t)):
Delta(V):=(KC^(-1)).PC:
V(t+h):=Delta(V)+V(t):
Delta(Vdot):=(2/h)*Delta(V)-2*Vdot(t):
Vdot(t+h):=Delta(Vdot)+Vdot(t):
Delta(V2dot):=M^(-1).((P(t+h)-P(t))-(K(t).Delta(V))):
V2dot(t+h):=Delta(V2dot)+V2dot(t):
for k from 1 to 9 do
a[k](t+h):=evalf(V(t+h)[k])
od:

u[t+h]:=evalf(x*(add(a[m](t+h)*x^(m-1),m=1..m1))-h/2*x*(add(a[m1+m3+m](t)*x^(m-1),m=1..m2))):
#plot(u[t+h],x=0..L,title="displacement u for time" ,axes = framed,legend = ["displacement u"],color="NavyBlue"):


w[t+h]:=evalf(x*(add(a[m1+m](t+h)*x^(m-1),m=1..m3))):
#plot(w[t+h],x=0..L,title="displacement w for time  ",axes = framed,legend = ["displacement w"],color="red"):

#phi[t+h]:=evalf(x*(add(a[m1+m3+m](t)*x^(m-1),m=1..m2)));
#plot(phi[t+h],x=0..L,title="displacement phi for time  ",axes = framed,legend = ["displacement w"],color="green");

t:=t+h:
x:='x';
od:

t:='t';w[0]:=0;u[0]:=0;

t

 

0

 

0

(4)

x:=1;

1

(5)

plots:-pointplot({seq([t,w[t]],t=0..0.3,h)});

 

``

NULL


Download plot.mw

okay so we have to, In the Notation of this question http://www.mapleprimes.com/questions/201278-Fix-A-Syntax-Error-In-My-Simple-Function-please-Help , let f0 be contained in V be any given function and de ne
a sequence (fn) with n contained in (N subscript 0) of functions fn contained in V by

(f subscript 0) := f0  and   fn+1=Afn  for all n contained in (N subscript 0).


prove that this sequence converges pointwise, i.e. that for all   i,j contained in [N]X[N],
(f subscript infinity)(i, j) := (lim n-> infinity  fn(i, j)) exists and that "triangle"(f subscript infintiy) =0
and that
f1 = 0 :

any help would be much appreciated guys <3  

I am using the following commands to create an animation of a simple cardioid being traced out. I would like to add to this an animation of the filling of the region bounded by the curve and the coordinate axes as the curve is being traced out. How might one achieve this?

with(plots):

animatecurve([2-2*sin(theta),theta,theta=0.. 2*Pi],coords=polar,axiscoordinates=polar,frames=250, numpoints=300, scaling=constrained)

Thanks.

I have a plot of a sinusoidal function(cos_phi) which changes with incrementing values of variable k. When this function goes above 1 or below -1 I would like to have the range of k for which it occurs outputed somehow. Would anybody know how to do this?

I am having trouble printing out a limit cylce on maple 16.  I have the attached file and if anybody could look at it and perhaps help me out it would be greatly appreciated.  The first limit cycle is supposed to look somewhat like the second one.  I'v tried many different things but nothing seems to be working.  an explenation would also be nice too.  if the file does not open correctly also let me know. thank you very much.  

 Math_4710_Hilbert_16.mw

Say I have this list,

tmp:=[[0, 1, 2], [1, 0, 2], [1, 1, 2], [1, 2, 0]];

and the sums of each element (list),

map(x->convert(x,`+`),tmp);

 

How do I quickly pickup the elements, where they sum to 3? Like this:

wanttohave:=[[0, 1, 2], [1, 0, 2],  [1, 2, 0]];

Thanks,

4 5 6 7 8 9 10 Last Page 6 of 992