3.162 total contributions since 2007

Interested in external applications, mex & engine applications with Fortran, C, C++, etc, particularly in the areas of speed improvement and memory management efficiencies.

Answered

How to sort a matrix in a random order

E.g., >> x = reshape(1:24,4,6) x = 1 5 9 13 17 21 2 6 10 14 18 22 3 7 ...

How to sort a matrix in a random order

E.g., >> x = reshape(1:24,4,6) x = 1 5 9 13 17 21 2 6 10 14 18 22 3 7 ...

alrededor de 20 horas ago | 0

| accepted

Answered

can someone help me to solve this error (This statement is not inside any function. (It follows the END that terminates the definition of the function "RK4".)

Looks like your plot( etc ) and legend( etc ) lines are stuck inbetween two functions. Did you intend them to be inside a scrip...

can someone help me to solve this error (This statement is not inside any function. (It follows the END that terminates the definition of the function "RK4".)

Looks like your plot( etc ) and legend( etc ) lines are stuck inbetween two functions. Did you intend them to be inside a scrip...

alrededor de 21 horas ago | 0

Answered

making array in for loop

Maybe prior to the loop define Alignment as a cell array to hold your results: Alignment = cell(10000,1);

making array in for loop

Maybe prior to the loop define Alignment as a cell array to hold your results: Alignment = cell(10000,1);

4 días ago | 0

| accepted

Answered

Why is my function generating a vector? Please help?

A uniform random number between a and b is: u = a + rand*(b-a) Generate two of these with a = -0.05 and b = 0.05 and then do t...

Why is my function generating a vector? Please help?

A uniform random number between a and b is: u = a + rand*(b-a) Generate two of these with a = -0.05 and b = 0.05 and then do t...

6 días ago | 0

Answered

Can someone explain the concept of B= null(A) in simple words?

The columns of B form basis vectors for the "null space" of A. Any linear combination of the B columns, when multiplied by A, w...

Can someone explain the concept of B= null(A) in simple words?

The columns of B form basis vectors for the "null space" of A. Any linear combination of the B columns, when multiplied by A, w...

6 días ago | 1

Answered

Unable to perform assignment because the size of the left side is 1-by-1 and the size of the right side is 5-by-1.

Looks like Cd and rho are vectors. So the result of those calculations on the right hand side will be vectors. But you are tryin...

Unable to perform assignment because the size of the left side is 1-by-1 and the size of the right side is 5-by-1.

Looks like Cd and rho are vectors. So the result of those calculations on the right hand side will be vectors. But you are tryin...

6 días ago | 0

| accepted

Answered

What can I use instead of VPA?

See the following FEX submissions by John D'Errico: https://www.mathworks.com/matlabcentral/fileexchange/36534-hpf-a-big-decima...

What can I use instead of VPA?

See the following FEX submissions by John D'Errico: https://www.mathworks.com/matlabcentral/fileexchange/36534-hpf-a-big-decima...

7 días ago | 1

Answered

Two equal datetime values showing as unequal when compared?

This is often just a display issue. E.g., >> s = datetime(2017,3,24,9,5,0.1) s = 24-Mar-2017 09:05:00 >> t = datetime(20...

Two equal datetime values showing as unequal when compared?

This is often just a display issue. E.g., >> s = datetime(2017,3,24,9,5,0.1) s = 24-Mar-2017 09:05:00 >> t = datetime(20...

7 días ago | 0

Answered

How can I multiply square submatrices more efficiently?

Here is the naive mex code. Could probably be made faster by doing the for-loop in parallel or trying to optimize cache hits, b...

How can I multiply square submatrices more efficiently?

Here is the naive mex code. Could probably be made faster by doing the for-loop in parallel or trying to optimize cache hits, b...

7 días ago | 1

| accepted

Answered

big difference in execution time in approximately the same code

The matrix multiply operator * in MATLAB calls highly optimized compiled BLAS library code in the background. The BLAS library ...

big difference in execution time in approximately the same code

The matrix multiply operator * in MATLAB calls highly optimized compiled BLAS library code in the background. The BLAS library ...

7 días ago | 2

| accepted

Answered

This is a question from numerical methods asking to explain what the code does for 12 marks.

https://en.wikipedia.org/wiki/Newton%27s_method

This is a question from numerical methods asking to explain what the code does for 12 marks.

https://en.wikipedia.org/wiki/Newton%27s_method

8 días ago | 0

| accepted

Answered

how to use passing by value and passing by reference

In the background, MATLAB effectively passes all variables into functions as shared data copies. This has the same effect as pas...

how to use passing by value and passing by reference

In the background, MATLAB effectively passes all variables into functions as shared data copies. This has the same effect as pas...

13 días ago | 0

Answered

matrix dimensions must agree error

Well, tau is only two elements and w is a lot more than that, so you can't multiply them element-wise like you are attempting to...

matrix dimensions must agree error

Well, tau is only two elements and w is a lot more than that, so you can't multiply them element-wise like you are attempting to...

13 días ago | 0

| accepted

Answered

Monte Carlo integration of sin(x)

This is what you are currently doing with the "random counting in an area" method: To make it go from -1 to +1 instead, you c...

Monte Carlo integration of sin(x)

This is what you are currently doing with the "random counting in an area" method: To make it go from -1 to +1 instead, you c...

13 días ago | 0

Answered

How can i split the values of vector into group of 5?

Assuming your vector is divisible by 5, you could do this: v = your vector g = reshape(v,5,[]); % reshape v into a 5 x N matri...

How can i split the values of vector into group of 5?

Assuming your vector is divisible by 5, you could do this: v = your vector g = reshape(v,5,[]); % reshape v into a 5 x N matri...

14 días ago | 0

Answered

Runge-Kutta function with a second order ODE

I'm still confused about your initial conditions and what y0 and y1 are. But let's back up a bit. For a 2nd order ODE, your st...

Runge-Kutta function with a second order ODE

I'm still confused about your initial conditions and what y0 and y1 are. But let's back up a bit. For a 2nd order ODE, your st...

14 días ago | 0

| accepted

Answered

how to call struct array in subroutine

See this related thread: https://www.mathworks.com/matlabcentral/answers/480829-what-is-the-best-way-to-insure-that-all-of-my-f...

how to call struct array in subroutine

See this related thread: https://www.mathworks.com/matlabcentral/answers/480829-what-is-the-best-way-to-insure-that-all-of-my-f...

14 días ago | 0

Answered

For and if loops

RT([3,6],v) is the syntax to use for the 3rd and 6th rows of the v'th column of RT. But this is of course going to be two eleme...

For and if loops

RT([3,6],v) is the syntax to use for the 3rd and 6th rows of the v'th column of RT. But this is of course going to be two eleme...

15 días ago | 0

Answered

Matlab Solving Differential Equation using Runge Kutta

In this line: k_1 = F_xy(x(i), xy(i)); You've got xy(i) as an argument, but there is no xy variable. I think you meant y(i) h...

Matlab Solving Differential Equation using Runge Kutta

In this line: k_1 = F_xy(x(i), xy(i)); You've got xy(i) as an argument, but there is no xy variable. I think you meant y(i) h...

18 días ago | 0

Answered

How can i write MacLaurin Sequence for 1/1+x2

An outline of your code could look something like this: x = _____; % you fill in the blank here, or get an input from the user ...

How can i write MacLaurin Sequence for 1/1+x2

An outline of your code could look something like this: x = _____; % you fill in the blank here, or get an input from the user ...

18 días ago | 0

Answered

Multiplication of complex matrices

There are two different operators: ' is the complex conjugate transpose .' is the straight transpose (no complex conjugate), n...

Multiplication of complex matrices

There are two different operators: ' is the complex conjugate transpose .' is the straight transpose (no complex conjugate), n...

19 días ago | 1

| accepted

Answered

why MATLAB 2019b MEX file is slower

Impossible to say for sure without knowing what your mex function is doing. One thing to note is that in R2018a the storage for...

why MATLAB 2019b MEX file is slower

Impossible to say for sure without knowing what your mex function is doing. One thing to note is that in R2018a the storage for...

19 días ago | 3

Answered

Different complex log values displayed for a "single" input

This is simply the result of floating point calculation artifacts ... not a bug. For instance, your assertion that all of your ...

Different complex log values displayed for a "single" input

This is simply the result of floating point calculation artifacts ... not a bug. For instance, your assertion that all of your ...

19 días ago | 0

| accepted

Answered

using ode45 for coupled equations

You are mixing syntax here: @(t,x)myfun Either use the syntax @(t,x)expression, where expression is the derivative or use the...

using ode45 for coupled equations

You are mixing syntax here: @(t,x)myfun Either use the syntax @(t,x)expression, where expression is the derivative or use the...

20 días ago | 0

| accepted

Answered

For several values of x, use MATLAB to confirm that sinh x = (ex − e−x)/2.

Note that sinh x and ex and e-x in your title would be sinh(x) and exp(x) and exp(-x) in MATLAB code.

For several values of x, use MATLAB to confirm that sinh x = (ex − e−x)/2.

Note that sinh x and ex and e-x in your title would be sinh(x) and exp(x) and exp(-x) in MATLAB code.

20 días ago | 0

Answered

How can I plot sin(x),sin(2x),...,sin(nx) for an input n all on the same graph? Matrix dimensions must agree error

Modifying your original code without a loop: x=[-2*pi:pi/64:2*pi]; i = input('Enter a positive whole number: \n'); z=[1:i]'; ...

How can I plot sin(x),sin(2x),...,sin(nx) for an input n all on the same graph? Matrix dimensions must agree error

Modifying your original code without a loop: x=[-2*pi:pi/64:2*pi]; i = input('Enter a positive whole number: \n'); z=[1:i]'; ...

20 días ago | 0

Answered

Dice rolling & loops

The loop could look something like this: n = numel(a); got5or6 = 0; for k=1:n if( _______ ) % you fill in the blank here...

Dice rolling & loops

The loop could look something like this: n = numel(a); got5or6 = 0; for k=1:n if( _______ ) % you fill in the blank here...

20 días ago | 0

| accepted

Answered

Matlab code for Euler method help

z(i+1) = z(i) + f*h; T(i+1) = T(i) + h; plot(T,z);

Matlab code for Euler method help

z(i+1) = z(i) + f*h; T(i+1) = T(i) + h; plot(T,z);

20 días ago | 0

Answered

how to plot euler forward method

Don't change the iteration variable withing the loop. Remove these lines from your code: n=1:11;

how to plot euler forward method

Don't change the iteration variable withing the loop. Remove these lines from your code: n=1:11;

21 días ago | 0

Answered

Creating a loop for game

The name "input" is an existing MATLAB function. It would be best to change this to something else, e.g. n. And use a differen...

Creating a loop for game

The name "input" is an existing MATLAB function. It would be best to change this to something else, e.g. n. And use a differen...

21 días ago | 0