Items tagged with maple maple Tagged Items Feed

Hello, 

is there a way I can use data (variables) from Maple environment in the Maplesim environment. 

I have a scirpt in maple that generates the robots joints angles and need to use them in the 3D robot built in maplesim. I know I can export/Import data, but this sounds redundant. Is there a way to simply use an input block as a source of the data in maplesim and have the variable name generated in maple used int. Similar to what Matlab/Simulink does.. 

 

 

thanks.

I'm attempting to plot the Nyquist plot for a complex system such that Maple cannot determine the frequency limits automatically, and suggested I use 'range' to specify the frequency limits. How do you do this? For example if you have a system G=1/(s+1), how do you plot the Nyquist plot in the frequency range 0 to 1rad/s?  Thanks

hi friends

i have a problem in maple with an error

  

Read'dsnumsort.map';
dsnumsort(Ns(0),[x,y,z]);
for i from 0 to 1000 do;
T:=i/25;
NsT:=Ns(T):
X[i]:=rhs(NsT[C1]); Vx[i]:=rhs(NsT[V1]);
Y[i]:=rhs(NsT[C2]); Vy[i]:=rhs(NsT[V2]);
Z[i]:=rhs(NsT[C3]); Vz[i]:=rhs(NsT[V3]);
KepVec[i]:=convert(crossprod([X[i],Y[i], Z[i]],[Vx[i],Vy[i], Vz[i]]),list);
KepAbs[i]:=norm(KepVec[i],2);
od:

pacecurve({[seq(X[i],Y[i],Z[i],i=0..1000)],[[-1/(2),0,0],[1/(2),0,0]]},labels=['x','y','z']);

but i see:

Error, invalid input: seq expects 2 arguments, but received 4

can you helpe me?

Thanks

Maplesim 7.0(1) installer looks for Maple 18 installation.  Will Maplesim work with Maple 2015?

If not, what is the expected date of making both products compatible?

Regards,

G

Here is a sample code:

 

restart;
N := 3;

for i to 2 do

   f[i] := unapply(a[0][i]+Sum(a[j][i]*t^j, j = 1 .. N), t) end

do;

sys : = {diff(x(t),t) = y(t), diff(y(t),t) = x(t) }

syssub := subs([x(t) = f[1](t), y(t) = f[2](t)],sys)

simplify(syssub)

   solve(simplify(syssub) union {simplify(f[1](0))=2, simplify(f[2](0))=1}),
     {a[0][1], a[0][2], a[1][1], a[1][2], a[2][1], a[2][2]}
);

 

Now for some reason solve returns nothing. This doesn't make sense to me since the system is determinable.

Here is a proof,

 ans: = dsolve(sys union {x(0) = 2, y(0) = 1},{x(t), y(t)}, type = series)

 

will actually return the result.

 

The reason I am doing this is because I am currently working on a bigger nonlinear DE system and it couldn't be done with just dsolve and since there isn't a package in Maple that makes a series representation and plug and chug, I have to write out the necessary steps.

 

Something tells me that Maple has a lot of trouble getting around the variable "t". But it still makes no sense why it can't even return a[0][1] = 1,a[0][2]=2

hi friends

i have a problem in maple with an error

> Us:=subs(G=1,m=1,L=1,U):
> D2r:=[diff(x(t),t,t),diff(y(t),t,t),diff(z(t),t,t)]:
> g:=subs(x=x(t),y=y(t),z=z(t),grad(Us,[x,y,z])):
> IniC:=x(0)=1,D(x)(0)=0,y(0)=0,D(y)(0)=1,z(0)=3/4,D(z)(0)=0:
> Ns:=dsolve({seq(D2r[i]=g[i],i=1..3),IniC},{x(t),y(t),z(t)},numeric);

after this i see :

Error, (in f) unable to store 'grad(U, [1., 0., .750000000000000000])[1]' when datatype=float[8] 

meaning it can only store 8 byte hardware floats.
Is this statement true? And how can i solve this problem? Or could/should i use a different data type?

