Hi! I'm trying to make a point plot and they must have different colors, like the CIE diagram. I wish they varied from red to blue in 31 steps, as shown in the picture (they are all black :\). Does anyone have any idea how to do this? The red line is for something else.

This the code for the dots:

Hi,

Can anyone help me to solve this system of equations in Maple?

solve({
(-(1*Rfd*(1-(Ladssec/Lfd)+((((1.007033 +1*(Laqssec+Ll))*Ladssec)/(((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll)))*Lfd)) *Ladssec)))/(Lfd))=-0.0069
,(-(1*Rfd*(-(Ladssec/Lfd)+((((1.007033 +1*(Laqssec+Ll))*Ladssec)/(((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll)))*L1d)) *Ladssec)))/(Lfd))=0.002689
,(-(1*Rfd*(-((Ra+0.04527646)/((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll))))*(Laqssec/L1q))*Ladssec)/(Lfd))=0.00002647
,(-(1*Rfd*(-((Ra+0.04527646)/((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll))))*(Laqssec/L2q)) *Ladssec)/(Lfd))=-0.00001362

,(-(1*R1d*(-(Ladssec/Lfd)+((((1.007033 +1*(Laqssec+Ll))*Ladssec)/(((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll)))*Lfd))*Ladssec)))/(L1d))=0.1052
,(-(1*R1d*(1-(Ladssec/Lfd)+((((1.007033 +1*(Laqssec+Ll))*Ladssec)/(((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll)))*L1d)) *Ladssec)))/(L1d))=-0.2585
,(-(1*R1d*(-((Ra+0.04527646)/((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll))))*(Laqssec/L1q))*Ladssec)/(L1d))=0.0009122
,(-(1*R1d*(-((Ra+0.04527646)/((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll))))*(Laqssec/L2q)) *Ladssec)/(L1d))=-0.0005093

,(-(1*R1q*(((Ra+0.04527646)/((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll))))*(Ladssec/Lfd))*Laqssec)/(L1q))=-0.000292
,(-(1*R1q*(-((Ra+0.04527646)/((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll))))*(Ladssec/L1d))*Laqssec)/(L1q))=0.0008507
,(-(1*R1q*(1-(Laqssec/L1q)+((((1.007033 +1*(Ladssec+Ll))*Laqssec)/(((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll)))*L1q))*Laqssec)))/(L1q))=-0.04423
,(-(1*R1q*(-(Laqssec/L2q)+((((1.007033 +1*(Ladssec+Ll))*Laqssec)/(((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll)))*L2q))*Laqssec)))/(L1q))=0.0003831

,(-(1*R2q*(((Ra+0.04527646)/((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll))))*(Ladssec/Lfd))*Laqssec)/(L2q))=-0.001785
,(-(1*R2q*(-((Ra+0.04527646)/((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll))))*(Ladssec/L1d))*Laqssec)/(L2q))=-0.03134
,(-(1*R2q*(-(Laqssec/L1q)+((((1.007033 +1*(Ladssec+Ll))*Laqssec)/(((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll)))*L1q))*Laqssec)))/(L2q))=4.239
,(-(1*R2q*(1-(Laqssec/L2q)+((((1.007033 +1*(Ladssec+Ll))*Laqssec)/(((Ra+0.04527646)^2+(1.007033 +1*(Ladssec+Ll))*(1.007033 +1*(Laqssec+Ll)))*L2q))*Laqssec)))/(L2q))=-4.868

}, {L1d, L1q, Ladssec, Laqssec, Lfd, R1d, R1q, Rfd, Ra, Ll},useassumptions)

assuming Lfd ~= 0.1632, R1d~=0.0269, L1d~=0.0750, L1q~=-0.1350, Ladssec~=0.0500, Laqssec~= -0.1500, Rfd~=0.0016, R1q~= 0.0215, Ra~=0.005, Ll~=0.15

which theory can explain the interaction of composition of group for combination of composition of group

if succeed to search a list of groups, what is the next step research them?

