Hello,

i need to represent a circle using the command [geometry] circle because i have a problem and i have to calculate an 

intersection between two circles whit command intersection() which is ruled by Package geometry.

For example, with the first circle  i did like the help page says

 But it's Error:

centre of circle(3300;2200)

radius=110

I don't understand the eror. Can you help me?

Thanks in advance.

Hello!

I have a question about densityplot.

I want to get this colormap:

As I understand, I should use plots:-densityplot(...) comand and "HSV" colorscheme, where H = [240..0], S=1, V=1.

But if i use zcoloring with this options, I will get wrong result.

My code:

restart;

a := 0; b := 1; #first and final point

plots:-densityplot(z, dummy = 0 .. 1, z = a .. b, grid = [2, 10], size = [90, 100], colorscheme = ["zcoloring", [proc (z) options operator, arrow; -240*(z-a)/(b-a)+240 end proc, proc (z) options operator, arrow; 1 end proc, proc (z) options operator, arrow; 1 end proc], colorspace = "HSV"], style = surface, axes = frame, labels = ["", ""], axis[1] = [tickmarks = []], title = "Test", titlefont = ["Calibri", "Bold", 16], axesfont = ["Arial", 14], size = [100, 500]); #densityplot with zcoloring option.

My result (zero must be equal to "Blue"):

I can't find mistake....

I am trying to model a disease. The equation is as follows:

S*X - f(X,S,Sp) = 100

I have data for S, I have data for X and I have data for f(X,S,Sp) however I want to find an equation for f(X,S,Sp) that has the best fit with the data because I need to use it later on in my calculation. If anyone is intrested S is the sensitivity of the blood test, Sp is the specificity of the blood test. This means that X*(1-Sp) is the number of false positives.

I currently dont have access to Maple hence I am doing all my modeling in excel 2016. So I am severly limited because excel is useless at algebra.

Maple.xlsx

