MaplePrimes Questions

hi.i am a problem with following dsolve.please help me....thanks alot

dsys3 := {10*f2(x)+12*(diff(f1(x), x))+14*f3(x) = 0, 2*(diff(f1(x), x, x))+4*(diff(f2(x), x))+6*(diff(f3(x), x)) = 0, 16*(diff(f3(x), x, x, x, x))+19*(diff(f3(x), x, x))+22*(diff(f1(x), x))+25*f2(x)+27*f3(x)+29*f3(x)+31+32 = 0, f1(0) = 0, f1(1) = 0, f2(0) = 0, f2(1) = 0, f3(0) = 0, f3(1) = 0, ((D@@1)(f1))(0) = 0, ((D@@1)(f1))(1) = 0, ((D@@1)(f2))(0) = 0, ((D@@1)(f2))(1) = 0, ((D@@1)(f3))(0) = 0, ((D@@1)(f3))(1) = 0}; dsol5 := dsolve(dsys3, 'maxmesh' = 500, numeric, range = 0 .. 1, abserr = .1, output = listprocedure); fy3 := eval(f3(x), dsol5); fy2 := eval(f2(x), dsol5); fy1 := eval(f1(x), dsol5)ERROR.mw

ADIABATIC FREE EXPANSION

 

Suppose an ideal gas expands to four times its initial volume. From experience for this process, the initial and final 

temperature are the same.

  1. Using a macroscopic approach, calculate the entropy change for the gas
  2. Using statistical considerations, calculate the change in entropy for the gas and show that it agrees with the answer 
    you obtained in part 1.

What is the probability that the total of two dice will be greater than 9, given that the first die is a 5?

 

Let N be an integer. 

 

For each pair of integers (n,m) where 1<= n,m <= N, we have a variable f_{n,m}(t). 

 

Then for these we have a system of ODEs 

 

d/dt f_{n,m}(t) = \sum_{n', m'} f_m'n' * f_m''n'' * (m'n'' - m''n') 

 

where m''=m-m', n''=n=n', and the sum is simply over for all pairs (n',m'). 

 

I simply do not know how to put these set of equations into Maple in a nice way. 

 

I will really appreciate any help!

  1. Sketch the two vectors listed after the formula for r (t).
  2. Sketch, on the same plane, the curve C dterminated by r (t), and indicate the orientation for the given values of t.

    r (t)=t*i +4*cos (t)j +9*sin (t)k,

    r (0), r(Pi/4),  t>=0

Cause CodeTools:-Test will not eval the paramater:

```
f := x -> x+1;

y := 2;

CodeTools:-Test( f(1), 2); # this will work as normal

CodeTools:-Test( f(1), y); # this will not work because `y` will not eval as 2

```

So current my solution is

```
test := subs(y=2, () -> CodeTools:-Test( f(1), y));

test();

```

 

I don't know is there any more proper solution for this.

More general, is there a way to force evaluate an `uneval` parameter?

 

 

 

 

the final result I want to get should be b*x

Thank you.

How to convert an  algebraic equation to a parametrized one ?

 

How to sketch the curve of a circular cylinder , x^2+y^3^=a^2 with the parametric equations lying on it ?

x=a*cos't) , y=a*sin(t) , z=b*t , if t varies from 0 to 2*Pi, the point P starts at(a,0,0) and moves upward.

How to show that the parametric equations are indeed lying on the cylinder ?

How to plot the whole things ?

 


restart:

m:=2; k:=75.11; L:=3; a=k-m*L; b:=k+m*L-1;

m := 2

k := 75.11

L := 3

a = 69.11

b := 80.11

z:=(k*m)/snr;

z := 150.22/snr

z1:=z^((b+1)/2);

z1 := 0.1893190041e89*(1/snr)^40.55500000

x:=z1/(0.693*GAMMA(m*L)*GAMMA(k));

x := 0.4283082808e-21*(1/snr)^40.55500000

me:=MeijerG([[1-a/2,1+a/2,1+(b+1)/2,1+(b+1)/2],[]],[[1+(b+1)/2],[(b+1)/2]],(1/z));

me := MeijerG([[1.-.5000000000*a, 1.+.5000000000*a, 41.55500000, 41.55500000], []], [[41.55500000], [40.55500000]], 0.6656903209e-2*snr)

cap:=x*me;

cap := 0.4283082808e-21*(1/snr)^40.55500000*MeijerG([[1.-.5000000000*a, 1.+.5000000000*a, 41.55500000, 41.55500000], []], [[41.55500000], [40.55500000]], 0.6656903209e-2*snr)

with(plots):

Warning, the name changecoords has been redefined

plot(cap,snr=1..20);

Warning, unable to evaluate the function to numeric values in the region; see the plotting command's help page to ensure the calling sequence is correct

Error, empty plot

 

 

 

 

Just a begineer in maple and dont know how to rectify this error.

any help..???

Download c_ora.mw

Hi,

