Hello,

I often need to copy equations from maple with mathml option so as to copy it in my equations in mathtype in microsoft word.

Unfortunetaly, for a bit long equation, i often have troubles et i can't copy it in my mathtype equations.

Have you some ideas so as the copy in mathtype equations runs everytime ?

Thanks a lot for your help

Hi MaplePrimers,

I've written a function that needs to be evaluated at a bunch of different points.  Evaluating it in a loop works. I'm wondering if there was a faster way to do this because I'm evaluating a list of 400k+ sets of points.  

Right now I'm using a loop with the following code.

for i from 1 to 500000 do 
     Results(i):= f(L[i][])[1];
od;

 

If I have a function f, which has 5 arguements, f(y, w, x,y,z).

I also have a list of those arguements:  L = [[1,2,3,2,3],[4,5,6,2,3],[7,8,9,2,3]]

What would be the fastest way to get a list of results?  Also, is there a way that I could preallocate memory for this list?

Ideally, I would like to get the output as as list.  ie:

Results:= [[f(1,2,3,2,3)], [f(4,5,6,2,3)], [f(7,8,9,2,3)]]

Hi all,

As we know that the differential of conjugate(z) is non-analytic.

But it seems that the diff gives the solution about the complex function conjugate(z)

What does it mean?

Hello guys, i have a system of equations ( dynamical system ) which i have its critical points but when i compute its critical points with maple i get different points , i dont know what is wrong . thank you for your time.

critical.mw

Is there any one to help me in calculating computational order of convergence of derivative free method. I am trying to calculate programs runs with immature error and does not print it.

Please is there any 1 to help.

Hello

Is it not possible to view page breaks in ypyr document in Maple? 

It will make it easier to at up at nice document??

Jakob

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:

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

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

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