I have a piecewise defined function iL1. Now I would like to define a function iL2 = ( iL1(t+2/3*T) for 0 <= t <= 1/3*T iL1(t-1/3*T) for 1/3*T < t <= T ) How can I do this in Maple? How can I change the boundaries in a piecewise defined function in order to be able to perform a time shift? The following linked worksheet shows the definition for iL1 and what I would like to achieve. Download 305_test_piecewise_timeshift.mws

I experienced a problem with OpenMaple: I want to use it in an application together with GMP (and other libraries) under C++. Now, when I just create gmp-integers from 1 to 100 and print them out, there is a strange behaviour if I do the same, but with starting a Maple-instance before: The numbers 8, 9 and everything greater equal 64 are not printed, instead some junk appears on the screen. But it looks like the number are internally correct, one can still calculate with them. Maybe MAPLE overwrites some global constants of GMP? Has somebody observed such problems and knows a fix? I'm using MAPLE 10, gmp-4.1.4 and gcc-3.3

Recently I was skimming on the Internet for algorithms for alignment of biological sequences (e.g. DNA sequences, protein sequences). The usual purpose of such comparisons is to determine the evolutionary genetic history of two sequences, which in the case of DNA you can think of as strings over the alphabet {A,G,C,T} (the nucleobases). Differences between sequences can arise though substitutions (A → G, G → A, C → T, or T → C), insertions, or deletions. In the last two cases, the resulting string differs in length from the original, and dealing with these size differences is the chief problem that alignment algorithms face. The Needleman-Wunsch algorithm is a very straightforward global sequence alignment algorithm, first proposed in 1970. It's a good example of the applicability of dynamic programming towards biological problems. Following is a Maple implementation based on this Ruby implementation. The algorithm depends on a similarity matrix which measures the "mutational distance" between two nucleobases. It takes two input strings, and returns two strings with dashes ("-") inserted to indicate insertion or deletion events between the sequences.

