@Carl Love 

Thanks, I mixed it up with dolittle. No problem, simple transpose is all what is needed. Will fix now

@ANANDMUNAGALA 

updated answer for Crout

What does

    Pi= infinity Σ n=0

mean?

Is this Maple code?

@Carl Love 

Thank you, I did try that also. but it was giving me a strange result.

restart;
F:= w-> cos(w*L)*cos(w*x)/(w^6+w^2*k);
with(inttrans);
simplify(fouriercos(F(w),w,x));

@Axel Vogt 

Thanks. But it will really help to give an answer using a self contained code using the function I gave. I still do not see what you are doing.

This is what I typed.

-------------------------------
restart;
F:= w-> cos(w*L)*cos(w*x)/(w^6+w^2*k);
with(inttrans);
fouriercos(F(x),x,w);
---------------------------------

It does not evaluate as you can see. I do not know what the convert(%,Int) is for, but here it is:

convert(%,Int);

Mathematica has Inverse and Forward Fourier Transform functions for both Sin/Cos. I do not understand why Maple does not have these.

@Alejandro Jakubi 

thanks. I tried it. It works but the final result is hard to read. Needs lots of manually formating, and I am not going to waste my time formating code by hand so I can read it.

I prefer plain normal text based programming. All this 2D stuff and document mode and menus and click this and click that and such is making Maple hard to use for me.

@digerdiga 

I'd like to help, but I can't stand document mode. If you can make a version using worksheet that one can open and use that will help. I have not figured out how to open a file like yours, written in document mode as a worksheet. If I save the file as SAVE AS->Maple classical worksheet, then many things get mangled up.

May be there is a easy way, but I have not figured it out.

@Joe Riel 

Would you by any chance know to start a new "worksheet" or whatever the correct name in Maple, such that it uses this editor for the whole thing? I find having to insert code region inside worksheet an awkward way to use this. I'd prefer to have the whole worksheet be the editor, just like with Matlab, where one works in the editor all the time, and then there is the console where output goes to, and one can also still type commands on the console if needed.

When I saw the web page for this code editor page here http://www.maplesoft.com/products/maple/new_features/codeeditor.aspx it gave the impression that it is a whole new separate window/editor and not "code region" inside a worksheet. (may be that was just a marketing way to saying it).

When I run your code on my Maple 17.02 on windows, the plot is text based, not graphics. I wonder why. It does not look like anything like yours

@Alejandro Jakubi 

I can't even read the help pages. Here is an example:

http://www.maplesoft.com/products/maple/new_features/signal_processing.aspx

Hard to read the Maple commands. I am using firefox on windows. Terrible way to present something when the fonts/images are so bad and hard to read. And a page on signal processing no less.

May be Maplesoft managers do not read these pages and so they do not know. someone should tell someone at Maple to have a look one day.

@Joe Riel 

I see it now. thanks

@brian bovril 

I pasted some code into worksheet. Now each line had a ">" on its left.  Then selected all the code including the ">". Then did CTRL-F, put > in the search box, click replace all, then closed the find/replace window. The ">" are still there.

Maple 17.02, windows 7.

change "If" to "if" and add "()" around boolean expressions and add "*" between letters for multiplication (you also seem to have defind functions F and G but not using them).

----------------------------------

F:= x-> (-c-a*x)/b;
G:= x-> (-d-e*x)/f;
if (-a/b = -e/f) then
   print (false)
else
   print (true);
end if;

-------------------------------

does not Maple tell you on your system where the syntax errors are?

@Adri van der Meer 

thanks. But I am still a little confused.

according to help, it says on map, the following:
     "The map commands apply fcn to the operands or elements of expr."
my expression is this:

      lst:=[[12,8],[16,3]];

Now op(lst) returns the operands:

op(lst);
      [12,8] , [16,3]

So, there are two opeands, or 2 elements in the expression lst. So far so good.

Since map says the fcn is supposed to be applied automatically to each operand, why then map(igcd,lst) did not work?  Why one had to ask op for help first in extracting the operands? It seems to me that what op did, which is extracting the operands, is what map should have done?

@Andriy 

Yes, once the idea is clear, now one can find many ways to do it. Like this:

lst:=[[12,8],[16,3]]:
igcd~(op(lst));