I am trying to use Optimization[LSSolve] to fit the solution to a differential equations to data. I can solve my problem using Matlab, but I'd like to be able to use Maple as well. This is Maple 10. The proc is not getting the values of the parameters. > data := [[0,95], [11,425], [22, 928], [33,1358], [44,1589], [56,1683], [67,1724]]: > try2 := proc(K,alpha,r,IC) local DE1,R; print(K,alpha,r,IC): # for debugging DE1:=diff(y(t),t)=r/alpha*y(t)*(1-(y(t)/K)^alpha); R:=dsolve({DE1,y(0)=IC},numeric); map((d) -> rhs(R(d[1])[2])-d[2],data): end: > sol2 := Optimization[LSSolve](try2(K,alpha,r,IC), initialpoint = {r=.09, K=1750, IC=95, alpha=.3});

Hi, Sorry but i ve repeated the same question a few days ago but because it was quite down the line i thought to give it another shot I would like to make a 3D plot of some density functions(F(X)) where each lets say cdf, corresponds to a time value.So far Ive tried using the transform((x, y) command.The example below where I show how I have worked out for 4 cdfs corresponding to 10,20,30 and 40 ms, hopefully it will make it clearer of what i mean.(I have a thousand of different distributions up to 10 secs). with(stats): > C001:=plot(statevalf[cdf,lognormal[0.3792,0.9556]], 0..4, colour=green):

I am using the latex2edu conversion service at qa.brownstone.net/latex2edu to convert latex files into MapleTA questions. However, for the last 4 days this service has not been available. Mail to support@brownstone.net and my local MapleSoft representative has not been answered. Has anyone else experienced this problem? Does anyone know whether the service is available elsewhere? Arne Jensen Aalborg University Denmark

And now for something completely different - something I don't know how to do in Maple! [ :-) ]. Given a system of linear inequalities, say like {a1>5, a2<100, a1 >< a2}, how do I find sample values (integers in this case) in the solution space? For example, I would be happy with {a1=6,a2=3}.

Hi,

we've been happily using Brownstone's latex -> edu conversion service for some time now, but always with that slightly queasy feeling of "what happens if it's unavailable?". Well that's the case at the moment: I've submitted a ticket with their support folk, but times in Australia don't match up well with US times! In the meantime we have a number of assignments due to be made available to students, but which only exist as latex files. Conversion by hand is, to say the least, painful and time consuming (our Level I coordinator has already converted two by hand, and doesn't need another few hours of such nonsense), so I'm hoping that there may be another way in which this can be done. Any suggestions?

I'm currently outputting some animations, but I'm not using the animate command - I'm using the command display([seq(tmp_alpha_plot||kk,kk=1..nt)],insequence=true) having previously created tmp_alpha_plot||ii:=display([ plot([seq([r_arr[i],alpha1_arr[i,ii]],i=0..nr)]), plot([seq([r_arr [i],alpha2_arr[i,ii]],i=0..nr)], color=magenta)]): What would be useful for diagnostics purposes is to be able to print in the corner of the animation both the timestep/frame currently displayed (I know this is automatic when you use the animate command, and I'll change if I need to), and the actual numerical value of the variables alpha1_arr and alpha2_arr. It would also be useful if I could also show the numerical values of some other variables I'm working with, as they're defined at that time.

Hey I was wondering which of the many Root() commands I would use to obtain an Integer root from a large polynomial for example X^2 - 9 would return 3 I will be using it for very large polynomials with big coefficents so if there is a more efficent one then that would be helpfull. Cheers in advance

I want to calculate the normal derivative with Maple?What's the command?Waiting for your help!

Hi - I want to graph an ode and the following are my commands: > restart; > with(plots); > with(DEtools); > plot1 := DEplot((diff(y(t), t, t))+6*(diff(y(t), t)) = 80*exp(2*t), {y(t)}, t = -1 .. 1, y = -1 .. 1, [[y(0) = 4], [(D(y))(0) = 3]], linecolor = ([BLUE, GREEN, BLACK]), arrows = MEDIUM); The error I get is : Error, (in DEtools/DEplot/CheckInitial) Too few initial conditions: [y(0) = 4] Can anybody tell me what I'm doing wrong? Thanks.

Hello, I wonder if someone could explain why, in the attached file, when I collapse the execution groups connected with the 3 equations the label (2) disappears but not the labels (1) or (3). Is there a workaround for this? Thank you View 2292_Label (2) disappears.mw on MapleNet or Download 2292_Label (2) disappears.mw
View file

I have this problem with the sum of a recursive seq: > restart; > with(student): > with(plots): Warning, the name changecoords has been redefined > a:=5.5*10^(-3): > R:=10*10^(-3): > d:=n-> ((-1)^n)*(R-(a^2)/(R+(-1)^(n-1)*d(n-1))): > d(0):=R: > afs:=seq(d(n),n=0..100): > q:=n->(-1*a*q(n-1)/(R+(-1)^(n-1)*d(n-1))): > q(0):=a*4*3.14*8.85*10^(-12): > lad:=seq(q(n),n=0..100): > sum(q(n)/(sqrt(y^2+(z-d(n))^2))/(4*3.14*8.85*10^(-12)), n=1..100); Error, (in q) too many levels of recursion How can I avoid this, is the syntax wrong? An additional question, if I want to plot the seq's how do I do this? I'm a bit over my head here when it comes to maple knowlegde.

Hi - Its me again. Let me try to explain my question with the following: y := proc (t) options operator, arrow; a*exp(2*t)+b*t*exp(2*t)+c*t^2*exp(2*t) end proc solve({3*c = -9, 3*b+12*c = 8, 3*a+6*b+2*c = 4}, {a, b, c}) It would be really cool if there was some way I could take the output from solve (values for constant a, b and c) and automatically assign these values in expression y. I can manually take these and type out y all over again but I think there's got to be a way that Maple recognizes that, alright, here's a, b and c and I'll just assign their values to a, b and c in expression y.

I have a class with about 400 students. For the new term I opened a new class which will have about the same audience like the one before. Therefore I would like to copy all student accounts but start with empty gradebooks. As far as I saw it's possible to upload a class roster but it's not possible to download it. Could anybody tell me an easy solution for this problem?

Hi, how can I stop a procedure of maple a use the global variables create before stop the code? Thanks Marco

Hi.. I'm a quite new user of maple, so I've got a question: I'm trying to solve a ode numeric with maple. I've got a solution, a procedure where I a can get a value of y(x) for any x.. (and also plot the "function") Additionally I would like to have a procedure which can give me the slope of the graph - y'(x) for any x. Can anyone help me to solve this problem? My code so far is: lignsub:=diff(y(x), x, x) = 0.7904982116e-3*sqrt(1+(diff(y(x), x))^2)/(1+0.1012020000e-2*sqrt(1+(diff(y(x), x))^2)) `løsning` := dsolve({lignsub, y(0) = 0, (D(y))(0) = 0}, type = numeric, output = listprocedure)