Hello,

I would like to plot coupled elliptic parametric curves.

I have a problem in the for loop in my code. The index i of the phi variable is not incremented.

Here my code:

for i to 4
do
u[i]:=t->L/2*cos(w*t);
v[i]:=t->H*sin(w*t+phi[i]);
end do;

The result looks like this :

u[1](t) -> - L cos(w t)
v[1](t) -> H sin(w t + phi[i])

u[2](t) -> - L cos(w t)
v[2](t) -> H sin(w t + phi[i])

u[3](t) -> - L cos(w t)
u[3](t) -> H sin(w t + phi[i])

u[4](t)  -> - L cos(w t)
u[4](t)  -> H sin(w t + phi[i])

And i would have prefered :

u[1](t) -> - L cos(w t)
v[1](t) -> H sin(w t + phi[1])

u[2](t) -> - L cos(w t)
v[2](t) -> H sin(w t + phi[2])

u[3](t) -> - L cos(w t)
v[3](t) -> H sin(w t + phi[3])

u[4](t)  -> - L cos(w t)
v[4](t)  -> H sin(w t + phi[4])

Do you have ideas how I can obtain iteration on my phi variable ?

Thank you for your help.

Can some one help me " How to conert a three or four coupled equation by using non dimensions variable "

Hi,

I have an issue calculating an electronics circuit with Maple, using units. I have a current source that I know, and I want to determin the voltage in a capacitor by solving an ODE (except that the current source is defined piecewise). And to make sure I have all the units and scales right, I use the standard unit package. All my variables have their units defined.

Except that Maple doesn't want to solve the equation. It seems to me that it assumes that the function I am trying to solve is unitless, and therefore refuses to solve. 

V__out := 3*Unit('kV');

C__out := 2*Unit('nF');
R__blead := 520*Unit('`kΩ`');

I__fly := proc (t) options operator, arrow; Unit('A')*piecewise(t < 3.25*Unit('us'), (1+(-1)*t/(3.25*Unit('us')))*.2, 0) end proc;

dsolve({I__fly(t*Unit('s'))-V__C(t*Unit('s'))/R__blead = C__out*(diff(V__C(t*Unit('s')), t)), V__C(0*Unit('s')) = V__out}, V__C(t*Unit('s')));
Error, (in Units:-Standard:-+) the units `A` and `S` have incompatible dimensions

Is there a way to make Maple assume the unit of what it's trying to solve ? I need it to understands that V__C is in Unit('V') ...

Thanks

Hi,

what i have to calculate in Maple for q=Mises-Stress and the third invariant of deviatoric stress as r???

It is for:

Link:
http://www.egr.msu.edu/software/abaqus/Documentation/docs/v6.7/books/stm/default.htm?startat=ch04s04ath116.html

hi...amount of Determinant  is infinity?how i can remove this bad calculation ?

thanks...mode_shape2.mw

hi .by changing amount of variable such as '' f '' or 'D5 ' in input data, not changing  in final result!!!!please chek attached file 

thanksX4.mw

Hello,

I try his code:

for i to 6 do u[i, j] := (u[i-1, j]+u[i+1, j]+lambda*(u[i, j-1]+u[i, j+1])+sigma(h, r[i+1])*(u[i+1, j]-u[i-1, j]))/mu end do;

 
Error, too many levels of recursion.

I want to the  results of the analytic equations, for example:

                                                          
 u[1,j]= (u[0, j] + u[2, j] + lambda (u[1, j - 1] + u[1, j + 1])                                                       

      + sigma(h, r[2]) (u[2, j] - u[0, j]))/mu

.

.

.

.
Regards.

Suppose that we have the following curve on a unit sphere

r^2+z^2=1

phi=pi/3

-1<=z<=1

which is the implicit equation of a circle in "cylindrical coorinates" (r,phi,z) on the sphere.

How can I plot this curve without solving the implicit equation for r or z ?

Also, I don't want to make any parameterization.

See attached file and code

0. This is the differential equation I'm trying to do:

http://www.intmath.com/differential-equations/6-rc-circuits.php

https://i.imgur.com/zlVIisR.png

1. After you look at the above imgur link, could you help me with this error 

Error, (in Units:-Standard:-+) the units ``&Omega;`` and `1/F` have incompatible dimensions

2. Why does my solve ODE fail? 

See my code: 

test_maple.mw

Spiral on the cone. 







Yes, of course, in Maple.  The same source

I have Maple package "DESOLVE" having two different types of files namely in format .mpl and .mws, where code is written in .mpl file and .mws is just a example worksheet, this package is available at following link;

http://cpc.cs.qub.ac.uk/summaries/ADYZ_v1_0.html

My problem is how to work with package?

One thing that I able to do is that I have opened the file type .mpl  in Maple worksheet and after executing whole whorsheet I able to to start working with package using with(desolv).

But when I type with(desolv) in new worksheet an error shows up.

How I can store this .mpl file in Maple library so as to use this package using "with(desolv)" in new worksheets

Desolv_worksheet_paper.mws?

Hello,

Remark by Markiyan Hirnyk:

The part of the body of this question disappears after my edit of its title only.

Consider the well-known Euler's formula  

 e^ix = cos x + i sin x   

When we calculate that for  x = π  we get:

e^iπ = cos π + i sin π   or

e^iπ = −1 + i × 0   (because cos π = −1 and sin π = 0)  or

e^iπ = −1  or  e^iπ + 1 = 0

It seems absolutely magical that such a neat equation combines  5  fundamental constants: e ,  i ,  π , 1 , 0

The purpose of this post - to give a simple visualization of equality  e^iπ = −1  (statical and animated) by expanding  e^iπ  in a series of complex numbers. These numbers we represent as vectors in the plane. We will see that the partial sums of this series are broken lines like a spiral, twisting around the point -1 steadily approaching to it.

Euler procedure has one required parameter  n is positive integer - the number of displayed terms of the series  for  e^iπ  

Optional parameter  a  is any symbol (by default  a=NULL). We use this option if  instead of a static spiral want to see an animated spiral. 

Procedure code can be found in the attached file  Euler.mw

Examples of use.

The first example shows  8 terms of the series (broken line of 8 units):

Euler(8);

The terms of the series where  n> = 10  on the same plot can not be seen as very small. In this case, we use  the second plot with magnification of  100 : 1 .  

The second example:

Euler(14);

In the third example, we see an animated broken line. It's  first 9 units represented  on the left plot, and then for n> = 10 on the right plot:

Euler(13, a);

Euler.mw

I want to have a matrix such that 

Crosslinking surfaces by a spiral.



My fantasy is a source.