Attached is a photo with the code I am working for.  

On the top is practice code with a simpler ODE to help with trouble shooting, on the bottom is the ODE I am working with.

I was hoping to gain insight about the _z1 symbol in the solution, I haven't been able to find much help on other threads.  I would like to know how I can go about working with it - if it is something on my end or if it is the nature of the equation I am working with.

Thank you for any help,

Josh

Hello there

I'm quite an amature so please don't judge.  I'm trying to use fsolve to solve a system of non-linear equations but Maple is just "spitting" on me the equations with no intention to solve them:

> delta5 := P*(1+mu5)*((1-2*mu5)*x/(sqrt(x^2+zeq^2)*(sqrt(x^2+zeq^2)*x))+x*zeq/sqrt(x^2+zeq^2)^3)/(2*Pi*E5);
print(`output redirected...`); # input placeholder
> shrinkage := P*(1+mu5)*((1-2*mu5)*x/(sqrt(x^2+Zb^2)*(sqrt(x^2+Zb^2)*x))+x*Zb/sqrt(x^2+Zb^2)^3)/(2*Pi*E5)-P*(1+mu5)*((1-2*mu5)*x/(sqrt(x^2+Za^2)*(sqrt(x^2+Za^2)*x))+x*Za/sqrt(x^2+Za^2)^3)/(2*Pi*E5);
> eq10 := subs(x = 1800, delta5)+subs(x = 1800, Zb = z2, Za = z1, shrinkage)+subs(x = 1800, Zb = z3, Za = z2, shrinkage)+subs(x = 1800, Zb = z4, Za = z3, shrinkage)+subs(x = 1800, Zb = z5, Za = z4, shrinkage) = 36.7*10^(-3);
print(`output redirected...`); # input placeholder
> eq9 := subs(x = 1500, delta5)+subs(x = 1500, Zb = z2, Za = z1, shrinkage)+subs(x = 1500, Zb = z3, Za = z2, shrinkage)+subs(x = 1500, Zb = z4, Za = z3, shrinkage)+subs(x = 1500, Zb = z5, Za = z4, shrinkage) = 47.2*10^(-3);
print(`output redirected...`); # input placeholder
> eq8 := subs(x = 1200, delta5)+subs(x = 1200, Zb = z2, Za = z1, shrinkage)+subs(x = 1200, Zb = z3, Za = z2, shrinkage)+subs(x = 1200, Zb = z4, Za = z3, shrinkage)+subs(x = 1200, Zb = z5, Za = z4, shrinkage) = 63.8*10^(-3);
> eq7 := subs(x = 900, delta5)+subs(x = 900, Zb = z2, Za = z1, shrinkage)+subs(x = 900, Zb = z3, Za = z2, shrinkage)+subs(x = 900, Zb = z4, Za = z3, shrinkage)+subs(x = 900, Zb = z5, Za = z4, shrinkage) = 91.1*10^(-3);
print(`output redirected...`); # input placeholder
> eq6 := subs(x = 600, delta5)+subs(x = 600, Zb = z2, Za = z1, shrinkage)+subs(x = 600, Zb = z3, Za = z2, shrinkage)+subs(x = 600, Zb = z4, Za = z3, shrinkage)+subs(x = 600, Zb = z5, Za = z4, shrinkage) = 137.9*10^(-3);
> eq5 := subs(x = 450, delta5)+subs(x = 450, Zb = z2, Za = z1, shrinkage)+subs(x = 450, Zb = z3, Za = z2, shrinkage)+subs(x = 450, Zb = z4, Za = z3, shrinkage)+subs(x = 450, Zb = z5, Za = z4, shrinkage) = 175.2*10^(-3);
> eq4 := subs(x = 300, delta5)+subs(x = 300, Zb = z2, Za = z1, shrinkage)+subs(x = 300, Zb = z3, Za = z2, shrinkage)+subs(x = 300, Zb = z4, Za = z3, shrinkage)+subs(x = 300, Zb = z5, Za = z4, shrinkage) = 230.9*10^(-3);
print(`output redirected...`); # input placeholder
> sys := {eq10, eq5, eq6, eq7, eq8, eq9};
print(`output redirected...`); # input placeholder
> fsolve(sys, {E1 = 1000 .. 2000, E2 = 0 .. 2000, E3 = 0 .. 2000, E4 = 0 .. 2000, E5 = 0 .. 2000, h4 = 100 .. 400});

