Items tagged with maple maple Tagged Items Feed

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?

 

 

 

Hi, I am completely new to Maple, and I need to use it to optimize my equations in order to make my PLC codes more compressed. I am calculating forward kinematics with the Denavit-Hartenberg method and as such I get long expressions. After a lot of google'ing and frustration, I thought I'd ask here in the hope that one of you might be able to assist me.

I have the following equations;

X := L10*cos(q5) - L16*(sin(q10)*(sin(q5)*sin(q8) - cos(q8)*(cos(q5)*cos(q6)*cos(q7) - cos(q5)*sin(q6)*sin(q7))) - cos(q10)*(sin(q9)*(cos(q8)*sin(q5) + sin(q8)*(cos(q5)*cos(q6)*cos(q7) - cos(q5)*sin(q6)*sin(q7))) + cos(q9)*(cos(q5)*cos(q6)*sin(q7) + cos(q5)*cos(q7)*sin(q6)))) - d2*(cos(q10)*(sin(q5)*sin(q8) - cos(q8)*(cos(q5)*cos(q6)*cos(q7) - cos(q5)*sin(q6)*sin(q7))) + sin(q10)*(sin(q9)*(cos(q8)*sin(q5) + sin(q8)*(cos(q5)*cos(q6)*cos(q7) - cos(q5)*sin(q6)*sin(q7))) + cos(q9)*(cos(q5)*cos(q6)*sin(q7) + cos(q5)*cos(q7)*sin(q6)))) + L15*(sin(q9)*(cos(q8)*sin(q5) + sin(q8)*(cos(q5)*cos(q6)*cos(q7) - cos(q5)*sin(q6)*sin(q7))) + cos(q9)*(cos(q5)*cos(q6)*sin(q7) + cos(q5)*cos(q7)*sin(q6))) - L11*cos(q5)*sin(q6) + d1*cos(q5)*cos(q6) - L13*sin(q5)*sin(q8) + L14*cos(q9)*(cos(q8)*sin(q5) + sin(q8)*(cos(q5)*cos(q6)*cos(q7) - cos(q5)*sin(q6)*sin(q7))) + L13*cos(q8)*(cos(q5)*cos(q6)*cos(q7) - cos(q5)*sin(q6)*sin(q7)) - L14*sin(q9)*(cos(q5)*cos(q6)*sin(q7) + cos(q5)*cos(q7)*sin(q6)) + L12*cos(q5)*cos(q6)*cos(q7) - L12*cos(q5)*sin(q6)*sin(q7);

Y := L10*sin(q5) - L9 + L16*(sin(q10)*(cos(q5)*sin(q8) - cos(q8)*(sin(q5)*sin(q6)*sin(q7) - cos(q6)*cos(q7)*sin(q5))) - cos(q10)*(sin(q9)*(cos(q5)*cos(q8) + sin(q8)*(sin(q5)*sin(q6)*sin(q7) - cos(q6)*cos(q7)*sin(q5))) - cos(q9)*(cos(q6)*sin(q5)*sin(q7) + cos(q7)*sin(q5)*sin(q6)))) + d2*(cos(q10)*(cos(q5)*sin(q8) - cos(q8)*(sin(q5)*sin(q6)*sin(q7) - cos(q6)*cos(q7)*sin(q5))) + sin(q10)*(sin(q9)*(cos(q5)*cos(q8) + sin(q8)*(sin(q5)*sin(q6)*sin(q7) - cos(q6)*cos(q7)*sin(q5))) - cos(q9)*(cos(q6)*sin(q5)*sin(q7) + cos(q7)*sin(q5)*sin(q6)))) - L15*(sin(q9)*(cos(q5)*cos(q8) + sin(q8)*(sin(q5)*sin(q6)*sin(q7) - cos(q6)*cos(q7)*sin(q5))) - cos(q9)*(cos(q6)*sin(q5)*sin(q7) + cos(q7)*sin(q5)*sin(q6))) + L13*cos(q5)*sin(q8) - L11*sin(q5)*sin(q6) + d1*cos(q6)*sin(q5) - L14*cos(q9)*(cos(q5)*cos(q8) + sin(q8)*(sin(q5)*sin(q6)*sin(q7) - cos(q6)*cos(q7)*sin(q5))) - L13*cos(q8)*(sin(q5)*sin(q6)*sin(q7) - cos(q6)*cos(q7)*sin(q5)) - L14*sin(q9)*(cos(q6)*sin(q5)*sin(q7) + cos(q7)*sin(q5)*sin(q6)) + L12*cos(q6)*cos(q7)*sin(q5) - L12*sin(q5)*sin(q6)*sin(q7);

