Torsten

Last seen: Today | Active since 2013

Followers: 1 Following: 0

Statistics

Feeds

Answered
How to get the formal expression of the expansion of power of a polynomial with symbolic calculation
syms x n = 10; k = 2; a = str2sym(compose("a" + (0:n-1))) f = x^n + sum(a.*x.^(0:n-1)) g = f^k; coefficients = coeffs(coll...

2 hours ago | 0

Answered
Using MATLAB to do Numeric Calculus
Your code is correct. Here is a simpler way to define s and v. Just substitute tcrash by tmax. tcrash = 1596.2/(2*4.905)+sqrt...

14 hours ago | 0

Answered
probability that values are belong to an interval
It depends on how x is distributed. If f is the (continuous) probability density function for x, you get the probability p as ...

16 hours ago | 0

| accepted

Answered
How to use ode15s to solve a stiff ode with mass?
tf = 200; x0 = [0 0 0 0 2500 0 0.8 0 50 50 50 50 5]; M = zeros(13,13); M(1,1) = 1; M(2,2) = 1; M(3,3) = 1; M(4,4) = 1; ...

1 day ago | 0

Answered
I am getting error:Warning: Failure at t=0.000000e+00. Unable to meet integration tolerances without reducing the step size below the smallest value allowed (7.905050e-323)
If you remove the semicolon behind %% mass balances dy_dz=[1/u*(-r3) 1/u*(...

2 days ago | 0

Answered
vectorized operations on symbolic functions
rng("default") syms x1 x2 x3; % symbolic variables y = x1^3/3 + x2^2/2 - x3; % symbolic function y X = rand(500,3) % each row...

2 days ago | 0

Answered
Dot product does not work for symbolic vectors
Since you cannot assume functions to be real-valued, you must use dot(r,r) - (r_x*r_x' + r_y*r_y' + r_z*r_z' ) instead of dot...

2 days ago | 1

Answered
NaN in converting a sym to double
exp(lambda1*c0*eta0) and exp(lambda3*c0*eta0) evaluate to Infinity for s big enough (s slightly bigger than 0 already suffice...

3 days ago | 0

| accepted

Answered
How to plot a three variable complex function
You cannot plot the complete result of a function that depends on three independent variables. But you can plot it on slices (i...

3 days ago | 0

Answered
PDEPE solver failure (spatial discretization) ; partial differential/derivative temperature modelling
Please include the constants defined as globals so that we can test your code. One obvious error is the boundary condition for...

3 days ago | 0

Answered
Packed bed storage modeling probelm
rhoS = 2750; % Solid density [kg/m^3] cpS = 850; % Solid specific heat capacity [J/(kg*K)] ks = 1.28; ...

3 days ago | 0

| accepted

Answered
Extracting values from optimization variable and setting constraints
Defining vectors QTYup and QTYdown as integer vectors with 0<=QTYup(i)<=5 and -5<=QTYdown(i)<=0 (i=1,...,3 for tomatoes, onions ...

4 days ago | 0

Answered
I am trying to form a conditional loop that shows the growth of money, I put in $1000 every year and it grows by 85 for ten years but keep returning one value
money = zeros(10,1); supersaver = 1000; Return = 0.08; money(1) = supersaver*(1+Return); % Money at end of year 1 for i = 2...

4 days ago | 0

| accepted

Answered
Getting error while trying to solve boundary condition equations.
For M1 = 0-0.4, bvp4c cannot find a solution. So I started with M1 = 0.5. close all; clear all; clc; K1=0; eta1=0.2; eta2=0....

4 days ago | 0

Answered
Meijer G-function MATLAB Implementation Problem
Cbar = 1.4361*meijerG([0],[1,1.22,1.72],[0.7200,1.2200,1.4510,1.9510,1.2550,1.7550,0,0],[],331.6059)

4 days ago | 1

| accepted

Answered
Where can I learn how to do regression in Matlab
First learn how to approach your problem from the statistical side, then look how to implement this theoretical approach in MATL...

5 days ago | 0

Answered
Trying to use Newton Raphson Method to solve circuit design
w = 1; T = 2*pi/w; tn = 0:T/40:T; I = 1e-15; Vt = 25.9e-3; Vo = 10; V1 = 6; V2 = 6; R = 50e3; x0 = [0.7;0.7]; e_d = 1...

5 days ago | 1

| accepted

Answered
rror using odearguments @(T,Y)ODE_LR(T,Y,KF_L,KB_1) returns a vector of length 2, but the length of initial conditions vector is 10002. The vector returned by @(T,Y)ODE_LR(T,Y
Use Non_Active_Receptor_concentration0 = 100; Active_Receptor_concentration0 = 0; instead of Non_Active_Receptor_concentrat...

5 days ago | 0

Answered
plot the function which is dependent on x, y and z with x, y and z on three axis.
You can plot the function on slices (i.e. 2d-objects (e.g. planes)) through the volume of interest. Of course a full plot over ...

5 days ago | 0

Answered
How do I feed ODE45 a cell array ?
Use state_dot as the function that you pass to the integrator. I followed your coding - the problem setup could have been much ...

6 days ago | 0

| accepted

Answered
Left and Right side have different elements
Replace t by t(i) in the loop.

6 days ago | 0

Answered
Getting integrand error from triple integral code code
f must return a vector of ones of the same size as the inputs x, y or z (which all have the same sizes). Thus replace f = @(x,...

6 days ago | 1

Answered
need to compile outputs from for loop into array
In most cases, outputarray = zeros(1,nnz(mod(yeararray,4)==0)) will work.

7 days ago | 0

| accepted

Answered
Error ---> Brace indexing is not supported for variables of this type.
clc clear g = 9.81; L = 0.1; m = 0.5; zeta = 2; x0 = 0.25; v0 = 0; tspan = [0 10]; step_sizes = [0.1, 0.05, 0.025...

7 days ago | 0

Answered
error lsqcurvefit for integrate function
Replace % Define the model function (vectorized for efficiency) function y = model(theta, x) k0 = theta(1); E0 = theta(2...

7 days ago | 0

Answered
How to set Drichlet and Neumann boundary conditions for a set of ode in BVP solver
Is there a possibility that I can set three boundaries (mix between Drichlet and Neumann boundary) for a set of two ODE's? No....

7 days ago | 0

Answered
Trying to determine roots of this polynomial (Det)
k = [1750 -750 0; -750 1250 -500; 0 -500 500]; m = [75 0 0; 0 75 0; 0 0 50]; eig(k,m)

7 days ago | 1

Answered
Civil engineering coding problem ( I dont know where the problem in my code is and how to solve it, thanks in advance )
"coord" has only three rows, but you try to access coord(4,1) and coord(4,2) which do not exist (see above).

7 days ago | 0

Answered
Why I am getting this error message? How to fix it?
Maybe you mean this function: function v = vec(m) %Helper function to turn a matrix of any size into a column vector using (:)...

7 days ago | 1

Answered
How to draw a three-dimensional function image corresponding to the analytic expression of a given piecewise function image
f = @(x,y)(4*y-4*x.*y-2*y.^2).*(y<x).*(x+y<1)+2*(1-x).^2.*(y<x).*(x+y>=1)+(4*y-4*y.^2-2*x.^2).*(y>=x).*(x+y<1)+(2-4*x+4*x.*y-2*y...

7 days ago | 0

| accepted

Torsten

Statistics

MATLAB Answers

MATLAB Answers

File Exchange

MATLAB Answers

MATLAB Answers

MATLAB Answers

MATLAB Answers

MATLAB Answers

Feeds