Using maple 16. When I run a worksheet it adds in extra rows of the red >s. The attached example was run 3 times. How do I stop this happening?

Download worksheet_issue.mw

Hi, i have problem with this code: 

Spacer := proc (N, p, q, sa, sb, M) local A, B, r, PA, PC, PB, stan, krok, droga; A := -N; B := N; r := 1-p-q; PA := empirical[0., sa, 1-sa]; PC := empirical[q, r, p]; PB := empirical[1-sb, sb, 0.]; stan := 0; droga := stan; to M do if A < stan and stan < B then krok := random[PC](1)-2 elif stan = A then krok := random[PA](1)-2 elif stan = B then krok := random[PB](1)-2 end if; stan := stan+krok; droga := droga, stan end do; {}*plots[pointplot]([seq([i, droga[i+1]], i = 0 .. nops([droga])-1)], axes = BOXED) end proc

---------------------------------------------------------------------------------------------------------------------------------

i have error when i start procedure like: 

error, (in Spacer) cannot determine if this expression is true or false: 0 < random[empirical[.5, 0.5e-1, .45]](1)+8 and random[empirical[.5, 0.5e-1, .45]](1) < 12

My run command: Spacer(3,0.5,0.5,0.8,0.8,20);

Please help me.

Hello Hello everybody 
   I have to solve the following differential equation numerically 

I don't know why it says : Error, (in f) unable to store '[0.]/(0.17571268341557e16+[-0.25659510610770e15])' when datatype=float[8]

Help pleaase!

thank you !!!

Download systéme_complet.mw

Hi,

I have a file with arbitrary points which I want to connect with a surface. I have seen similar problems on this site but I can't fix my problem with the answers that are given, mainly because my points are just plain vectors. Can anybody help me? 
These are my points/vectors:

v1 := `<,>`(15, 6, 14.89);
v2 := `<,>`(15, 8, 14.72);
v3 := `<,>`(15, 10, 14.79);
v4 := `<,>`(15, 12, 14.95);
v5 := `<,>`(15, 14, 15.17);
v6 := `<,>`(15, 16, 15.43);
w1 := `<,>`(15/2, 4, 9.78);
w2 := `<,>`(15/2, 6, 9.95);
w3 := `<,>`(15/2, 8, 9.85); 
x1 := `<,>`(15*(1/3), 4, 8.71);
x2 := `<,>`(15*(1/3), 6, 8.44);
y1 := `<,>`(15/4, 4, 8.31); 

I can only plot the different points with 

with(plots);
pointplot3d({v1, v2, v3, v4, v5, v6, w1, w2, w3, x1, x2, y1}, scaling = CONSTRAINED, symbolsize = 20, color = red, axes = normal, labels = ["L/h", "n", "kg/m"]);

I hope someone can help me!

Hello,

I am trying to solve the equation (non-linear) with one variable t2 (please see bellow):

(35595.29412*(52040.0-2400.0*t2))*(11-t2)*exp(ln(-360.0*t2+29600.0)-(0.9803921570e-5*(-360.0*t2+59200.0))*t2)*exp(ln(52040.0-2400.0*t2)-8.804313725+.5176470590*t2)*(360000*(exp(ln(-360.0*t2+29600.0)-(0.9803921570e-5*(-360.0*t2+59200.0))*t2))^2+(620*(-360.0*t2+29600.0))*(1200*exp(ln(-360.0*t2+29600.0)-(0.9803921570e-5*(-360.0*t2+59200.0))*t2)-2.232000*10^5*t2+1.83520000*10^7))*exp(ln(-360.0*t2+29600.0)-(0.9803921570e-5*(-2760.0*t2+59200.0))*t2)/(360000*exp(ln(52040.0-2400.0*t2)-8.804313725+.5176470590*t2)*(exp(ln(-360.0*t2+29600.0)-(0.9803921570e-5*(-360.0*t2+59200.0))*t2))^2+(372000*((52040.0-2400.0*t2)*exp(ln(-360.0*t2+29600.0)-(0.9803921570e-5*(-2760.0*t2+59200.0))*t2)+(-360.0*t2+29600.0)*exp(ln(52040.0-2400.0*t2)-8.804313725+.5176470590*t2)))*exp(ln(-360.0*t2+29600.0)-(0.9803921570e-5*(-360.0*t2+59200.0))*t2)+(384400*(52040.0-2400.0*t2))*(-360.0*t2+29600.0)*exp(ln(-360.0*t2+29600.0)-(0.9803921570e-5*(-2760.0*t2+59200.0))*t2))^2 = 1/150

