Hi,

I want to parallel the following code.it is good that part1 solved by one thread and part2 solved by other thread.

please help me.it must be noted that part1 and part 2 are independent "at each" iteration.

dtheta[m](x):= theta[m-1](x)+phi[m-1](x);
dphi[m](x):= theta[m-1](x)-2*phi[m-1](x);

for i from 2 by 1 to 10 do;

# part1

theta[i](x):=subs(m=i,dtheta[m](x)):
theta[i](x):=value(%):
eval(-theta[i](x),x=0);

Question :

Write a procedure that determines the solutions of a quadratic equation from inputs ,  and  by using the discriminant and the quadratic formula.

The quadratic equation procedure should be able to solve all cases: invalid input, linear case, real and complex roots. The procedure should also plot the given equation.

I dont understand how to do this can someone please please do it this for me urgently required.

today, i stuty a differential equation,  firsy, i conver the equation from cartesian to polar as  follow:

restart: 

sys := diff(x(t), t) = -y(t)+x(t)*(1-2*x(t)^2-3*y(t)^2), diff(y(t), t) = x(t)+y(t)*(1-2*x(t)^2-3*y(t)^2):

then  using transformation x=r*cos(theta),y=r*sin(theta) convert it to polar coordinate,and get:

sys1:=diff(r(t), t) = r(t)*(1-2*r(t)^2-r(t)^2*sin(theta(t))^2), diff(theta(t), t) = 1

then  plot the phase portrait, i use to method ,but all of them failed:

first ,use DEplot

  with(DEtools): 

 DEplot(sys1, [theta(t), r(t)], t = 0 .. 10, [[theta(0) = 2, r(0) = 2]], coords = polar);

and get error:
Error, (in DEtools/DEplot) can only plot in cartesian co-ordinates, got invalid option coords = polar

second,use odeplot

   with(plots):

> q1 := dsolve([diff(r(t), t) = r(t)*(1-2*r(t)^2-r(t)^2*sin(theta(t))^2), diff(theta(t), t) = 1, theta(0) = Pi, r(0) = 2], numeric, [theta(t), r(t)]);

> with(plots);

> odeplot(q1, [r(t), theta(t)], t = 0 .. 1, numpoints = 500, axiscoordinates = polar, coords = polar);

and also get a wrong result, so how to plot the phase portrait in polar coordinate using DEplot or odeplot

How implementing in Maple this MATLAB's sequence?

timing_data = 0:0.01:5;
X = timing_data; %(s)
delta_t = mean(diff(X))

delta_t =

    0.0100

Gracias

i want to write a program to simulate a walk of a drunken sailor who walks out of a bar at the centre of a city. the city has a regular pattern of 24\24 square blocks. the sailor doesnt know where he is going and how to exist teh city. the question is how many blocks will he have to walk to exit the city?

let him start at the centre of the city and assume that he exists the city when he reaches either north, east, south or western borders. it can happen that he wlks...

Is someone else noticing that the last ubuntu update (12.04 64 bits) is crashing some maple 13 commands (like asympt)? In the 32 bits version of 12.04 and in ubuntu 11.10 or windows 7 (both 32 and 64 bits) commands like "asympt(3x/2 + 3/2, x)" works fine, but after update it fails...

When was the first versionof MAPLE released?  I shall appreciate if someone answers this question.  Thanks.

Kodanat B.M. Nambudiripad

hi, i want to write the following algorithm in maple with two threads, please help me. for m from 1 by 1 to 100 if thread1 then result1:= some equation if thread2 then result2:= some equation end if end do:

Hi, I have an equation that i need to integrate. Howver, due to the absolute equations I have in it, an analytical integration is not possible. How could i use a numerical method instead?

In matlab, instead of using "int", i use "quad". Is there such a thing in maple?

Please advice. Thanks

I have a problem with inputting a new package in my version of maple (15) for mac. The packes makes it possible to do calculations with cos/sin/tan in degrees instead of radians. I need to use this for my Engeneering studies in Denmark.

In pc you just place the package in the labrary subfolder under the Maple 15 folder, but i can't find that subfolder anywhere on my mac??
Where should i put the package to make it work?

I have 2 files; "Trig.ind" and "Trig.lib" - whichone to use?

Hi everyone,

I have a problem with DEplot.

Here is the code:

restart;with(DEtools), with(plots), with(linalg):

z:=exp(I*((t/(2*Pi))+(2*Pi*j/3)));

f1:=Re((I/(2*Pi))((1/(subs(j=1,z)-(x(t)+I*y(t))))+(1/(subs(j=2,z)-(x(t)+I*y(t))))+(1/(subs(j=3,z)-(x(t)+I*y(t))))));

g1:=-1*Im((I/(2*Pi))((1/(subs(j=1,z)-(x(t)+I*y(t))))+(1/(subs(j=2,z)-(x(t)+I*y(t))))+(1/(subs(j=3,z)-(x(t)+I*y(t))))));

sys1:=[diff(x(t),t)=f1,diff(y(t),t)=g1]:

Is there a command (possibly in a package) which can help in reading parameters from ini text files like

    [section1]

    option1=1.0

    option2=2.0

    [section2]

    option3=3.0

Thanks in advance!

A := Matrix([[5, 4, 5, 4, 8, 4, 6, 1, 4], [2, 8, 4, 7, 5, 1, 9, 6, 1], [8, 4, 9, 8, 8, 8, 6, 2, 8], [7, 1, 8, 6, 6, 5, 3, 7, 7], [8, 7, 9, 4, 7, 8, 2, 7, 4], [5, 1, 7, 8, 5, 7, 2, 2, 8], [7, 4, 1, 7, 7, 1, 7, 8, 6], [2, 9, 9, 4, 7, 5, 2, 6, 8]]);

Create the vector b that is column 7 from A and the matrix C that is made from columns 1 to 6 and 8 to 9 of A (in the same order as the columns of A

E := proc (phi) options operator, arrow; int(1/(1-(1-1/kappa^2)*sin(phi)^2)^(1/2), phi = 0 .. phi) end proc

E(pi/2)

Is there a way to get this function to work properly without having to simplify the integral manually ?

kappa is a constant (with a value of 12.4)

Thanks in advance!

There is a utilde on the Latin pallette, among others. How do I get xtilde as a variable?