how I can convert this maple code to matlab ones?...

how I can convert this maple code to Matlab ones?

 >
 >

Can this apparent contradiction be explained?...

This worksheet displays an intersection between two spheres based on a test which seems unrelated to the display.

How can this be explained?

Intersecting_spheres.mw

evalf/Int issue with Maple 2020...

Hi,

This chunk of code gives me almost instantly the desired result with Maple 2015 but fails returning a result after 5 minutes when ran from  Maple 2020.
Changing the method or trying to tune their options doesn't fix the issue (for a larger value of the absolute tolerance the result, close to 1, is even obviously false).

TIA

 > restart;
 > interface(version)
 (1)
 > with(Statistics):
 > Model_1 := (X, Y, Z) -> 1.3479120558270161+(-1)*.3390713110446567*X+1.0332964944900824*Y+(-1)*.8046331250569234*Z+0.9028899142595827e-1*X^2+(-1)*.28022910356641223*Y^2+1.3698461513469624*Z^2+0.6023630210658247e-1*Y*X+(-1)*.1988111077193626*Z*X+.6782902277463169*Z*Y+(-1)*0.7589729001497135e-1*X*Y*Z: Model_2 := (X, Y, Z) -> .7215754021067852-.961682838202105*X+.4890842364596477*Y-.8214961524335734*Z+.15745942215497866*X^2-1.8023544416693338*Y^2+.36598799463016246*Z^2+1.3957614823018496*Y*X+.725398415577742*Z*X+1.9474707393604542*Z*Y-1.1780201448596594*X*Y*Z: ExpMod_2 := unapply(expand(Model_2(Model_1(X1, Y1, Z1), Y2, Z2)), (X1, Y1, Z1, Y2, Z2)):
 > fy1 := PDF(Normal(0.1055, 0.0297), Y1): fz1 := PDF(Normal(1, 0.2/3), Z1): fy2 := PDF(Normal(0.17, 0.0476), Y2): fz2 := PDF(Normal(1, 0.2/3), Z2): C   := evalf(          (1-eval(CDF(Normal(0.1055, 0.0297), Y), Y=0))          *          (1-eval(CDF(Normal(0.17, 0.0476), Y), Y=0))          *          (eval(CDF(Normal(1, 0.2/3), Z), Z=1.2)-eval(CDF(Normal(1, 0.2/3), Z), Z=0.8))^2        ): J  := Int(         fy1*fz1*fy2*fz2*(1+tanh(10^4*(ExpMod_2(2.14, Y1, Z1, Y2, Z2)-1.25)))/2,         Y1=0..0.4,         Z1=0.8..1.2,         Y2=0..0.4817,         Z2=0.8..1.2,         method = _CubaDivonne, methodoptions=[absepsilon=1e-8]       ): Prob(X__2 > 1.25) = nprintf("%1.3e", evalf(J)/C);
 (2)
 >

"Error while printing" message when trying "Print ...

When trying to print a Maple Document as a PDF on Mac, I am now getting an "Error while printing" message. A 1-page pdf is created that's empty.   This is a new error, one that I've never seen when printing as PDF for any other Maple Document I've created.

I'm running Mac OS 12.3.1, the latest version of Mac OS and using Maple 2020.

Help with assumptions and complex numbers...

I'm trying to solve a coupled inductor problem.  All I'm trying to do is solve for the currents (which I could do), then solve for resonance by setting the imaginary part of I1 to 0 and solving for omega.  I then want to plug this resonant frequency back into I1 and evaluate it for some component values and coupling coefficient.

The problem is that even though I define the various variables as real,  Maple doesn't seem to recognize them as such.  I suspect I've failed to define one of them, and that is why it is responding like it is, but I can't find it.

I'm also getting some weird behavior.....for example if I type :  wres0, I get the value I assigned to it --that's fine.  But then, when I type : wres0/(2*pi), I get an error, "Error, unable to parse 'mverbatim"   --- what does this mean?

This should be a pretty straight forward calculation on Maple, but at this point, I think I may be better off doing it by hand!  I appreciate any help you can provide.  Thank you.

Two questions regarding a sliding ladder...

This worksheet animates part of the motion of the classic ladder sliding down a wall.

Respondents will need to establish their own link to the DirectSearch package

Re: How to treat differentiaion operator as a coef...

Hello all,

Would you allow me to ask this one question?

Is there any possible way to treat the differentiation operator (e.g., 'D') as an ordinary coefficient of a polynomial?

For example, as described in the attached worksheet, I was wondering if there is a way to re-write the LHS of the expression 'eq_e5_9b' as the expression 'desired'.

In Kwon Park

 > restart;
 > with(DEtools):
 > eq_e5_9b := psi__d0*Delta__delta(t) + psi__q0 * D(Delta__delta(t))/omega__0 = p_(Delta__psi__d(t))/omega__0 - Delta__psi__q(t);
 (1)
 > collect(lhs(eq_e5_9b), Delta__delta(t));
 (2)
 > desired := Delta__delta(t) * (psi__d0 + psi__q0 * D/omega__0);
 (3)
 >

Is this O(N^8)?...

I am trying to decide O of times fo a computer process. I thought it was n^2*log(N) or something like that.

With 12  points, it may be hard to determine if O(N^4), etc.

See big-o.mw

Why do these integration results differ?...

The uploaded worksheet attempts to determine arc lengths of an ellipse by two methods. Why do their results differ?

ArcLength_of_Ellipse.mw

Integrating int(1/(A^2 - B^2*cos(t))^(1/2), t = 0 ...

I want to integrate the function int(1/(A^2 - B^2*cos(t))^(1/2), t = 0 .. 2*Pi), I get the correct solution, but I want to prevent all the "if" statements from appearing. I have tried to use "assume" but the if statement still does not disappear.

How to save explore plot ?...

Is there any way to save "explore plot" in such a way that the slider of constants should appear in the extension file?

So that I can see the constant changing in the extension file too.

How do I find the center of mass and moment of ine...

For example for an ellipsoid with principal axes 3, 2 and 1 along the x, y and z coordinate axes respectively, with center at the origin and truncated at x equals 1.

What math will reveal its center of mass and its moment of inertia for an axis of revolution through the center of mass and parallel to the y axis?

Search a list of permutations...

Hi

I wanted to know how to search a list of permutations given criteria: print the permutation(s) If the 4th element of each perm was equal to [C,Z} and the seventh was equal to [B,Y]

perm.mw

how to Rearrange an Equation to a transfer functio...

I tried to simplify this equation:

where Ts,Cf,Rf,Lg & L are Variables to equivelent equation in the form of:

(.....)Z3+(....)Z2+(.....)Z+(.........)      /   (.....)Z3+(....)Z2+(.....)Z+(.........)

please (i need code to calculate brownian motion ...

Hi professors,

i want code of  function in maple to calculate values  of B(t),

i have this matlab code

randn('state',100)
T = 1; N = 5; dt = T/N;

dW = sqrt(dt)*randn(1,N);
W = cumsum(dW)

but does not calculate B(0.25) and B(0.5)

Thanks,

 5 6 7 8 9 10 11 Last Page 7 of 56
﻿