During evaluation of the solve command I received a warning that solutions may have been lost. How can I overcome this problem? Also, I need that t2>0.

Thanks in advance,

Dmitry 

Hi,

I have a linear system to solve.

mm:=proc(a,x,h,i)
local A,Z1,Z2,Z,F,result;  # to declare the local variable
A:=array(1..2,1..2,[[1,1],[a,a+h]]);
Z1:=evalf(int(1/(abs(y-x)+.000000001),y=a..a+h));
Z2:=evalf(int(y/(abs(y-x)+.000000001),y=a..a+h));
Z:=array([Z1,Z2]);
F:=evalm(inverse(A)&*Z);
result:=F[i]
end:

My questions: 

1) My exact Z1 is Z1:=evalf(int(1/(abs(y-x)),y=a..a+h)); but I ask if can I put

Z1:=evalf(int(1/(abs(y-x)+.000000001),y=a..a+h));

the same for Z2.

2) Can I writte in a simple form the vector Z.  Because, later, il have a second system contains Z1,Z2, Z3, Z4,Z5.  The difference between Z1 and Z2 is the variable "y" added in the integral of Z2.

Many thinks.

I want to create a NxN circulant matrix of a system with fixed boundaries. This is an idea of what I'm trying to do (but it doesn't work):

> N := 8;

> A := Matrix(N, shape = Circulant[[2, VectorCalculus[`-`](1), ZeroVector[row](VectorCalculus[`+`](N, VectorCalculus[`-`](3))), VectorCalculus[`-`](1)]]);

>A[1,N]:=0;

>A[N,1]:=0;

