fsolve(x^2/20-10*x-15*cos(x+15),x);

1.274092075

I know to get maple to give me the other root I have graphically examine the funtion to see that it's around 200.  I need to use fsolve using x=190..210 to find it.  But that's a bit of a pain, I don't want to have to graph the function everytime, if my limits aren't far enough I may not find it. 

Wondering if anyone can help me with this long winded problem...

 

So... What I have to do is randomly select 6 numbers from the interval [0,1], not hugely important how many decimal places there are but I've used 10.

Ive written some stuff to do this

T := rand(0..1000000000000):

t := proc()
  evalf(T()/1000000000000)
end proc:

C := proc(n)
local S,i:
S := NULL:
  for i from 1 to n do
    S := S, t(i)
  end do
end proc:
 

how can i get the inetersection point (a real number of x axis) of a complex plot and a real plot? and how to assign them respectively to an array?

for example the junctions of lines:

d1:=(sec(I*(x^2-2)))(exp(I*x)*sin(2*x)+cos(I*(2*x-3)));       and

d2:=0.5;  by  plot([d1, d2], x = -5 .. 5);

 

I'm programming really much right now in Maple. To do so, I simply write foo :=proc(...) ... end proc statements directly into a Maple-document.

However, you neither have syntax-highlighting (save Standard-Keywords such as 'proc' or 'for', 'do' etc.) nor auto-indent (you actually have to use the space bar, tab won't work).

You certainly agree that that is very annoying. Are there other options? Is there e.g. some sort of Editor where I can program most of my stuff without entering maple and will it load packages (with(...) in Maple)?

Thanks in advance!

This is similar to the roots of trig function thread that was just posted.  

When I Set_EnvAllSolutions:=true; and then use cos(x)=0 for example are the _B3~ and _Z3~ just place holders?  

If I want the actual numeric values over a range and I'm trying to use fsolve, I can't seem to get it to work. 

SetEnvAllSolutions:=true; fsolve(cos(x)=0, x=0..3*Pi);

I get 4.712388980 as the only answer.  I should have 3 values shouldn't I?

 

I'm currently trying to plot the Ferrers graph of a plane partition, i.e. a 3D structure of stacked boxes. My input would be a plane partition p_ij with the information ( i, j, k ) for each 'stack of boxes' where i, j are the coordinates in the matrix and k is the value of the respective entry. As an example, the following input 3 2 1 1 2 2 1 1 1 should produce a plot that looks like this: Thanks a lot in advance! Cheers Nils

I've been using a combination of Fortran and gnuplots to do my plots. my research required an intensive amount of data plotting. I've been trying to use maple to do scatterplots and curve fitting for my data. here is the code: > K := readdata("/Users/xxx/Work/dt40s5l10.110/ALL/dt40s5l10.110.1.DFIR.AVE", 8) # the file contains 2972 lines and 8 columns of data > unassign('G', 'X'); > seq(assign(G[i], (K[i+1, 5]-K[i, 5])/(K[i+1, 1]-K[i, 1])), i = 1 .. 2970); # to calculate the derivative > seq(assign(X[i], K[i, 1]), i = 1 .. 2970); when I tried to plot X vs G using scatterplots

Hi!

I'm trying to solve a differential equation of the form y'=cos(pi*t*y).

I wrote:

<de:= diff(y(t),t)=cos(pi*t*y(t))

And that works fine.

Then I tried to use dsolve(de) but that shows nothing and just goes to the next line

So I figured maybe it couldn't solve it and decided to try the transform methods. For instance,

with(inttrans)

ics:=y(0)=0

dsolve({de,ics},y(t),method=laplace) <- with regards to this, I also tried Fourier

How can i solve this complicated problem????restart:
R:=3:
kappa:=2.4048/R:
a:=-mu^2/kappa^2:
sigma1:=Int(t^2*BesselJ(0,t)*BesselJ(1,t)*exp(a*t^2),t=0..2.4048):
sigma2:=Int(t^3*BesselJ(0,t)^2*exp(a*t^2),t=0..2.4048):
sigma3:=int(t*BesselJ(0,t)^2*exp(a*t^2),t=0..2.4048):
sigma4:=Int(t*BesselJ(0,t)^2*BesselJ(0,(q*t)/kappa)*exp(a*t^2),t=0..2.4048):
s[q]:=(sigma4/sigma3)^2:
epsilon:=lambda^2/2+kappa^2+mu^2*(2-(2*sigma1-(mu/kappa)2*sigma2)/sigma3):
h:=(s[q]*exp(q^2/(2*lambda^2))*erfc(q/(sqrt(2)*lambda))):

Hello,

I'm trying to do a scatter plot but with a logarithmic vertical scale.  Does anyone know how to do this?  I'm currently using the pointplot() function and giving it a list of points.

 

Thanks,

Andy

How to convert a data surface into a region:

surfdata([[1, 1, .69], [1, 2, .48], [2, 1, .37], [2, 2, .44]], axes = frame, labels = [x, y, z], filled = true)

I've tried to fill the volume below using the option "filled " without any results. I will appreciate any suggestion.

Hello,

can you tell me, please, how I can create labels when plot command (and its labels option) doesn't work?

 

Try these examples:

Ex1:  plot(1/(3+exp(x)), x = 0 .. 10, labels = ["x", "y"]);

Ex2:  plot(1/(3+exp(-x)), x = 0 .. 10, labels = ["x", "y"]);
 

In first case everything is OK, in second "y" label doesn't appear and I can't manage to put it into plot.

Thank you.

 

PS: I'm using Maple 13.

I have defined a sparse matrix and verified that it is space by using the command- sparsematrixplot(Kd,'matrixview');

first, is there another way to check wether the matrix is stored in the sparse form?

second, when trying to sovle an algebraic set of equation using the procedures for hendling the sparse matrix i encountered the following error:

Error, (in HWcall[1]) column vector dimension less than length of sparse data
 

what is the meaning of this error and how can I repair it?

Dvir
 

 

Hi, there were already some topics about labeling plots but I have still a question:

I would like to label my axis with e.g.: rho [kg/m^3]. roh should be displayed symbolic (2D output) as well as m^3. On the other hand I would like to avoid that Maple builds a (2D) fraction with kg/m^2. Does anybody know how to do this?

The other thing: Is it possible to add an legend to a 3D plot?

Thanks!

how can i calculate all roots(or some roots in a given interval, or the J th root) of a (real or complex) trigonometric periodic equation?