Thanks for your time

Best regards

Could anyone please hepl me? I have the following system

e1 := exp(F(r)/phi_0)*L*A(r) = (1/2)*(2*(diff(A(r), r, r))*B(r)*A(r)*r*C(r)+2*B(r)*A(r)*(diff(A(r), r))*(diff(C(r), r))*r-(diff(A(r), r))^2*B(r)*r*C(r)-(diff(A(r), r))*(diff(B(r), r))*A(r)*r*C(r)+4*B(r)*A(r)*(diff(A(r), r))*C(r))/(B(r)^2*A(r)*r*C(r));
e2 := alpha*(diff(F(r), r, r))+(alpha^2+omega)*(diff(F(r), r))^2+(1/4)*(4*(diff(C(r), r, r))*B(r)*A(r)^2*C(r)*r+2*(diff(A(r), r, r))*A(r)*B(r)*r*C(r)^2-2*B(r)*A(r)^2*(diff(C(r), r))^2*r-(diff(A(r), r))^2*B(r)*r*C(r)^2-2*A(r)^2*C(r)*(diff(C(r), r))*(diff(B(r), r))*r-(diff(A(r), r))*(diff(B(r), r))*A(r)*r*C(r)^2+8*B(r)*A(r)^2*C(r)*(diff(C(r), r))-4*A(r)^2*C(r)^2*(diff(B(r), r)))/(r*A(r)^2*B(r)*C(r)^2)-(1/4)*(2*(diff(A(r), r, r))*B(r)*A(r)*r*C(r)+2*B(r)*A(r)*(diff(A(r), r))*(diff(C(r), r))*r-(diff(A(r), r))^2*B(r)*r*C(r)-(diff(A(r), r))*(diff(B(r), r))*A(r)*r*C(r)+4*B(r)*A(r)*(diff(A(r), r))*C(r))/(B(r)*A(r)^2*r*C(r)) = 0;
e3 := (1/4)*(-2*(diff(C(r), r, r))*B(r)*A(r)*r^2-B(r)*(diff(A(r), r))*(diff(C(r), r))*r^2+A(r)*(diff(C(r), r))*(diff(B(r), r))*r^2-8*B(r)*A(r)*(diff(C(r), r))*r-2*B(r)*(diff(A(r), r))*C(r)*r+2*A(r)*C(r)*(diff(B(r), r))*r+4*B(r)^2*A(r)-4*B(r)*A(r)*C(r))/(B(r)^2*A(r)) = -(1/4)*(2*(diff(A(r), r, r))*B(r)*A(r)*r*C(r)+2*B(r)*A(r)*(diff(A(r), r))*(diff(C(r), r))*r-(diff(A(r), r))^2*B(r)*r*C(r)-(diff(A(r), r))*(diff(B(r), r))*A(r)*r*C(r)+4*B(r)*A(r)*(diff(A(r), r))*C(r))*r/(B(r)^2*A(r)^2);
e4 := -(alpha^2+2*omega)*(diff(F(r), r))*(-(1/2)*(-(diff(A(r), r))*B(r)*r^4*C(r)^2-A(r)*(diff(B(r), r))*r^4*C(r)^2-4*A(r)*B(r)*r^3*C(r)^2-2*A(r)*B(r)*r^4*C(r)*(diff(C(r), r)))/(A(r)*B(r)*r^4*C(r)^2)-(diff(B(r), r))/B(r)+(diff(F(r), r, r))/(diff(F(r), r))+alpha*(diff(F(r), r)))/B(r) = -exp(F(r)/phi_0)*V_0*(alpha-1/phi_0);