In other words, I want to make a circulant matrix of the elements 2, -1, 0, 0, ..., -1 (with N-3 zeros). Then I want to set the top right and bottom left elements to zero (that's the fixed boundary portion). The first command almost works for N=4 -- except with extra square brackets in the matrix -- but doesn't work for N=8 or 16 (N is always a power of 2). The second and third commands seem to be bad syntax (I haven't really tried to figure them out yet).

Google translate is HTTPS.  Https protocol is not supported in the HTTP package.  I can't seem to get it to work using the Sockets package.

Any ideas?

i will convert matlab code to maple please help me

clc
clc
clear all
E1=181*10^9;
E2=10.30*10^9;
G12=7.17*10^9;
nu_12=0.28;
nu_21=(nu_12*E2)/E1;
fprintf('Please enter a numerical value for Total Height.')
htotal=input('htotal=');
fprintf('Number of entries for ply orientations should equal the number of entries for aspect ratios .','Please enter values separated by spaces in [].')
O(1,:)=input('ORIENTATION=');
fprintf('Please enter values separated by spaces in [] for the Independent Ply Aspect Ratios(=Thickness of each ply/Total Height).')
AR=input('ASPECT-RATIO=');
[ARrows,ARcols] = size(AR(1,:));
[orows,ocols] = size(O(1,:));
if ocols~=ARcols
fprintf('Check Number of Independent Ply Orientation and Aspect Ratio.')
end
plycount = ocols;
k=0;
sumAR = sum(AR);
H = htotal*(AR/sumAR)
% Define matrix with ply orientations as row one and corresponding height
% as row two for beter display
OHmatrix = [O;H]
Z=zeros(1,plycount+1);
hindex = -(htotal)/2 ;
for i=1:plycount+1
if i==1
Z(:,i) = hindex;
else
Z(:,i) = Z(:,i-1)+ H(i-1);
end
end
s11=1/E1;
s12=-nu_12/E1;
s16=0;
s21=s12;
s22=1/E2;
s26=0;
s61=s16;
s62=s26;
s66=1/G12;
S=[s11 s12 s16;s21 s22 s26;s61 s62 s66];
Q11 = E1/(1- nu_12*nu_21);
Q22 = E2/(1- nu_12*nu_21);
Q12 = (nu_21*E1)/(1- nu_12*nu_21);
Q66 = G12;
% note that Q21 = Q12
Qp = [ Q11, Q12 , 0 ;Q12, Q22 , 0 ;0 , 0 , Q66 ];
% Intializing ABD as a 3x3 zero matrix
A = zeros(3,3) ;
B = zeros(3,3) ;
D = zeros(3,3) ;
% Calculating the A, B, D Matrices for each ply
Qp(3,3) = Qp(3,3)*2 ;
for l = 1 : plycount
thetar = (O(l)/180)*pi;
% define "m" and "n" as
m = cos(thetar);
n = sin(thetar);
% 2D Transformation matrix T :
T = [m^2 , n^2 , 2*m*n ;n^2 , m^2 , -2*m*n ;-m*n , m*n , m^2 - n^2 ];
% Correcting for engineering strain-true strain
Q = T\Qp*T ;
Q(:,3) = Q(:,3)*.5 ;
l
Q
for i=1:3
for j=1:3
A(i,j)= A(i,j) + (Q(i,j))*( Z(l+1) - Z(l) );
B(i,j)= B(i,j) + 0.5* ((Q(i,j))*( (Z(l+1))^2 - (Z(l))^2 ));
D(i,j)= D(i,j) + (1/3)*((Q(i,j))*( (Z(l+1))^3 - (Z(l))^3 ));
end
end
for i=1:3
for j=1:3
if abs(A(i,j))< 1.0e-4
A(i,j) = 0;
end
if abs(B(i,j))< 1.0e-4
B(i,j) = 0;
end
if abs(D(i,j))< 1.0e-4
D(i,j) = 0;
end
end
end
end
%Finding AB-BD matrix and its inverse
ABDmatrix = [A B ; B D];
A = ABDmatrix(1:3, 1:3)
B = ABDmatrix(1:3, 4:6)
D = ABDmatrix(4:6, 4:6)

Hi;

I've written a Maple's procedure but when I execute it some warning messages have been displayed as follows:

read "My_Proc.mpl":
Warning, `L` is implicitly declared local to procedure `My_Proc`
Warning, `P` is implicitly declared local to procedure `My_Proc`
Warning, `V` is implicitly declared local to procedure `My_Proc`

So, how can I get rid of such warnings?
Thanks in advance.

I wish to  plot 2D animate for the soltion of this equation here is the code

restart;
with(PDEtools):
with(ArrayTools):
with(plots):

f:=u->sech(u):
g:=v->sech(v):
h:=1/10:
N:=20:
M:=20:
V:=x->x^2:
psi:=Array(0..N/h+1,0..M/h+1):
for i from 0 to N/h do
psi[i,0]:=evalf(f(i*h)):
od:
for j from 0 to M/h do
psi[0,j]:=evalf(g(j*h)):
od:
for i from 1 to N/h do
for j from 1 to M/h do
psi[i,j]:=-psi[i-1, j-1]+(1-(1/8)*h^2*V((1/2)*h*(j-i-1)))*psi[i, j-1]+(1-(1/8)*h^2*V((1/2)*h*(j- i+1)))*psi[i-1,j]:
od:
od:
ls:=[seq([seq([i*h,j*h,psi[i,j]],i=0..N/h)],j=0..M/h)]:

surfdata((ls),axes=boxed,labels=[`u`,`v`,`psi(u,v)`],shading=zhue,style=patchcontour);

Group Theory,

Given a Group G and a normal Subgroup H.

I can get " map(Representative,RightCosets(H,G)) ".

Is ist possible to get the factor group  G/N ?

Best regards

Kurt Ewald