Mainly retired from Eastman Kodak. (Of course, Kodak itself is now semi-retired. I don't think I had any influence in that.) I still write MATLAB code as I find something interesting, and I do attempt to write new tools to put on the File Exchange when I think I can make a contribution.

However, I DON'T answer direct e-mail questions, and I do NOT do homework. In general, your homework is YOUR problem, not mine. Please do NOT e-mail me with your homework problems or student projects. If I was willing to answer your questions, I would very rapidly become overwhelmed, because one question is never just one simple question.

When I'm not doing something with MATLAB, you might find me playing bridge, either in the club or online on BBO.

Professional Interests: MATLAB, numerical analysis, mathematical modeling

Avocational interests: Bridge, woodworking, woodturning

Answered

How can I plot a sin (x^2) function

First, there is ABSOLUTELY NO reason to predefine x as a sym. So this line is completely irrelevant: syms x; When you do creat...

How can I plot a sin (x^2) function

First, there is ABSOLUTELY NO reason to predefine x as a sym. So this line is completely irrelevant: syms x; When you do creat...

3 hours ago | 1

| accepted

Answered

Fixing code to do bisection properly

Um, look closely at the code that you wrote. Hey, it is your code. You should know what is in it. First, what does your functio...

Fixing code to do bisection properly

Um, look closely at the code that you wrote. Hey, it is your code. You should know what is in it. First, what does your functio...

4 hours ago | 0

Answered

difference between convolution code and 'conv' function

No difference at all, IF you are asking about the difference between a hand written code to do a numerical convolution computati...

difference between convolution code and 'conv' function

No difference at all, IF you are asking about the difference between a hand written code to do a numerical convolution computati...

19 hours ago | 0

| accepted

Answered

clock speed and cores

In general, this is impossible to answer. For ultimate speed, you want many cores, a fast clock speed, fast disk access, and a g...

clock speed and cores

In general, this is impossible to answer. For ultimate speed, you want many cores, a fast clock speed, fast disk access, and a g...

1 day ago | 0

Answered

Bisection method not working. Is my code heading in the right direction?

I had to laugh, as I recognized your coding style from a previous question. Then I remembered your name. Note that not(abs((end...

Bisection method not working. Is my code heading in the right direction?

I had to laugh, as I recognized your coding style from a previous question. Then I remembered your name. Note that not(abs((end...

1 day ago | 0

Answered

Variance between two vectors

The variance between two vectors is something that seems to be a bit ambiguous in meaning. I can think f at east a couple of way...

Variance between two vectors

The variance between two vectors is something that seems to be a bit ambiguous in meaning. I can think f at east a couple of way...

1 day ago | 0

Answered

why is my Fixed point iteration method only giving me first iteration. (i included info at the bottom of code and the next few iterations answers)

Think about it. Does p EVER change? Why not? How would you make it change? Inside the loop, we see only this: N=N+1; F=p; p...

why is my Fixed point iteration method only giving me first iteration. (i included info at the bottom of code and the next few iterations answers)

Think about it. Does p EVER change? Why not? How would you make it change? Inside the loop, we see only this: N=N+1; F=p; p...

2 days ago | 0

| accepted

Answered

Conversion to function_handle from double is not possible.

LOOK AT THE ERROR MESSAGE. Think about what it tells you. Error in RegulaFalsi_Mod (line 21) f(a) = (1/k)*f(a); What ...

Conversion to function_handle from double is not possible.

LOOK AT THE ERROR MESSAGE. Think about what it tells you. Error in RegulaFalsi_Mod (line 21) f(a) = (1/k)*f(a); What ...

2 days ago | 1

Answered

optimization toolbox

The optimization toolbox has existed for MANY years before 2006. I had a copy in my 1988 release. (Back then, it was called a ve...

optimization toolbox

The optimization toolbox has existed for MANY years before 2006. I had a copy in my 1988 release. (Back then, it was called a ve...

3 days ago | 1

Answered

Can PLSregress handle no intercept

It is NEVER a good idea to hack toolbox code. The exceptions are where you know enough to be able to write the modified algorit...

Can PLSregress handle no intercept

It is NEVER a good idea to hack toolbox code. The exceptions are where you know enough to be able to write the modified algorit...

3 days ago | 0

Answered

How can I count the largest number of repeated numbers in a double?

First, this is NOT a double vector. 00000000101100000011110000110000 It might be a string. But if you tried to write that vec...

How can I count the largest number of repeated numbers in a double?

First, this is NOT a double vector. 00000000101100000011110000110000 It might be a string. But if you tried to write that vec...

3 days ago | 1

Answered

A particular differential equation... why can't I solve it?

Easier than you might think to solve, even with pencil and paper. But sometimes a computer won't see the trick, at least, not wi...

A particular differential equation... why can't I solve it?

Easier than you might think to solve, even with pencil and paper. But sometimes a computer won't see the trick, at least, not wi...

3 days ago | 0

| accepted

Answered

could anyone help me to solve the issue

Easy enough. Just use setdiff. A=[2 3; 1 4; 2 4; 3 4]; That is, ... B = A; cind = [3 4]; for i ...

could anyone help me to solve the issue

Easy enough. Just use setdiff. A=[2 3; 1 4; 2 4; 3 4]; That is, ... B = A; cind = [3 4]; for i ...

4 days ago | 0

| accepted

Answered

Main diagonal operations problem

Easy enough, it seems. First, determine the diagonal elements. x = diag(rand(1,4)*.7 + .3); Next, you need to choose the other...

Main diagonal operations problem

Easy enough, it seems. First, determine the diagonal elements. x = diag(rand(1,4)*.7 + .3); Next, you need to choose the other...

4 days ago | 2

Answered

The accuracy of cumtrapz (numerical integration)????

This is something covered in any basic class that talks about numerical methods. (And why are you using methods about which you ...

The accuracy of cumtrapz (numerical integration)????

This is something covered in any basic class that talks about numerical methods. (And why are you using methods about which you ...

4 days ago | 2

| accepted

Answered

How to make a variable an even number?

If you just want to check if a variable X is an even number, then just check the value of mod(x,2). If the result is zero, then ...

How to make a variable an even number?

If you just want to check if a variable X is an even number, then just check the value of mod(x,2). If the result is zero, then ...

5 days ago | 1

Answered

Fitting of raw data

You don't have any simple model that you can use, at least, you have not suggested any. And, even if somebody suggests using pol...

Fitting of raw data

You don't have any simple model that you can use, at least, you have not suggested any. And, even if somebody suggests using pol...

5 days ago | 0

Answered

I don't know if matlab is calculating this properly

MATLAB is calculating what it calculated properly. The issue is, you told it to calculate the wrong thing. Computers are sooooo ...

I don't know if matlab is calculating this properly

MATLAB is calculating what it calculated properly. The issue is, you told it to calculate the wrong thing. Computers are sooooo ...

6 days ago | 3

Answered

Linear interpolation in loop

I totally agree with Bob. (By the way, this is not linear interpolation. Interpolation passes exactly through ALL points.) You ...

Linear interpolation in loop

I totally agree with Bob. (By the way, this is not linear interpolation. Interpolation passes exactly through ALL points.) You ...

7 days ago | 1

| accepted

Answered

Plotting a 2D crystal lattice from two primitive lattice vectors

First, you need to understand that MATLAB does not understand what you intend by this notation: v1= -0.5i + -sqrt(3)/2j Both i...

Plotting a 2D crystal lattice from two primitive lattice vectors

First, you need to understand that MATLAB does not understand what you intend by this notation: v1= -0.5i + -sqrt(3)/2j Both i...

8 days ago | 0

Answered

Can someone help me

Why do you assume an analytical Laplace transform exists? There are many functions for which this is not true. Even then, there ...

Can someone help me

Why do you assume an analytical Laplace transform exists? There are many functions for which this is not true. Even then, there ...

9 days ago | 2

Answered

Help on trigonometric operations using syms variable.

You probably need to use rewrite. syms t4 S = sind(t4-90) S = sin((pi*(t4 - 90))/180) rewrite(expand(S),'sincos') ans ...

Help on trigonometric operations using syms variable.

You probably need to use rewrite. syms t4 S = sind(t4-90) S = sin((pi*(t4 - 90))/180) rewrite(expand(S),'sincos') ans ...

10 days ago | 2

| accepted

Answered

How to improve the accuracy of lu decomposition?

Sorry, but it is often the case that people are sloppy in their work. That is my expectation here. You gave us A and x_sol, but...

How to improve the accuracy of lu decomposition?

Sorry, but it is often the case that people are sloppy in their work. That is my expectation here. You gave us A and x_sol, but...

10 days ago | 0

Answered

how can I perform curve fitting for 3 degree exponential equation?

Fit will do it, sort of. You will need to pose an explicit model for the three terms. But please consider if you really can do ...

how can I perform curve fitting for 3 degree exponential equation?

Fit will do it, sort of. You will need to pose an explicit model for the three terms. But please consider if you really can do ...

11 days ago | 0

| accepted

Answered

which way to call function is better?

Neither way is "better". Both feel just a little clumsy, since I'm not sure why you feel the need to package the variables toget...

which way to call function is better?

Neither way is "better". Both feel just a little clumsy, since I'm not sure why you feel the need to package the variables toget...

12 days ago | 2

| accepted

Answered

secant method trouble shooting

When you have a problem like this, plot the function. In fact, plot the function anyway! Next, it is a bad idea to use variable...

secant method trouble shooting

When you have a problem like this, plot the function. In fact, plot the function anyway! Next, it is a bad idea to use variable...

12 days ago | 0

| accepted

Answered

Rolling Beta For Multiple x and y variables simultaneously

Is the x vector equally spaced? If so, then my movingslope code (found on the File Exchange) will do it trivially and efficientl...

Rolling Beta For Multiple x and y variables simultaneously

Is the x vector equally spaced? If so, then my movingslope code (found on the File Exchange) will do it trivially and efficientl...

12 days ago | 1

Answered

How do I plot an equation with an imaginary quantity?

Look at what you have. t=[0:0.01:10]; y=exp(-5*i*t); y(1:5) ans = 1 + 0i 0.99875 - 0.049979i ...

How do I plot an equation with an imaginary quantity?

Look at what you have. t=[0:0.01:10]; y=exp(-5*i*t); y(1:5) ans = 1 + 0i 0.99875 - 0.049979i ...

13 days ago | 0

Answered

Generate random Bistochastic matrix with zero diagonals

Hmm. I notice that the re-normalization trick offered by dpb sometimes seems to not converge well. That begs the question of wh...

Generate random Bistochastic matrix with zero diagonals

Hmm. I notice that the re-normalization trick offered by dpb sometimes seems to not converge well. That begs the question of wh...

13 days ago | 0

Answered

matalb convolution with memory

This is not a job for convolution, since h lies on both sides of the equality, with a lag on h. Use filter instead. You should...

matalb convolution with memory

This is not a job for convolution, since h lies on both sides of the equality, with a lag on h. Use filter instead. You should...

14 days ago | 0