# James Tursa

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

**Statistics**

RANK**17**

of 243,877

REPUTATION**12,396**

CONTRIBUTIONS

**17** Questions

**3,973** Answers

ANSWER ACCEPTANCE **
41.18%
**

VOTES RECEIVED**2,224**

RANK**165** of 17,057

REPUTATION**7,174**

AVERAGE RATING**4.60**

CONTRIBUTIONS**32** Files

DOWNLOADS **144**

ALL TIME DOWNLOADS**67163**

CONTRIBUTIONS**0 Posts**

CONTRIBUTIONS**0** Public Channels

AVERAGE RATING

CONTRIBUTIONS**0** Highlights

AVERAGE NO. OF LIKES

**Content Feed**

help with my 4th order runge kutta code

My guess without seeing the rest of your code is that phi is a column vector, so when you add phi*h to the row vector Y(i,:) it ...

11 hours ago | 0

| accepted

how to change from for loop to while loop

You never change i or Errors within the while loop, so how is the loop ever supposed to exit? Maybe you want to add this line: ...

11 hours ago | 0

Is mxCopyPtrToInteger4() not compatible with Interleaved Complex API ?

See this post: https://www.mathworks.com/matlabcentral/answers/93860-why-does-the-fortran-mex-api-function-mxcopyptrtointeger4-...

10 days ago | 0

| accepted

What should go in a next-generation MATLAB X?

Symmetric variables and Hermitian variables. MATLAB could implement bit flags in the mxArray header to indicate this and they c...

12 days ago | 0

Why indexing vs not indexing the 2D matrix will lead to different result?

Some links related to this known behavior: https://www.mathworks.com/matlabcentral/answers/599818-computational-complexity-of-m...

12 days ago | 0

Properly resizing mxCreateStructMatrix post-2018a?

TMW doesn't publish the complete rules for how the data areas of structs are organized, but historically they have always put ea...

19 days ago | 0

How to get the ceofficients of a expression derived by ODE or integral ? For exampli,C1 and C2 in u

Typically these would be obtained by combining the solution with initial values for u. Do you have initial values?

2 months ago | 0

Is it possible to vectorize this 'for loop' involving multiple matrix inversions

You might look at this FEX submission by Bruno Luong: https://www.mathworks.com/matlabcentral/fileexchange/24260-multiple-same-...

2 months ago | 0

| accepted

Mulitplying 3D matrix by 3D matrix - result being 4D.

result = pagemtimes(matrix1,reshape(matrix2,3,512,1,512));

3 months ago | 0

| accepted

Matrix dimensions must agree error problem

Use element-wise division with the dot: Y=acosd(sin(altitude)*sin(latitude)-sin(declination)./(cos(altitude)*cos(latitude))); ...

3 months ago | 0

Combining two Matrices every other row

You could do direct assignment. E.g., [m,n] = size(A); C = zeros(2*m,n); C(1:2:end,:) = A; C(2:2:end,:) = B;

3 months ago | 0

Gauss-Elimination method

See the rank( ) function: https://www.mathworks.com/help/matlab/ref/rank.html?s_tid=doc_ta

3 months ago | 1

Looping num2str

x(i) references a single element of x. However, the expression [num2str(y(i)),'%'] generates multiple characters. In essence, ...

3 months ago | 1

Mean computaion in matlab

You can specify the dimension to use with the mean( ) function. E.g., maybe this is the computation you want: data = your 3D m...

3 months ago | 0

What is the mean?

B' is the complex conjugate transpose of B. If B is real, then it is equivalent to just the tranpose and the assignment simply ...

3 months ago | 0

| accepted

How can I model equation of motion (2nd order ODE) when system matrices are in terms of state variables?

Yes you can have derivative functions that depend on state variables. This is quite common. E.g., define a 4x1 state vector y th...

3 months ago | 0

empty double matrix in case of nul value

Do you mean you want to do something like this: if( isempty(malicious) ) malicious = 0; end

3 months ago | 0

To allocate character array='2a'

'2a' is a two element character string. You can't assign two elements to a single element. I.e., you can't assign a 1x2 vector...

3 months ago | 0

How to get the desired quaternion representation from a rotation matrix?

As you already know, both q and -q represent the same rotation. Which sign to pick is entirely up to the underlying algorithm. ...

4 months ago | 0

| accepted

Implementation of Runge Kutta Numerical Solution for system of ODE's

The main technical problem (other than your syntax not being MATLAB) is that you are trying to solve your four 1st order equatio...

4 months ago | 0

| accepted

How to convert complex float to complex integer in MEX gateway function?

I would guess you can just use the appropriate data types. E.g., mxComplexInt32* Data1 = mxGetComplexInt32s(prhs[0]); mxComp...

4 months ago | 0

| accepted

Can you help me with this loop?

Can you use num2str? Message = [num2str(TotalMinutes) ' minutes are equal to ' num2str(Hours) ' hours and ' num2str(Minutes) ' ...

4 months ago | 0

How to separate a string(1x1 cell) into a 1x4 cell

If the hex codes are always 8 characters, why can't you just pick off the characters you want? E.g., {result(1:8),result(10:17...

4 months ago | 1

how to generate a code for solving 2st order differential equations using improved runge kutta 4th order method

The posted code actually has a bug. Also it is almost set up properly for vector solutions, but not quite. Make these changes:...

4 months ago | 0

Computer crashed when calculating matrix (mixed with sparse and full) multiplication and summation

Remember that each element of the B*C result is simply the dot product of a row of B with a column of C. If that column of C ha...

4 months ago | 0

| accepted

Fast vector reshaping/permutation

Don't do the permute( ) operation. Just use pagemtimes( ) downstream in your code with the appropriate 'transpose' option. Thi...

4 months ago | 1

Can someone explain this Loop?

Of the numbers 1 ... 10 that get produced by the randperm( ) call, only three of them are less than 4: 1, 2, 3. 1 will cause n...

4 months ago | 0

| accepted

How to generate unique random integers between 1 to n (with no possibility of Sequence)

If 1:n is the only sequence you don't want, just use randperm(n) and reject the one case you don't want. I.e., if you get it the...

4 months ago | 0

How do solve four variables with one equation using Runge-kutta?

You have three state variables, namely v, lc, and z. But you only have code for updating v. You need code for updating lc and ...

4 months ago | 2

| accepted