Risposto
Hello, I am trying to solve the ODE but I have an error I do not understand how to solve. Would appreciate some help. Thank you!
Like this? %% ODE45_estimation % ------ Constants ------ % m = 1248.5; g = 9.81; W = m*g; S = 17.1; % ------ Time int...

quasi 4 anni fa | 0

Risposto
plotting Solid of revolution of a function
How about: X = 0:0.1:pi; R = sin(X); [z,x,y] = cylinder(R); surf(x,y,z)

quasi 4 anni fa | 2

Risposto
Split an array using specific points
Like this? Y = [5;2;3;6;7;9;5;5;6;3;8;7;74;9;63;47;10;3]; diiv = [5; 9; 15]; y1 = Y(1:diiv(1)); y2 = Y(diiv(1)+1:diiv(2)); ...

quasi 4 anni fa | 1

Risposto
Algorithm for Fractional power calculation
How about using the Newton-Raphson algorithm. Here's the basic idea: % x^n = b % Let f(x) = x^n - b % dfdx(x) = n*x^(n-1) %...

quasi 4 anni fa | 0

Risposto
Stochastic differential equation Gompertz plotting a graph
You have x(1) = 0, so log(x(1)) is -Inf, which means that subsequent values of x will be NaNs.

quasi 4 anni fa | 0

| accettato

Risposto
Defining integration function in matlanb
Like this? % If Cp is a constant then the integral is just Cp*T % If it is a function of T, then define the function: Cp = @(...

quasi 4 anni fa | 0

Risposto
Plotting and coding problems
More like this perhaps: T = 300:0.5:500; w = 0.224; Tc = 304.2; ...

quasi 4 anni fa | 0

Risposto
I am getting error, "array indices must be positive or logical integer values". Please help
You need multiplication signs between n and (w1+k) in calculations of u and z: [X,Y] = meshgrid(-4:0.2:4); n=2; k=1; %n=1,...

quasi 4 anni fa | 0

Risposto
Need help in finding all the solutions of non linear ellipse equations using Newtons Method using following code.
Check your Jacobian equations, especially the term: 2*x(1)+2*x(1)+4; I think this should be 2*x(1)+2*x(2)+4;

quasi 4 anni fa | 1

Risposto
Displayin command just once for the n number of iterations
What about if k == 0 disp("Hello World!") elseif k == 480 disp("Bye Bye World!") end

quasi 4 anni fa | 0

Risposto
Index in position 1 is invalid. Array indices must be positive integers or logical values.
In line V(i,j+1) = ((1-2*s-k).*V (i,j))+ (s.*(V (i-1,j) + V (i+1,j))) - (k.*(Rm.*Iinj+Vrest)); you have V (i-1,j). When i is ...

quasi 4 anni fa | 0

| accettato

Risposto
What is the way of solving the noise added differential equations in matlab?
Try to avoid the (mis)use of global! %% ----------------INPUT PARAMETERS--------------------------- f=0.33; alpha=0.5; w0=sqrt...

quasi 4 anni fa | 0

| accettato

Risposto
Unable to perform assignment because the left and right sides have a different number of elements.
Like this? clear %Vmax values for each reaction step v1=0.01:0.01:2; v2=1; %Km values for each reaction step K1=1; K2=1; ...

quasi 4 anni fa | 0

| accettato

Risposto
How do I create a stacked plot like attached figure using subplot command ?
Try help stackedplot

quasi 4 anni fa | 0

| accettato

Risposto
Second order equations. Cauchy task
More like this: % Write equations as follows: % dudx = v u(0)=-1 % dvdx = -v/4-6u v(0)=0 % % Then you hav...

quasi 4 anni fa | 1

| accettato

Risposto
Runge Kutta 3 ODE
You can't have t(i)=t(0)+i*h; indices must be a positive integer - you have t(0), Matlab doesn't like this!

circa 4 anni fa | 0

| accettato

Risposto
Why it is saying this? Index in position 2 is invalid. Array indices must be positive integers or logical values. Error in Homework (line 21) surf(x,y,F(i,j))
Like this x=[-5:5]; y=[-5:5]; for i=1:numel(x) for j=1:numel(y) if (x(i)>=0)&&(y(j)>=0) F(i,j)=2.*x(i)+y(j)^2; en...

circa 4 anni fa | 0

| accettato

Risposto
Definite integral with parameter
You could try something like this: a = 0.043478; q = 100; % Ha W = 40; t = 4; c = @(x,Fa )(5.*x.^2-Fa.*x)./(((a.*x+q).*W-((...

circa 4 anni fa | 0

| accettato

Risposto
Integration of a joint pdf with for loop?
When k = 1and i = 1 then k-i = 0, therefore the second index in fun(i,k-i) is 0, but indices must be positive integers in Matlab...

circa 4 anni fa | 0

Risposto
How would you I write a matlab function for theta that represents the 2nd order differential of the torque equation?
Does this help? % Theta seems to have two components, so: D(theta) = [p1+p2+2*p3*cos(theta(2)) p2+p3*cos(theta(2)); ...

circa 4 anni fa | 1

| accettato

Risposto
Error message with arrays
You only need h = [2:20:220]; t_a = 293; q_flux = 800; n = length(h); y = q_flux./h+t_a; i = i+1; plot(h,y)...

circa 4 anni fa | 0

Risposto
How to solve two coupled differential equations using ode45.
"Why is my previous code is not working properly by treating the problem to be a DAE." I don't know! "Also I noticed you wrote...

circa 4 anni fa | 1

Risposto
Multiple Anonym Function return wrong values
Here's one solution: C = [4622 5480 18923]*1e-12; U = [454e3 454e3-20e3 20e3].*sqrt(2); W = Energy(C,U); disp(W) function W...

circa 4 anni fa | 0

| accettato

Risposto
Hi, is it possible to have variables in an ode45 that varies according to another function?
Put something like [T, p, rho, speedsound] = atmos(height); immediately after height = x(1); in function dynameqn. Not sur...

circa 4 anni fa | 0

Risposto
root finding and plot of graph
Are you looking for something like this? hi = 1e-6; lo = 1e-9; n = 500; dt = (hi-lo)/n; t3 = lo:dt:hi; y = f(t3); subplo...

circa 4 anni fa | 0

Risposto
Calculator result and MATLAB result are different
Have you checked that both are using degrees, or both radians? The Matlab as you've written it is using radians (use sind(24) a...

circa 4 anni fa | 1

Risposto
How can I rotate a set of points by an angle around its origin?
Do you mean like this? ITM = [461400.0758 617474.3059 461400.8069 617475.4609 461400.8069 617475.4609 461402.3069 617477.633...

circa 4 anni fa | 1

| accettato

Risposto
Unable to perform assignment because the indices on the left side are not compatible with the size of the right side.
Try replacing x(4) = Velocities(3608.9,5249.3); by x(4,:) = Velocities(3608.9,5249.3); Same for x(5).

circa 4 anni fa | 0

Risposto
Removing zeros from matrix
Like this? A = [1 2 3 4 5 6; 1 2 3 4 5 6; 1 2 3 4 5 0; 1 2 3 4 0 0; 1 2 3 0 0 0]; B = A'; B=B(:); B(B==0)=[]

circa 4 anni fa | 0

Risposto
Error in ODE arguments
I suspect you mean like this (notice the way M divides, using the back-slash): % Numerical solution of IVP % M*xddot + C*x...

circa 4 anni fa | 2

| accettato