Suppose I want to know whether or not z is in the range of 1..10. Is there any way of finding this out without writing z twice, e. g. if (z>=1) and (z

I have recently upgraded to Maple 11 and I have been having a recurrent problem when I run documents with embedded plots that use "fieldplot". These documents functioned perfect using Maple 10. My problem is that I get the error "Invalid Plot Structure" when I run from a button the code restart: with(plots): plotxy:=fieldplot([x,y],x=-2..2,y=-2..2): SetProperty('Plot0', 'value', plotxy): end use; for the first time after the document is first opened. If the code is run again (after acknowledging the error, and without closing out the document), the field plot will be displayed as it should.

I defined two matrices; vectorC:=Matrix(1,4,symbol=c); vectorE:=Matrix([[e[0]],[e[1]],[e[2]],[e[3]]]); so I wrote the equation: eq3:=evalm(vectorC&*vectorE)=0; I solved e[1] as follows: solve(eq3,e[1]); However, it shows: Error, (in solve) invalid arguments, How can I solve the problems? Thanks

Can someone help me to plot the familiar bell shaped normal distribution curve ? I'm using the NormalPlot command with a RandomVariable, but all I'm getting is a scatter plot of points along the line, y = x. Can anyone help me to get the gaussian, bell shaped curve ? v/r,

I'm trying to enter a piecewise function with three parts, but can't seem to figure out how to make an extra line (the "expression" selections have one, but it only allows two lines to be entered). I read somewhere else to do CTRL+Shift+R, but that didn't seem to work. What I'm trying to enter is
y=x^2-5 x<2
-3x-7 -x<-2 and x<5
abs value(x-2) 5<x

I can get the first two lines, but can't find a way to include the third line. I'm supposed to graph the function so that all 3 are shown (may not be connected), and then I have to evaluate f(-3),f(3), and f(7). I haven't tried the evaluate part, because I can't seem to get all three components to graph. I have Maple 10. I'm new to Maple, any help would be greatly appreciated, this is quite the learning experience!

Hi all, I was just wondering if there is a way to change the window settings for the display function. I am making a picture using Maple 9 for a school project and have many plots on the same graph, but I only see portions of it. Thanks for your help.

I have a problem that I'm somewhat confused about and I'm looking for help if possible. The problem states: Given the equation (x^2-1)/(x-1)=0, graph on an appropriate scale. What is the domain, what is the solution, and what are the x and y intercepts? I'm not sure how to put this one into Maple, if anyone can help I'd greatly appreciate it! This is my first time using Maple and I'm a little confused. I have Maple 10. Thank you! Stacey

Help please! I have a problem that states "use implicit plot to graph 3x^2 + 2xy + y^2 = 4 on the interval [-5,5]". I'm not sure how to do an implicit graph, or how to identify the interval to plot. If anyone can help me, please let me know, I've been working on this one for hours (I have Maple 10). Thank you, Stacey

Is the Statistics package for Maple 11 able to handle, and or deal with conditional probabilities ? I may have to start using Bayes' Theorem really soon. Can you also create your own sample spaces and assign probabilities to the simple events in that space ? I would prefer accomplishing the latter using Maple's pre-defined packages instead of having to create my own, if it is possible to do so. thanks, v/r,

Hi, I have numerous terms all of the type opl[i]:= ... with i an integer . Now I want to calculate the total sum of all these terms. I've already worked with a commando like: "for i from 1 to 4 do opl[i]:=... end do;" but I can't find a way to make the sum of multiple terms. By the way is there a site of a good handbook where I could learn more about programming in Maple and make animated curves or objects?

I thought I would take up a little time to write about some of the new plotting features in Maple 11. This post will give an overview of what is in Maple 11 plotting - subsequent posts will give more details on the individual features. Let's look at each of the features mentioned on the plotting enhancements page. All of these are only for 2-D plots.

The attached Maple10 worksheet solves a system of differential equations. A plot of the solutions y1(t) (red curve) and y2(t) (green curve) appears below. Download 2353_fsolve-avoid.mws
View file details As a check, I want to use fsolve and the "avoid" option to find both times at which y1 has the same value as when y1=y2, but I'm having trouble. I would appreciate any advice on how I can get the "avoid" option to work for me. fsolve finds the earlier time easily enough (t=.9036852930e-1), but when I use the "avoid" option (highlighted statement below) to find the later time (t=.5225499740), I get the error "Error, (in fsolve) avoid = {.9036852930e-1} is an invalid option." The code is appended below. Thanks. Glenn ======== > restart; > with(plots): > sys:= Diff(y1(t),t)=-3*y1(t) + 2*y2(t), Diff(y2(t),t)=y1(t)-3*y2(t); > ic := y1(0)=1, y2(0)=5; d d sys := -- y1(t) = -3 y1(t) + 2 y2(t), -- y2(t) = y1(t) - 3 y2(t) dt dt ic := y1(0) = 1, y2(0) = 5 > odesol:=dsolve({sys,ic},{y1(t),y2(t)},type=numeric,output=listprocedure); odesol := [t = (proc(t) ... end proc), y1(t) = (proc(t) ... end proc), y2(t) = (proc(t) ... end proc)] > Y1:=eval(y1(t),odesol); Y2:=eval(y2(t),odesol); Y1 := proc(t) ... end proc Y2 := proc(t) ... end proc > Teq:=fsolve('Y1'(t)='Y2'(t),t); Teq := 0.5225499740 > Yeq:=Y1(Teq); Y2(Teq); Yeq := 1.45974175440983366 1.45974175447049736 > T1:=fsolve('Y1'(t)=Yeq); T2:=fsolve('Y2'(t)=Yeq); T1 := 0.09036852930 T2 := 0.5225499740 > fsolve('Y1'(t)=Yeq, t, avoid={T1}); fsolve('Y1'(t)=Yeq, t=T2); Error, (in fsolve) avoid = {.9036852930e-1} is an invalid option 0.5225499740 > Y1(%); 1.45974175440983366 > plot({Y1(t),Y2(t)},t=0..3);