Risposto
Projectile motion without drag
You don't define the Vx(k) value before using it, hence the error. Since your Vx values don't change you can fix this error by ...

circa 5 anni fa | 1

| accettato

Risposto
Define x,y,z and evaluate them.
You can define variables simply by typical assignment. E.g., x=5; y=4; etc. You can use * for multiply and ^ for exponentiat...

circa 5 anni fa | 0

Risposto
Fortran MEX returns only some elements of output array.
The problem isn't using allocatable variables ... MATLAB mex functions can handle that just fine. The problem is that you are u...

circa 5 anni fa | 1

| accettato

Risposto
System of ODEs with interdependent variables
I would advise writing a derivative function for this. E.g., dydt = myderivative(t,y,gna,gk,gl,Vna,Vk,Vl,I) V = y(1); n = y(...

circa 5 anni fa | 0

| accettato

Risposto
Outlet Limit: water depth and flow rate
You are being asked to write a function that takes H and Q as inputs, and outputs the diameter D. So the function outline would...

circa 5 anni fa | 0

Risposto
writing out a quadratic
Use the element-wise operators that include the "dot" (such as .^ instead of just ^), and use the multiply operator * between th...

circa 5 anni fa | 0

Risposto
How can I select three random sample of a matrix with minimum separation between samples?
A rejection method is straightforward. Simply make your selections in a forever loop until you get something that works. The a...

circa 5 anni fa | 1

| accettato

Risposto
Cell to matrix under new variable
Do not do this! You are taking a method of storing your data that is easy to maintain and access (e.g., in a loop) and proposin...

circa 5 anni fa | 0

Risposto
Adding 2 matrices together in matlab, got errors instead
result = [a,b];

circa 5 anni fa | 1

| accettato

Risposto
Reduce size of multidimensional array
One way using summing: result = squeeze(sum(your_array,3));

circa 5 anni fa | 2

| accettato

Risposto
Euler's code for multiple ODEs
You have two states, so your result should have two values per step. So this code: y = zeros (1,N); % Initialize the Y vec...

circa 5 anni fa | 1

Risposto
Runge-Kutta 4th Order on a 2nd Order ODE
All of your derivatives should be with respect to the independent variable x. In the context of your problem, dy/dz does not ma...

circa 5 anni fa | 1

| accettato

Risposto
Something wrong with quat2angle() function?
I agree that this looks like a bug, and you should submit a bug report to TMW for this. I ran it in R2017a and R2020a and both ...

circa 5 anni fa | 0

| accettato

Risposto
want to solve (d^2 y)/(dx^2 )+dy/dx-6y=0 using 4th order Runge-Kutta method with y(0) = 3 and y’(0) = 1
Change your function handle from this f = @(x) (dy^2)/(dx^2 )+dy/dx-6*y; to this f = @(x,y) [y(2);-y(2)+6*y(1)]; The last pa...

circa 5 anni fa | 1

| accettato

Risposto
How to use Bisection Program to locate Roots of given function
If the derivative depended only on the time t then your strategy for finding the maximum/minimum points would be a good one. Ju...

circa 5 anni fa | 0

Risposto
ODE45, ODEfunc question
What you might want to do is pass a function handle for this. E.g., in your calling code create the function handle: C1 = what...

circa 5 anni fa | 1

| accettato

Risposto
Runge-kutta 4th order Method Help
All of your F's need to be functions that operate on the current full states. E.g., take this one F_C1 = (xc*yc*((C1.*R1)./(R1...

circa 5 anni fa | 0

| accettato

Risposto
Difference/usage of mwIndex, mwSignedIndex and mwSize in mex C-code
MATLAB provides multiple libraries that mex code links to. There can be both 32-bit and 64-bit versions of functions in the var...

circa 5 anni fa | 1

| accettato

Risposto
Matrix index is out of range for deletion.
You should step through this with the debugger to see what is going on. I.e., click to the left of a line and then run your cod...

circa 5 anni fa | 0

Risposto
How to stack k copies of a matrix
B = repmat(A,1,1,k);

circa 5 anni fa | 0

Risposto
how can i find the co-ordinates of intersections of two lines?
Suppose you have two line equations: 2x + 3y = 1 3x - 5y = 11 Then you can form the matrix A and the right hand side vector b...

circa 5 anni fa | 0

Risposto
difference between euler angles and rotation vector
The rotvecd( ) function gives a scaled rotation axis. I.e., what it would take to do a single rotation of an object on all thre...

circa 5 anni fa | 2

| accettato

Risposto
Cannot input function handles into array
When you do this: basisN0 = @(x)[ 1 ]; basisN1 = @(x)[ 1, x ]; basisN2 = @(x)[ 1, x, x^2 ]; basisN3 = @(x)[ 1, x, x^2, x^3 ]...

circa 5 anni fa | 0

| accettato

Risposto
Analytically finding the roots of polynomials with symbolic variables
Tell the solve( ) function that the max degree of the polynomial is 3 to force explicit solutions for the result: syms a len p...

circa 5 anni fa | 0

| accettato

Risposto
Matlab unable to find solution to cubic polynomial
Tell the solve( ) function that the max degree of the polynomial is 3 to force explicit solutions for the result: syms a p = (...

circa 5 anni fa | 0

| accettato

Risposto
Replacing a for loop with matrix multiplication
result = sum(x.*(inv(SIG)*x)); or result = sum(x.*(SIG\x));

circa 5 anni fa | 0

| accettato

Risposto
How to find the numerical intersection points for a sine and linear function to find solutions ?
You can't do floating point calculations like this and expect the results to match exactly. You would need to use a tolerance t...

circa 5 anni fa | 0

Risposto
Solving a Second Order Differential Equation
Here is what I get: >> syms y(t) >> eqn = 4*diff(y,t,2)+32*diff(y,t)+24*y^3 == 192; >> Dy = diff(y,t); >> cond = [y(0)==1.5,...

circa 5 anni fa | 0

Risposto
Runge Kutta4 for IVP : y' = ty
Good start! Some corrections noted below. This line y0 = exp(1); should be this according to your text y0 = 3; And this li...

circa 5 anni fa | 0

Risposto
Power Method Eigenvectors Code
The only obvious thing I would note is that the following line doesn't accomplish anything because you don't use err in your cod...

circa 5 anni fa | 1