Torsten
Statistics
0 Questions
5,145 Answers
RANK
19
of 277,209
REPUTATION
12,492
CONTRIBUTIONS
0 Questions
5,145 Answers
ANSWER ACCEPTANCE
0.00%
VOTES RECEIVED
1,392
RANK
of 18,645
REPUTATION
N/A
AVERAGE RATING
0.00
CONTRIBUTIONS
0 Files
DOWNLOADS
0
ALL TIME DOWNLOADS
0
RANK
of 126,789
CONTRIBUTIONS
0 Problems
0 Solutions
SCORE
0
NUMBER OF BADGES
0
CONTRIBUTIONS
0 Posts
CONTRIBUTIONS
0 Public Channels
AVERAGE RATING
CONTRIBUTIONS
0 Highlights
AVERAGE NO. OF LIKES
Content Feed
Why I got an error in "plot(X1,sqrt(2*g*X1)*tanh(sqrt(2*g*X1)/(2*L)*t),'o')" in the following code??
Seems you got the derivative wrong. % Constants g = 9.81; % m/s^2 v = 4; % m/s t = 3.45; % s L = 5; % m X0 = 3; % m Es = ...
9 hours ago | 0
Dsolve solves differential equation with a very long output, is something wrong?
syms lambda mu real positive syms t p1(t) p2(t) p3(t) eqn1 = diff(p1,t) == -5*lambda*p1 + mu*p2; eqn2 = diff(p2,t) == 5*lambd...
9 hours ago | 1
| accepted
Extracting an intermediate vector from ode45
Use s = zeros(numel(t),1); u = s; z = s; for i = 1:numel(t) [~, s(i),u(i),z(i)] = smc_obsv(t(i), x(i,:)); end instead o...
14 hours ago | 0
How to numerically solve system of equations and differential equations simultaneously?
MATLAB's ode solvers allow a mixture of differential and algebraic equations as in your case. These systems are called differen...
17 hours ago | 2
Plot a quadrilateral having 8 points (each side passes through two points)
P1 = [0 0 0]; P2 = [1 0 0]; P3 = [0 1 0]; P4 = [0 0 1]; P5 = [1 1 0]; P6 = [0 1 1]; P7 = [1 0 1]; P8 = [1 1 1]; hold on ...
1 day ago | 0
ode45 taking too long to solve set of ODEs.
Try to use ode15s instead of ode45. And prescribe the vector of output times. Otherwise ode45 will save all of its successful i...
1 day ago | 0
How to get a better fit with given equations that would approximately match given data output?
Use rng("default") as first line in your code.
2 days ago | 0
Finding a nonlinear equation root
fun = @(x)cosh(x).*cos(x)+1; icount = 0; for i = 1:2:21 x0 = i*pi/2; icount = icount + 1; xroot(icount) = fzero(fun,x...
2 days ago | 0
Problen with optimization of two parameters, and with differential equations
All curves that stem from your model tend to sqrt(keq/(1+keq)). Since your measurement data tend to 0.5, my guess is that you ha...
2 days ago | 0
Problen with optimization of two parameters, and with differential equations
The parameters you want to fit are not independent. Since T remains constant, the complete expression k0*exp(-Ea/(8310*T))*cc*...
2 days ago | 0
How to use fmincon's constraintviolation option
Can you tell me what I'm doing wrong? You try to betray yourself that your problem has a solution. But it has none.
3 days ago | 1
| accepted
Solving Integral in matlab
syms s a=sym('1'); b=sym('71/100'); c=sym('53/100'); f = 1/((s+a^2)*(s+a)*(s+b)*(s+c)); F = a*b*c/2*int(f) limit(F-subs(F,...
3 days ago | 0
im looking for the error in the following code, please help
DalphaS = CaputoDerivative(S(1:n),alpha,dt); DalphaI = CaputoDerivative(I(1:n),alpha,dt); DalphaR = CaputoDerivative(R(1:n),al...
3 days ago | 0
Systems of Trigonometric Equation does not return any value
You get contradicting results for theta10. Thus your system is not consistent. format long AB = 105.58445229; CD = 16.557818...
4 days ago | 0
I couldn't match the index no., can anyone please help me to resolve this issue.
yita_mn has size 6x6, but your loop in which you access yita(i,j) runs for i = 1:16 and j = 1:16. And dy should be allocated as...
4 days ago | 0
Drawing a heatmap of the max value of a state variable against two parameters in the system
tspan = 0:0.01:100; y0 = [100000;0;0;1;0]; p = 0.001:0.001:0.01; eta = 0.001:0.001:0.02; for i = 1:numel(p) for j = 1:num...
4 days ago | 1
| accepted
volume of a packed bed reactor
One main error in your code is that you use variables that you define later on. Programming works sequentially: If you define a...
4 days ago | 0
| accepted
Finding random normally distributed number above zero
Maybe taking abs(randn()) will give you something useful. It's the folded normal distribution (which in this case equals the tru...
4 days ago | 0
Need to remake this optimize function to be efficient takes way to long currently
You can easily solve when det(m) = 0 in advance and insert the formula in your code. This will enhance speed enormously: syms p...
5 days ago | 1
ode45 how to write differential equation (within function) that incorporates itself at previous time stamp
This is a delay differential equation. Use dde23 to solve. What are the variables you solve for in the equation dIdt = V*S+ep...
5 days ago | 0
calculate integral of this function
F1 = double(int(dM_dH^2/cos(theta),x,0,Pi)) instead of F1= integral(@(x) (dM_dH).^2/cos(theta),0,Pi)
5 days ago | 0
| accepted
Seventh order differential equation
% Set model parameters l = 1; P = 1; Ga = 1; Eatilde = 1; ha = 1; E1tilde = 1; h1 = 1; E2tilde = 1; h2 = 1; xmesh = ...
6 days ago | 1
| accepted
Ode45 calling a matrix and an array in a function
Your arguments to q_dot are inverted: Use function q_dot = q_dotf(ts,q) instead of function q_dot = q_dotf(q,ts) And note t...
6 days ago | 0
| accepted
Stability analysis of a non-linear ODE system
syms Sci C Sr Sh R Cf Cp Ce E HR H Sp P k1 k2 k3 k4 k5 k6 k7 k8 k9 k10 k11 k12 k13 k14 k15 k16 p1 p2 p3 mu eta alpha theta CL F...
7 days ago | 1
How to solve the parameters of a diffusion model?
Make a code to determine the roots of your second equation. Make a code that evaluates the infinite sum to determine q_t from y...
7 days ago | 0
The error is " Unrecognized function or variable 'Vmp' " on line 21, but i need to find the root of Vmp(which is -0.05). Can anyone help me with this??
A fixed point iteration for your function does not converge. Didn't you see this in your former question ? And if it worked, g ...
7 days ago | 0
Summing scattered data over a 2D grid
format long % Generate random coordinates and concentrations n = 100000; x = -0.5+rand(n,1); y = -0.5+rand(n,1); conc = 10*...
8 days ago | 0
I need the maximum power to be -0.05 but my output is -3709.4691. Does anyone have any tips?
% Given values Voc = 0.5; % V T = 297; % K q = 1.6022e-19; % C kB = 1.3806e-23; % J/K a=q/(kB*T); % Define the function f...
8 days ago | 0
solve the mass spring system where the mass matrix depends explicitly on time
Setting y1' = y3 and y2' = y4, you arrive at the following code: M = @(t) [t 0; 0 t]; K = [2 1;5 8]; MM = @(t)[eye(2),zeros(2...
8 days ago | 0
Write a matlab function for approximating the value of pi (certain equation)
You sum twice. And (2*n-1)^2*(2*n+1)^2 = (4*n^2-1)^2 Either use n = 20; p = 0; for k = 1:n p = p + 1/(4*k^2-1)^2; e...
8 days ago | 1
| accepted