Z := L15*(cos(q9)*(cos(q6)*cos(q7) - sin(q6)*sin(q7)) - sin(q8)*sin(q9)*(cos(q6)*sin(q7) + cos(q7)*sin(q6))) - L11*cos(q6) - L8 - d1*sin(q6) + L16*(cos(q10)*(cos(q9)*(cos(q6)*cos(q7) - sin(q6)*sin(q7)) - sin(q8)*sin(q9)*(cos(q6)*sin(q7) + cos(q7)*sin(q6))) - cos(q8)*sin(q10)*(cos(q6)*sin(q7) + cos(q7)*sin(q6))) - d2*(sin(q10)*(cos(q9)*(cos(q6)*cos(q7) - sin(q6)*sin(q7)) - sin(q8)*sin(q9)*(cos(q6)*sin(q7) + cos(q7)*sin(q6))) + cos(q8)*cos(q10)*(cos(q6)*sin(q7) + cos(q7)*sin(q6))) - L13*cos(q8)*(cos(q6)*sin(q7) + cos(q7)*sin(q6)) - L14*sin(q9)*(cos(q6)*cos(q7) - sin(q6)*sin(q7)) - L12*cos(q6)*sin(q7) - L12*cos(q7)*sin(q6) - L14*cos(q9)*sin(q8)*(cos(q6)*sin(q7) + cos(q7)*sin(q6));

 

I need to optimize these equations, but still keep them separate. I would like to use mutual expressions for the calculations within, but still as I said keep the outputs of X, Y and Z separate.

This is MATLAB code.

 

Thanks in advance for any help.

Hi,

 

I am new to Maple and have a problem when solving three equations with three variables. But when  I plug in into solve function then it gives no answer.

eqn1 := 24900 = A*exp(-X*1.293995859*10^22)+A*exp(-Y*1.293995859*10^22)+5852.27;

eqn2 := 6000 = A*exp(-X*1.293995859*10^22)+2422.929937;

eqn3 := 19100 = A*exp(-Y*1.293995859*10^22)+8275.199937;

Variables are [A,X,Y]

Hi!

I want to select 'Export using shapes for greater fidelity' under option, Export. But I can't find the tab ' Export to PDF-Format'. I'm using Maple 18 on a Mac-computer. 

Can anyone help?

Thanks

 

Esben 

For example, given a 3d point p(x,y,z), with (x,y,z) as its coordinates. Then it is transformed by rotation and translation, as 

p'=R(p)*p+t(p), where R(p) is a 3x3 rotation matrix that is a matrix of functions of p, and t(p) is a 3x1 vector function of p. 

My question is how to derive dp'/d(as a 3x3 matrix) using maple? 

To make it clear,I want to do it in a way that dp'/dp = ∂p'/∂p + ∂p'/∂R*∂R/∂p +  ∂p'/∂t*∂t/∂p

And I'd like to know each intermediate quantity, such as p'/∂R, R/∂p.


Anyone can help?

Thanks a lot. 

Using Maple 18, I solved for minimum and maximum price. Instead of using fsolve I wanna use procedure programming structure in order to get the same results. How can I do it?

min_sol := fsolve([bc_cond, slope_cond, x[G, 1] = w[aggr, 1]], {p = 0 .. 1, x[G, 1] = 0 .. w[aggr, 1], x[G, 2] = 0 .. w[aggr, 2]}); p_min := subs(min_sol, p); max_sol := fsolve([bc_cond, slope_cond, x[G, 2] = w[aggr, 2]], {p = 0 .. 1, x[G, 1] = 0 .. w[aggr, 1], x[G, 2] = 0 .. w[aggr, 2]}); p_max := subs(max_sol, p);
{p = 0.3857139820, x[G, 1] = 127.8000000, x[G, 2] = 38.99045418}
0.3857139820
{p = 0.8841007104, x[G, 1] = 44.30160890, x[G, 2] = 164.2000000}
0.8841007104

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