I have a list
restart;
L := [10*x+y-6*z+66 = 0, 26*x-10*y-3*z+96 = 0, 2*x-y+6 = 0, 2*x-y+6 = 0, 4*x-5*y+3*z-6 = 0, 30*x-22*y+7*z+48 = 0, 2*x-2*y+z = 0, 26*x-19*y+6*z+42 = 0, 22*x-17*y+6*z+30 = 0, 14*x-30*y+23*z-96 = 0, 6*x-10*y+7*z-24 = 0, 2*x-7*y+6*z-30 = 0, 4*x-y-z+10 = 0, 42*x-2*y-19*z+20 = 0, 38*x-2*y-17*z+20 = 0, 30*x-y-14*z+10 = 0, 18*x+2*y-11*z-20 = 0, 6*x+13*y-16*z-130 = 0, 18*x+7*y-16*z-70 = 0, 6*x+2*y-5*z-20 = 0, 14*x+15*y-22*z-150 = 0, 18*x+13*y-22*z-130 = 0, 6*x+4*y-7*z-40 = 0, 26*x-3*y-10*z+110 = 0, 10*x-6*y+z+52 = 0, 2*x-z+8 = 0, 26*x+6*y-19*z+92 = 0, 22*x+6*y-17*z+76 = 0, 2*x+y-2*z+6 = 0, 14*x-18*y+11*z+92 = 0, 2*x-6*y+5*z+20 = 0, 2*x+6*y-7*z-4 = 0, 10*x-3*y-2*z+22 = 0, 26*x-7*y-6*z+66 = 0, 2*x-y = 0, y-z+11 = 0, 4*x-9*y+7*z-77 = 0, 6*x-5*y+2*z-22 = 0, 24*x-13*y+z-11 = 0, 2*x-y = 0, 2*x-y = 0, y-z+11 = 0, 4*x-5*y+3*z-33 = 0, 6*x-4*y+z-11 = 0, 4*x-3*y+z-11 = 0, 18*x-19*y+10*z-110 = 0, 14*x-13*y+6*z-66 = 0, 2*x-2*y+z-11 = 0, 26*x-6*y-7*z+68 = 0, 10*x-2*y-3*z+24 = 0, 2*x-z+2 = 0, 6*x+y-4*z-1 = 0, 4*x+y-3*z-3 = 0, y-z-7 = 0, 4*x+3*y-5*z-17 = 0, 18*x+10*y-19*z-52 = 0, 8*x+17*y-21*z-111 = 0, 4*x+5*y-7*z-31 = 0, 4*x-3*y+z-2 = 0, 10*x-6*y+z+4 = 0, 14*x-10*y+3*z-4 = 0, 26*x-15*y+2*z+14 = 0, 22*x-13*y+2*z+10 = 0, 18*x-26*y+17*z-84 = 0, 2*x-6*y+5*z-28 = 0, 38*x-6*y-13*z-16 = 0, 6*x-y-2*z-2 = 0, 18*x-2*y-7*z-16 = 0, 6*x+y-4*z-22 = 0, 6*x+y-4*z-22 = 0, 14*x+11*y-18*z-138 = 0, 2*x+y-2*z-14 = 0, 2*x+y-2*z-14 = 0, 22*x+2*y-13*z+40 = 0, 14*x+3*y-10*z+22 = 0, 2*x-2*y+z+8 = 0, 2*x-5*y+4*z+14 = 0, 2*x+5*y-6*z-6 = 0, 2*x-z+2 = 0, 6*x+y-4*z+8 = 0, 22*x-12*y+z-2 = 0, 8*x-9*y+5*z-10 = 0, y-z+2 = 0, y-z+11 = 0, 2*x-y-9 = 0, 8*x-3*y-z-25 = 0, 34*x-18*y+z-164 = 0, 18*x-10*y+z-92 = 0, 14*x-10*y+3*z-96 = 0, 28*x-17*y+3*z-159 = 0, 8*x-5*y+z-47 = 0, 2*x+3*y-4*z+29 = 0, 4*x-3*y+z-23 = 0, 32*x-15*y-z-103 = 0, 20*x-9*y-z-61 = 0, 10*x-6*y+z-44 = 0, 14*x-10*y+3*z-76 = 0, 2*x-2*y+z-16 = 0, 4*x-3*y+z-23 = 0, 12*x-19*y+13*z-89 = 0, 16*x-19*y+11*z-87 = 0, 20*x-3*y-7*z-73 = 0, 14*x-3*y-4*z-43 = 0, 4*x+3*y-5*z-47 = 0, 16*x-3*y-5*z-53 = 0, 22*x-2*y-9*z-92 = 0, 12*x+5*y-11*z-105 = 0, 16*x+3*y-11*z-107 = 0, 12*x+y-7*z-13 = 0, 6*x+y-4*z-7 = 0, 10*x-14*y+9*z+4 = 0, y-z-7 = 0, 2*x-z-7 = 0, 2*x+9*y-10*z-70 = 0, 18*x+y-10*z-70 = 0, 12*x+11*y-17*z-119 = 0, 2*x+y-2*z-14 = 0, 2*x-4*y+3*z+15 = 0, 4*x+y-3*z-15 = 0, 2*x-11*y+10*z+50 = 0, 14*x+3*y-10*z-50 = 0, 8*x+13*y-17*z-85 = 0, 4*x+13*y-15*z-75 = 0, 2*x+5*y-6*z-30 = 0, 12*x-7*y+z-101 = 0, 2*x-y-15 = 0, 6*x-4*y+z-56 = 0, 6*x-7*y+4*z-47 = 0, 8*x-7*y+3*z-51 = 0, 2*x-4*y+3*z-24 = 0, 6*x+y-4*z-46 = 0, 2*x+3*y-4*z-10 = 0, 2*x-z-16 = 0, y-z+2 = 0, 6*x-2*y-z-52 = 0, 10*x-2*y-3*z-84 = 0, 10*x+y-6*z-78 = 0, 2*x+y-2*z-14 = 0]

Now I want to write the equation of the planes passing two points A and B and perpendicular to all the equations of the above list. I can't to write a proc. With each plane of List, I tried
with(geom3d);
point(A, 1, 5, 3): 
point(B, -2, -4, 3):
with(geom3d): 
Equation(plane(P, L[5], [x, y, z])): 
line(d, [A, B], t):
 n1 := ParallelVector(d):
 n2 := NormalVector(P):
 with(linalg):
 n3 := crossprod(n1, n2):
 k := sort(primpart(lhs(Equation(plane(R, [A, n3], [x, y, z])))));