I'm using my own Maple package to make questions in Maple T.A.. I'm working with finite automata so I'd like to draw graphs in the text of the questions. I have a function that draws automata, it uses plots:-display as it's output. Is there a way that I can use plots:-display in plotmaple, like this? (This is example doesn't work, I'd like to know why and how I have to modify it.)

$circ=maple("geometry:-circle('Circ', [geometry:-point('PP',0,0),3]);
DRA:=geometry:-draw(Circ);
DISP:=plots:-display(DRA);
DRA;");
$plo=plotmaple("plots:-display($circ),plotdevice='gif', plotoptions='height=250, width=250'");

Thank you for the answers.

test.mw

Hello. As you can see the minimum of set X is 2.

the time pair in set W whose difference is 2 corresponds to elements 54 and 55. How do i locate and print these automatically?

Digits:=4:

U := [1.010, 1.110, 1.210, 1.310, 1.410, 1.510, 1.610, 1.710, 1.810, 1.910, 2.020, 2.120, 2.220, 2.320, 2.420, 2.520, 2.620, 2.720, 2.820, 2.920, 3.030, 3.130, 3.230, 3.330, 3.430, 3.530, 3.630, 3.730, 3.830, 3.930, 4.040, 4.140, 4.240, 4.340, 4.440, 4.540, 4.640, 4.740, 4.840, 4.940, 5.050, 5.150, 5.250, 5.350, 5.450, 5.550, 5.650, 5.750, 5.850, 5.950, 6.060, 6.160, 6.260, 6.360, 6.460, 6.560, 6.660, 6.760, 6.860, 6.960, 7.070, 7.170, 7.270, 7.370, 7.470, 7.570, 7.670, 7.770, 7.870, 7.970, 8.080, 8.180, 8.280, 8.380, 8.480, 8.580, 8.680, 8.780, 8.880, 8.980, 9.090, 9.190, 9.290, 9.390, 9.490, 9.590, 9.690, 9.790, 9.890, 9.990, 10.01, 11.11, 12.21]:

W := select(proc (t) options operator, arrow; .59 >= frac(t) end proc, U);

[1.010, 1.110, 1.210, 1.310, 1.410, 1.510, 2.020, 2.120, 2.220, 2.320, 2.420, 2.520, 3.030, 3.130, 3.230, 3.330, 3.430, 3.530, 4.040, 4.140, 4.240, 4.340, 4.440, 4.540, 5.050, 5.150, 5.250, 5.350, 5.450, 5.550, 6.060, 6.160, 6.260, 6.360, 6.460, 6.560, 7.070, 7.170, 7.270, 7.370, 7.470, 7.570, 8.080, 8.180, 8.280, 8.380, 8.480, 8.580, 9.090, 9.190, 9.290, 9.390, 9.490, 9.590, 10.01, 11.11, 12.21]

(1)

X:=[seq(trunc(W[i+1])*60+frac(W[i+1])*100-(trunc(W[i])*60+frac(W[i])*100),i=1..nops(W)-1)];

[10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 11.00, 10.00, 10.00, 10.00, 10.00, 10.00, 2.00, 70.00, 70.00]

(2)

NULL

min(X),{W[54],W[55]};

2.00, {9.590, 10.01}

(3)

 



Download test.mw

Good afternoon sir.

 

I request your kind suggestion to my cited query.

 

 

With thanks & Regards

 

M.Anand

Assistant Professor in Mathematics

SR International Institute of Technology,

Hyderabad, Andhra Pradesh, INDIA.

 

Need to evaluate above summation in maple. For simplicity infinity can be replaced by some constant say 1000.

The value of parameters used are:

alpha=73.11

beta=76.11

E=4.75

help....

Hi,

I have been trying to solve 2D Diffusion Equation with zero Neumann BC over the unit disk. If I use Gaussian type function with a sharp peak as initial condition, I get huge errors between initial values. Let's say u(r,phi,t) is the solution of the PDE and f(r,phi) is initial value function. The expectation is for the point (r*,phi*) ,  u(r*,phi*,0)=f(r*,phi*), but it is not.

Is Numerical integration in Maple not able to handle such sharp peak? I tried some of the built-in methods such as MonteCarlo,CubaVegas but no difference.

It might be a good idea to specify some nodes arround the peak. There is a command called "peaks", but I could not use it, error message says "invalid arguments".

Thanks in advance.

Soln_2D_Gaussian.mw

I am a research student and I am writing up my thesis. I was reading one of the paper written earlier by authors in 2003.

The authors calculated a symbolic rank of a matrix A, and got 9, using Maple 6.

 

Then new methods developed in 2012 and proved that the "true" rank should be 8. The later paper that "Most  recent versions of Maple have this simplification built in and are able to return a rank of 8.".

 

I just noticed that in Maple 18.01, this matrix was still evaluted for a symbolic rank 9, if no simplifcaiton was done before using Rank(), see attached.

 

I didnt explore a lot, but just as a notice. I am a bit concerned as most of my research was trying to deal with exponentials.

 

Is that something to be fixed in future versions?

 

Thanks!

 

 Bug.mw

First 1365 1366 1367 1368 1369 1370 1371 Last Page 1367 of 2434