Comments and Ratings

This is much faster, and has no requirement as to types or ordering:

ii = 1;
while ii<size(A,1)
    tf = true(size(A,1),1);
    for jj = ii:size(A,1)
        if isequal(A(ii,:),A(jj,:)) && ii~=jj
            tf(jj) = 0;
        end
    end
    A = A(tf,:);
    ii = ii + 1;
end

For example, try it with:
A = {1,'red',magic(3);'blue',magic(4),'green';1,'red',magic(3)}

Cyberfrog,
GUI_34 has a custom menubar. None of the GUIs in this collection are programmed with GUIDE. Please email me if you have more questions.

Actually atan2 is fine. The source of the error in your example is this line:

[x2 y2 z2] = sph2cart(t2, p2, r2); % step 5b

The value of y2 should equal zero but there are round-off errors occurring. If you simply assign y2 = 0 after this line in your example, you will see that t1 is back at zero as expected. This has nothing to do with atan2.

Three comments.

1. I do not see where you mention that this file needs the symbolic toolbox to run. It does. I wouldn't have downloaded it if I had known that because I don't have the symbolic toolbox. You should update the description to reflect this requirement.

2. Rating your own files is not good behavior. Everyone expects that an author likes his own files or he wouldn't have submitted them. Please don't try to trick people into thinking your file is esteemed by the community. Let others rate your files, then your skills and file quality can grow, and people looking for a quality file won't be misled by false ratings.

3. A numerical alternative that works for Shahrir's problem is:
f = @(x) cosh(x).*cos(x) + 1;
rt = newtzero(f); % Found on the File Exchange.
x = -8:.1:8;
plot(x,f(x),'b',rt,f(rt),'*r')

Tanguy,
You are correct. I didn't mention resizable GUIs because I make all of my own GUIs non-resizable. Is there something you have in mind in particular? If so email me with your suggestion, and I may include it in the next update. Thanks.

Works well. One thing I would like to see added is an optional second output which stores whether the dialog terminated due to timeout or user acknowledgment. That way if this was called from a code which was running automatically, the user could be notified that a warning was thrown earlier. This could be useful for debugging.

Thanks someone,
BUG! I should be using listboxtop, not value. Thanks, I will update that too.
Also, the suggestions for added exercises are coming into my email. Keep them coming! Any other GUI suggestions are welcome, simple or complex.

Juliette,

You cannot index the output of npermutek(1:5,5) into an array (nuc) of length 4! The index exceeds dimensions because the output of the above call to NPERMUTEK will have many '5's in it. You will get the same error by doing this:

nuc = 'ABCD'
R = nuc(5) % or nuc(6) or nuc(7), etc. nuc only has 4 elements!

This is not an NPERMUTEK problem, it a user misunderstanding the output problem.

Something is still wrong.

[D,G,S] = LCS('fbce','abcde');S
S =
bce

Bruno,
In all fairness to Oleg, your code does not solve the same problem. Try this:

N = 10; c = {};
for ii = 1:N
    if round(rand)
       c = {num2cell(1:ii) c};
    else
        c = {mat2cell(ones(ii),ii,ii),c};
    end
end

tic
M1 = flatcell(c);
toc

tic
M2 = fcellrecurse(c);
toc

isequal(M1,M2) % Not necessarily.

To Oleg: In actuality, what Bruno's code does seems to me more useful. What is the purpose for flattening EVERYTHING in the cell?

Another reason I bypassed kron is illustrated here:

http://i217.photobucket.com/albums/cc229/spamanon/tester_output-2.png

Urs, I think you are correct about mentioning kron in the see also.

Mayowa,

Indeed that is faster for smaller A,B. However, for larger A and B, kronecker is several times faster.

Cool idea! Works well.

Siyi, Your code, as written, will work when there is one submatrix within A. However, it will need to be modified for multiple submatrices. Even without the modification, it is slower to use the approach you suggest. In the following, I made a copy of findsubmat, took out all of the NaN checking and the transposing of A, and put your code in the proper place. Here are some results:

A = round(rand(2000)*1);
B = [1 2;4 5];
A(30:31,60:61) = B; % Only one submatrix.

tic
C = findsubmat(A,B); % Elapsed time is 0.055974 seconds.
toc

tic
C = findsubmat_2(A,B); % Elapsed time is 0.502954 seconds.
toc

Also, this is really the kind of exchange that better takes place through email. Mine is listed in the help for findsubmat.
Feel free to email me about this code.
Thanks.

An update has been submitted which allows user to enter the set to be partitioned. For example: partitions(['a','b','c'])

Wolfgang, I realized that after I posted the screenshot. I believe that file was my first experience with that author. Funny to see it now and everyone knows where it came from :-).

nargin is a built-in Matlab function, as is numel. Take these questions to email please! My email address is in the file.

4.4 on my Dell XPS M1710, ver 2007a, Vista Home Premium, Invidia GeForce Go 7900 GS

Please don't rate your own files.
Everyone assumes the author likes the file.

This seems a little faster for large number of dimensions:

IDX = find(A);
L = size(A);
[x{1:length(L)}] = ind2sub(L,IDX);
sub = [x{:}];
if nargout==2
    v = A(IDX);
end

Something to consider.

I realized I should have a check that the input arg is a handle of the correct type. Update on the way.

In the lines where you calculate the quadratic formula, you need a parenthesis around the 2*a, like this;
x1=(-b+sqrt(b^2-4*a*c))/(2*a)
x2=(-b-sqrt(b^2-4*a*c))/(2*a)
Otherwise your solutions are only valid for a = 1.

This is a script that issues a clear command - NOT A GOOD IDEA. Please only submit functions to the file exchange. There is no need for 3 separate files to make the game you want. Everything you are using the ask.m and check.m functions to do can be done in a simple while loop!