k*signum(lcoeff(k)) = 0

I got true answer.

Hello,

How can I rotate vertically label of vertical axis ? Can you help me ?

Thank you.

I use maple 2016 on linux with X.  When minimized, the maple icon is very small.  Also, the help icon is very small.

Are larger icons available?

The maple icon is ./bin/Maple2016.png.  I can not find the help icon.  Any ideas?

I want to solve for the roots of a polynomial, such as a x^2+b x + c = 0, for which the output is only the positive root. All coefficients/variables in the polynomial are positive. 

Recently, someone posted an answer to a question where at some point they performed this task and their solution was really slick. But I can't find it. The answer used either solve, or eval or something like that. (Yes, I did perform a search via the MaplePrimes search before asking this question.) 

I am trying to calculate the integral

where

Maple cannot calculate the integral. I tried to expand theta in the series form and substitute in the integral, still cannot calculate it.

any suggestion to tackle this problem whould be helpful.

Thank you

I have a PDE with boundary conditions, from a NASA paper.  I always seem to have problems expressing the bc.

PDE := diff(u(x, t), t) = (1/4)*exp(2)*exp(-u(x, t))*(diff(diff(u(x, t), x), x))/(x^2+2);

The initial/boundary conditions are

@t=0, u(x, t) = 2-2*ln(-x^2+2)

@x=0, diff(u(x,t),x)=0  ## this is the bc I have problem expressing

@x=1, u(x,t) = 2+ln(1+t)

The exact solution given in the paper:

2 + ln(1+t) - 2*ln(2-x^2)

I tried

ics := u(x, 0) = 2*(1-ln(2-x^2));

bcs := D[1](u(0,t))=0, u(1,t)=2+ln(1+t);

PDEtools[Solve]([PDE,ics,bcs]); ## no solution

How do I do this?

Tom Dean

fn1 and gn1 are the 2x2 matrices. For example
 

fn1:=Matrix(2,2,[seq(f[i],i=1..4)])
gn1:=Matrix(2,2,[seq(g[i],i=1..4)])


eq1:=evalm(Matrix(2,2, [1, 1, 0, 3])&*fn1+gn1)=Matrix(2,2,[0]);

eq2:=evalm(fn1+gn1)=Matrix(2,2, [0]);

how can we find fn1 and gn1 matrices?

Consider the following code:

LM := [
   Matrix([[1,2],[3,4]]),
   Matrix([[5,6],[7,8]])
];
A := Matrix([[0,1],[1,0]]);
map(x -> A . x,LM);
A .~ LM;

where LM is a list of two matrices (just a test example), and A is some (test)matrix that I want to multiply onto each of these two matrices from the left, say. The map-construction works, as expected, but the elementwise operation .~ produces an error. Why?

I have a non-linear function to be optimized. It involves infinite sums. Maple plots the function so I can see where the minimum is. However the NLP solve keeps on evaluating without providing the solution. I have tried to write the function as a procedure but it does not work either.

I'd appreciate any suggestion

Greetings

How can I get LPSolve to output the unique subsets of {3,1,1,2,2,1} which sum to 5 (no recycling of set values)

partition.mw

(the code is a Yury/Love hybrid).

I have several maple worksheets (from the web) that have discussion blocks mixed within executable blocks.

All the executable blocks are delineated with a single '[' at the left while the discussion blocks do not.

How do I do this?

Tom Dean

So i got this code, im trying to iterate with jacobi and gaussseidel method.

H := HilbertMatrix(n, n, 1); b := Matrix(n, 1, proc (i) options operator, arrow; add(1/(i+j-1), j = 1 .. n) end proc); A := Matrix(n, 1, 1); Multiply(H, A); norm1H := norm(H, 1); norm2H := norm(H, 2); normHinf := norm(H, infinity); norm1b := norm(b, 1); norm2b := norm(b, 2); norminfb := norm(b, infinity); IterativeApproximate(H, initialapprox = Vector(n, 0), tolerance = 10^(-7), maxiterations = 10, method = gaussseidel)

But sadly no iteration gave me an answer, anyone knows wheres my mistake? i really help with this! 


thanks in advance