and this is what Maple gives after the fsolve

fsolve({(3937.500000*(.2/(202500+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(450*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(202500+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.3888888889e-2/E5+(3937.500000*(.2/(202500+(650+h4)^2)+(450*(650+h4))/(202500+(650+h4)^2)^(3/2)))/E5 = .1752000000, (3937.500000*(.2/(360000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(600*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(360000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.2187500000e-2/E5+(3937.500000*(.2/(360000+(650+h4)^2)+(600*(650+h4))/(360000+(650+h4)^2)^(3/2)))/E5 = .1379000000, (3937.500000*(.2/(810000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(900*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(810000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.9722222220e-3/E5+(3937.500000*(.2/(810000+(650+h4)^2)+(900*(650+h4))/(810000+(650+h4)^2)^(3/2)))/E5 = 0.9110000000e-1, (3937.500000*(.2/(1440000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(1200*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(1440000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.5468750000e-3/E5+(3937.500000*(.2/(1440000+(650+h4)^2)+(1200*(650+h4))/(1440000+(650+h4)^2)^(3/2)))/E5 = 0.6380000000e-1, (3937.500000*(.2/(2250000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(1500*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(2250000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.3500000000e-3/E5+(3937.500000*(.2/(2250000+(650+h4)^2)+(1500*(650+h4))/(2250000+(650+h4)^2)^(3/2)))/E5 = 0.4720000000e-1, (3937.500000*(.2/(3240000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)+(1800*(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3)))/(3240000+(146.0507832*(E1/E5)^(1/3)+197.1094212*(E2/E5)^(1/3)+295.6641318*(E3/E5)^(1/3)+1.*h4*(E4/E5)^(1/3))^2)^(3/2)))/E5-0.2430555555e-3/E5+(3937.500000*(.2/(3240000+(650+h4)^2)+(1800*(650+h4))/(3240000+(650+h4)^2)^(3/2)))/E5 = 0.3670000000e-1}, {E1, E2, E3, E4, E5, h4}, {E1 = 1000 .. 2000, E2 = 0 .. 2000, E3 = 0 .. 2000, E4 = 0 .. 2000, E5 = 0 .. 2000, h4 = 100 .. 400})

I have a nice procedure that is as follows, that explores the behaviour of a mapping in R^2.

InversePoincare := proc (v) options operator, arrow; v*Norm(v, 2)/(1-Norm(v, 2)^2) end proc

SphereVectorPlot := proc (T, radius, radiusIncrements, angleIncrements, lengthParam)
local listOfPairs;
listOfPairs := {seq(seq([[radius*r*cos(2*angle*Pi/angleIncrements)/radiusIncrements, radius*r*sin(2*angle*Pi/angleIncrements)/radiusIncrements], T(InversePoincare(`<,>`(radius*r*cos(2*angle*Pi/angleIncrements)/radiusIncrements, radius*r*sin(2*angle*Pi/angleIncrements)/radiusIncrements)))-InversePoincare(`<,>`(radius*r*cos(2*angle*Pi/angleIncrements)/radiusIncrements, radius*r*sin(2*angle*Pi/angleIncrements)/radiusIncrements))], r = 0 .. radiusIncrements), angle = 0 .. angleIncrements)};
if lengthParam = false then
 return arrow(listOfPairs, scaling = constrained, shape = arrow)
else
 return arrow(listOfPairs, length = lengthParam, scaling = constrained, shape = arrow)
end if
end proc

Essentially, we plot vectors v at location u. Since the vectors v in general vary wildly in magnitude, I would like to visualize the magnitude change by using color and/or transparency, preferably the latter. Is it possible to do so, and if so, how?

Is it possible to plot the equation 1=x^2+(y-1)^2 within Maple without rewriting it into a function or something like that? I can see that it is an easy to plot circle, but I just wonder whether it is possible.

Hi, 

i have this:

and i don't want the letter to be strings.

I want something like:

[z,d,i,p,s,...,y]

Thank you a lot

I am thinking about plotting a star chart of the whole sky. While it might not be difficult to get the (x,y)-coordinates of fairly good number of stars from star catalogues and to plot them accordingly in a 2d-plot, it might be more difficult to plot the stars according to their brightness. Also, I would like to plot the chart in cylindrical coordinate with astronomical coordinate axes for x and y. Moreover, only the brighter stars should be visible when the whole chart is shown, but when zooming in (while still having the axes visible) to have a closer look at some specific locations where I wand to make some square to show some fields, fainter stars should appear.

Is this possible in Maple ? I'd prefer to do a plot as I can really draw several squares where I want them which is not possible in my planetarium software.

Trying to solve the 1-dimensional heat equation with maple with constant boundary temperatures:

restart;

with(PDETools):
U := diff_table(u(x,t)):
pde := U[t]=U[x,x];
bc := u(0, t) =0, u(1, t) = 1, u(x,0)=x;
pdsolve([pde,bc]);

The solution of this equation is u(x,t)=x , but pdsolve(...) does not return anything at all! What is going wrong? Is it too hard PDE for maple? And if it is too hard, where can be found the types of equations, which are too hard and not too hard? Thank you.

how can you graph when piecewise function is not continuous?

example:

BB := piecewise (-1.57 < c and c < -1.56, h, -0.06< c and  c < -0.05, aa, -0.5< c and c < 0.04, bb);

so if i want to perform a function on elements from a list i do 

X:=[3, 5, 6, 7]

add(x*ln(x), x = X)

but what if i want to have two lists X and Y and perform the function x_i*ln(y_i) and add that up?

PS: keep in mind that im working with integers aswell as floats

thanks :)

how to graph in maple 

for example

-2 < x < -3, h

-1 < x < -2, b

why do I get the error Error, (in rtable/Sum) invalid arguments

In positive numbers, I get it ok

Hello

I am trying to do a do loop

X[1] := Statistics:-Sample(Binomial(200, .5), 1);

for i to 10 do X[i+1] := Statistics:-Sample(Binomial(200-X[i], .5), 1) end do;

And it shows error like this

Error, (in rtable/Sum) invalid arguments.

Did I do something wrong in my code?

I would like, for an arbitary point in R^2, to calculate the projection onto a curve, such as 1/x or 1/exp(x). Is there a "cheap" way to do so? If it makes it easier, I could mostly be interested in points strictly below the curve.

Dear all,

I'm trying to solve a system of two differential equations which depend on two parameters, which I call "K" and "polytrope" in my worksheet. The two differential equations are for the density rho(r) of the star and it's mass m(r). What I'm trying to do with dsovle is to change the values of K and polytrope once rho(r) drops below a certain value and letting dsolve continue the integration with the new values for the parameters. I was thinking about using events, but I'm not sure how to do it, since Maple doesn't seem to allow the use of "<" and ">" in the trigger of an action and I'm not sure how to let dsolve distinguish between the cases rho(r)<rho_crit and rho(r)>=rho_crit.

Thank you for your help.

Download TOV_with_2_EOS.mw

Change the coefficient of restitution to 0.80 and the ball takes off through the terrain. 

Why does it do that and what is the fix?

Hello,

I have modeled a slidercrank with absolute coordinates.

I have a slight problem when I want to mesure the Euler angles with Euler angle sensors.

In fact, the angle of the crank around the z axis (gamma1) always stay betwen -pi and pi. I would have like to see the value of my angle always increases with the evolution of the time since I impose a constant torque on the crank.

How can I do to observe the Euler angles of a body (parameterized with absolute coordinates) without any limit between [-pi,pi]?

Thank you for your help.

PS: I insert my model attached

slidercrank_3D_AC.msim