phi_0 := -alpha/(2*alpha^2+2*omega); L := V_0*(1-(alpha-1/phi_0)*alpha/(3*alpha^2+2*omega)); V_0 := -lambda*exp(-fc/phi_0); fc := ln((4*alpha^2+2*omega)/(G_0*(3*alpha^2+2*omega)))/alpha; m := (2/(1+g))^(1/2); n := g*(2/(1+g))^(1/2); P := (G_0*(3*alpha^2+2*omega)/(4*alpha^2+2*omega))^(-2*alpha/(n-m)); eta := 1.4*G_0*Ms*(2/(1+g))^(-1/2)/c^2; g := 1-alpha^2/(2*alpha^2+omega);

omega := -10^5; alpha := 1; G_0 := 6.67*10^(-11); lambda := 10^(-52); c := 2.9*10^8; Ms := 1.9*10^30;
ri := evalf(1000*eta);

ics := A(2.109660445*10^6) = 1, (D(A))(2.109660445*10^6) = 2.370091128*10^(-15)*sqrt(2)*sqrt(99998)*sqrt(199997), B(2.109660445*10^6) = 1, C(2.109660445*10^6) = 1, (D(C))(2.109660445*10^6) = 4.740182256*10^(-15)*(1-(99999/19999300006)*sqrt(2)*sqrt(99998)*sqrt(199997))*(1-1.000017501*10^(-8)*sqrt(2)*sqrt(99998)*sqrt(199997))^(-(99999/19999300006)*sqrt(2)*sqrt(99998)*sqrt(199997))*sqrt(2)*sqrt(99998)*sqrt(199997), f(2.109660445*10^6) = 23.43081116, (D(f))(2.109660445*10^6) = 4.749681180*10^(-15):

eta:=2109.660445: sys:=e1,e2,e3,e4; vars:=[A(r),B(r),C(r),F(r)];

dsn3 := dsolve([sys, ics], numeric, vars, range = 3*eta .. 50*eta);

Results in

Warning, cannot evaluate the solution past the initial point, problem may be complex, initially singular or improperly set up

Setting f(r)=Const,V_0=0 which is a physically relevant case, results in

Error, (in simplify/normal) numeric exception: division by zero

I suugest the problem is that the equation contain sqared derivatives, hence there are several solution branches corresponding to different signs of square root. Maple chooses the singular branch. How can I force it to choose another branch or calculete all of them?

Thanks in advance..





Hi all,

 

It's been a while since I have used Maple. To be honest I haven't used it for over six years.

 

I am trying to solve simple differential equations, however I have many issues.

 

I am trying to simulate what author of this paper did 06421188.pdf

 

My file looks like this (Pendulum.mw)

 

Can someone help me to simulate this system? I simply can't remember how to do it.

 

Cheers,

Bart

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})

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

Hi,

 

The default way a prompt shows on my maple is

 

[> blah blah blah

 

but i see in many worksheets including maple help pages that it shows as

 

> blah blah blah

 

My Question is how do i remove that left square bracket "[" from appearing ?

There was a command for torusknot(x,y,z),r,R.  Is it included in Maple 17? The index and help to do not list it.

Is there a newer version of a torus knot?

I am having issues opening my final year project, it was working an hour ago and now it will not open. When I try to open the file it brings up a box title TEXT FORMAT CHOICE with the options: MAPLE TEXT, PLAIN TEXT, MAPLE INPUT and CANCEL.

This is all my work and I need it to complete my year.

Any help with how I can rectify this would be extremely helpful.

Many thanks

Hi,

I'm, trying to write a nested loop in maple, but I keep getting the Unterminated Loop error. I'm sure the solution is quite simple, but I haven't been able to find it. Any suggestions would be much appreciated. The loop looks like this:

 

i:=0:
di:=25:
n:=1:  

while (n<=nLimit) do
sol1:=fsolve(eq1, y=i..i+di)
     if type(sol1, numeric) = true then
     lambda(n) :=sol1
     i:=i+di
     n:=n+1
     else  i:=i+di
     end if:
end do:

I have the following characteristic equation by use of maple. How do I find a condition on x, that will return real eigenvalues and complex eigenvalues?

 

 

 

1 2 3 4 5 6 7 Last Page 3 of 270