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)

9 mesi fa | 1

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)); ...

9 mesi 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) %...

9 mesi 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.

9 mesi 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 = @(...

9 mesi fa | 0

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

9 mesi fa | 0

Risposto
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,...

9 mesi 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;

9 mesi 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

10 mesi 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 ...

10 mesi 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...

10 mesi 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; ...

10 mesi fa | 0

| accettato

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

10 mesi fa | 0

| accettato

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

10 mesi 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!

10 mesi 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...

10 mesi 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-((...

10 mesi 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...

10 mesi 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)); ...

10 mesi 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)...

10 mesi 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...

10 mesi 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...

10 mesi 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...

11 mesi 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...

11 mesi 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...

11 mesi 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...

11 mesi 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).

11 mesi 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)=[]

12 mesi 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...

12 mesi fa | 2

| accettato

Risposto
How to Approximate The Solution for an Initial Value Problem?
Try f = @(t,x) [x(2); sin(1.9*t)-4*x(1)]; [t,x] = ode45(f,[0,150],[1,0]); plot(t,x)

12 mesi fa | 1

| accettato