Community Profile

# John D'Errico

### Retired

Last seen: Today Active since 2006

A retired applied mathematician, consultant, sometime mechanical engineer, numerical analyst, woodworker, bridge player. Please DON'T contact me, as too many students then decide I am their personal consultant if I allow it. Anyway. I won't answer questions. PERIOD. I don't do your homework. I won't do consulting by mail.

All

#### Content Feed

View by

Inverse computation with a starting guess
It is often the case that someone THINKs they need a matrix inverse, yet they don't. And since you don't know if you can do this...

5 hours ago | 0

Vectorization approach to find matching elements in a matrix.
The part of your code that is a serious problem is where you build the arrays by growing them at each step. This is a really bad...

6 hours ago | 0

Please how can I calculate the values of this finite sum of numbers in fully sequential , I'm finding it difficult
Easy peasy. sumof4thpowers = @(N) nchoosek(N+1,5) + 11*nchoosek(N+2,5) + 11*nchoosek(N+3,5) + nchoosek(N+4,5); For example: s...

9 hours ago | 0

optimization in smart energy management system.I have used MILP algorithm but i got infeasible solution no point satisfy the constraints.how could i sort out these issues.
No solution means exactly what it sounds like. NO solution. And for this class of problem, if it tells you that NO solution exis...

18 hours ago | 1

How to find derivatives using the diff command
There are two versions of diff in MATLAB. One operates on numeric vectors and arrays. It simply subtracts each number from the o...

1 day ago | 0

Sum two matrices shifted by a non integer number?
There is no magic, whereby you can do things that take computational time, and do them instantly. What you want to do is implic...

1 day ago | 0

Falkner-Skan code
You write: tau = (mu/(2*delta_y))*(4*un,2-un,3) But what des that mean to you? Is the fragment (4*un,2-un,3) intended to be...

1 day ago | 0

| accepted

sqrtm vs chol to produce draws the multivariate normal distribution
Don't bother. sqrtm also can have a similar problem, but you won't like what it does either! For example: A = rand(3,2); A = A...

3 days ago | 1

| accepted

How to get tangent line from a point not on the curve
The curve is: syms x y x0 F(x) = 945.75*x.^0.1002; A line tangent to the curve at any spcific location x0 would be simple to ...

3 days ago | 1

| accepted

Solving systems of trig equations
My gut says no symbolic solution will exist. But that does not mean my gut is accurate here. First, I would throw pencil and pap...

3 days ago | 0

Merging piecewise polynomial (pp) structures
Is there any method to "merge" two PP functions? Well, no. Should there be one? I've been using and writing splines tools for 40...

4 days ago | 1

Solve of implicity equations
alpha = 0.14; A = 3300; n = 0.01 ; b_2019 = 69100; rho = 0.02; uno_alpha_A = ((1-alpha)*A); dos_rho = (2+rho); uno_n = (1+n);...

4 days ago | 0

What am I doing wrong?
ex2d is a function. You call a function, even a function you wrote yourself, by passing in arguments. For example, when you use ...

4 days ago | 0

Why can't the index variable of a for loop be stored in a structure or other array?
Surely this is not a big problem? for ind = 1:10 status.ind = ind; % do Stuff end

4 days ago | 0

is there any way to get a 2nd order polynominal through 2 points(each end point)
You want to fit a 2nd degree polynomial through two points? Of course only a straight line would work there through two points, ...

4 days ago | 1

Is there an option to overwrite the old release when installing a new one?
I'll admit, I typically have 3 to 6 releases, and usually at the high end of that scale. They are useful to me in case I might w...

4 days ago | 0

Volume of retinal lesion
It seems clear. You have multiple slices. Compute the area of each slice. At that point, you could just sum the area of each sli...

5 days ago | 1

Plot complex signal with imaginary and complex
n = 2000; t = linspace(0,500,n); Signal = cos(t/10) + i*sin(t/10); % a simple signal, as a function of time. plot3(t,real(S...

5 days ago | 0

Simplex optimisation using fminsearch
You are still looking for help on this? There is ABSOLUTELY no need to use a simplex optimizer for this problem. NONE AT ALL. L...

5 days ago | 0

3rd order approximation solving for variable
x = (0.2:0.2:0.8); y = [.3168, .7404, 1.3736, 2.2692]; format long g p = polyfit(x,y,3) So, to within floating point trash...

5 days ago | 0

| accepted

General variable creation question
Don't generate a zillion variables. Instead, learn to use arrays. Cell arrays, structs, multi-dimensional arrays. This will leav...

5 days ago | 0

How can I make a 3d plot with two equations and one variable?
By the way, there is NO reason to define f as symbolic, and then assign the number 0.0405 to f. f is now double preicsion. syms...

5 days ago | 0

Solving an equation gives me three possible answers. How can I reference only one of these answers?
syms v f eta; f=0.0394; f = 0.0394; eqnvel = (9.801-3.802*v)*v*eta==(((f*200*(pi/4)*0.1^2)/6.434)+((27*(pi/4)*0.1^2)/64.34))*...

6 days ago | 1

| accepted

Expanded use of .* operator
R2016b No need for access. That is when it appeared. Before that, bsxfun was needed. bsxfun appeared in R2007a. Before that, w...

6 days ago | 0

| accepted

Problem with convergence using newtons method
First, you CANNOT prove that Newton's method has quadratic convergence using any existing function, or any such piece of code. ...

6 days ago | 0

How to specify a domain of x and y values for a graph?
Did you read the help for fsurf? What does the second argument do? help fsurf syms x y f = y^2*(...

6 days ago | 0

| accepted

I would like to plot the possible solutions (x, y, z) of the following equation system: dn^2=(x-xn)^2+(y-yn)^2+(z-zn)^2 where n=1 to 3 and dn, xn, yn, zn are known
You have 3 spheres, with known centers and radii. Expand those equations. Then subtract them from each other. So maybe Eq1-Eq2...

6 days ago | 0

f(x)=0 matlab code for fixed point
You have been given a function f(x), and you now want to write a fixed point solver that will first, automatically create g(x) b...

6 days ago | 0

| accepted