# Torsten

**Professional Interests:**

**Statistics**

**MATLAB Answers**

**0** Questions

**3,394** Answers

RANK**24**

of 257,635

REPUTATION**8,084**

CONTRIBUTIONS

**0** Questions

**3,394** Answers

ANSWER ACCEPTANCE **
0.00%
**

VOTES RECEIVED**941**

RANK

of 17,752

REPUTATION**N/A**

AVERAGE RATING**0.00**

CONTRIBUTIONS**0** Files

DOWNLOADS **0**

ALL TIME DOWNLOADS**0**

RANK

of 109,912

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**

Ho to perform Monte Carlo Simulation with a factor and a vector.

data_TSS=readtable('turb_danuube_15min.csv'); % data TSS & time TSS=data_TSS(1:1:353665,2); time=data_TSS(1:1:353665,1); %% T...

15 hours ago | 0

A simple error minimization (optimization) problem

A = [C1,C2,C3]; b = L1+L2+L3; C = A\b; C1 = C(1) C2 = C(2) C3 = C(3)

17 hours ago | 0

| accepted

while loop to count partail sum of series

sum_{i=1}^{i=N} 1/(i*(i+1)) = 1 - 1/(N+1) No need for such a difficult while construction with precision estimate.

21 hours ago | 0

How to get energy using power summation equation in MATLAB?

t = [0 t]; E = sum(P.*diff(t))

21 hours ago | 0

| accepted

Solve system of differential equations

x01=-5;y01=-7; x0=[x01,y01]'; beta=25; syms x1 y1 t x11=[x1,y1]'; c1=(1/2)*(x1-i*sin(t))^2+(3/2)*(y1-i*cos(t))^2; row=100*...

21 hours ago | 0

Loop to check if 2 conditions are true for values 1:n, (1:n)+1, (1:n)+2, etc. in a vector and save ID

for i=1:365-7+1 period = rainfall(i:i+6); sum_period = sum(period); n = numel(find(period>0)); if sum_period...

22 hours ago | 0

| accepted

Use ode45 when I have a constant that varies over time.

n=10; xmin_in = 15; % Minimum inlet temperature xmax_in = 16; % Maximum inlet temperature T_in = xmin_in+ra...

22 hours ago | 0

| accepted

Determine y for a function with a symbolic variable x

syms x F(x) = x.^2; F(0)

24 hours ago | 0

Rotate the coordinate system to align an existing plane with Y'Z' plane

https://stackoverflow.com/questions/9423621/3d-rotations-of-a-plane But are you sure that the parabola in the YZ plane is not a...

1 day ago | 0

How to solve a system of differential equations using ode45 and vary the magnetic field and plot the solutions as a function of this field?

b=[0:0.1:10]; d0=0.2 wb=1 wd=wb-d0 wl=wb db=0.18 dd=0.05 a=0.02 Ua=0.15 Ub=0.02 g1= 0.001 g2=g1 g3= 0.0001 g4=g3 %...

2 days ago | 0

Problems with declaring the variable right

xg= 7*2*5*1*6*3; xh= 3/10; F1=3; F2=100; T= 1/500; fun = @(k1,i) xg*sin(2*pi*F1*(k1-i)*T); fun(14,0:14)

2 days ago | 0

Spectral factorization of polynomial x(s)

Let a(s)*a(-s)+b(s)*b(-s) = p0 + p1*s + ... + p2n*s^(2*n) Determine the roots of a(s)*a(-s) + b(s)*b(-s) The 2*n roots ap...

2 days ago | 0

| accepted

Correctly evaluate exponential function for limit values

function result = f(x,g,h,N) result = N^2*ones(size(x)); K = g(x).*h(x) ~= 0; if (any (K(:))) result(K) = your exp...

2 days ago | 0

"Not enough input arguments, Error in my_ode_without_tld_1 (line 24) F = interp1(tdata,ydata)" why this error is showing?

If F has the same number of elements as i, the following should work: i = 0.1:0.1:30; y0 = [0 0]; F = load('shear force.dat')...

3 days ago | 1

How to find coefficient matrix?

I assume the last xi^n should be xi^m in sum_{i=1}^{n} y_i*x_i^n x = [1 3 5 7 9]; y = [3 -1 3 78 -0.2]; m = 4; xx = x.^(...

3 days ago | 0

| accepted

How to use Monte Carlo simulation to calculate the area of polygon?

Randomly create points inside the bounding box. Decide which of the points are inside (P_i) and which of the points are outside...

3 days ago | 0

FMINCON requires the following inputs to be of data type double

Although I don't understand your problem setup because the variables in the objective function are not influenced by your constr...

3 days ago | 0

Error using plot3 not enough input arguments

plot3(i_1, j_1, k_1) plot3(i_2, j_2, k_2) plot3(i_3, j_3, k_3) plot3(i_4, j_4, k_4) instead of plot3(i_1, j_1, k_1, tspan) ...

3 days ago | 0

| accepted

Problem with PDE resolution

function [c,f,s]=heatpde(x,t,u,dudx) c = [1 ;1]; f = [0.1*dudx(1);0.5*dudx(2)]; a = 0.2*u(1)*u(2)/(1+0.1*u(1)+0.02*u(2)+0.03*...

4 days ago | 0

Spatio-temporal plug flow reactor simulation with method of lines

function PFR_t_L_sim_v2_1 % function method_of_lines clear all; clc; close all function dCdt = reaction(t,C) C1 ...

4 days ago | 1

| accepted

Maximum Likelihood Estimation function

x = your data vector fun = @(a) n/a + sum(x) - 2*a*sum(exp(a*x)./(exp(a*x)+1); a = fzero(fun,1)

5 days ago | 0

| accepted

Hello, the code to enter and raise the integral does not run, what is wrongly written?

a=0.286; b=0.1851; c=1.272; f=@(t)t.^(0.5*c).*exp(-0.5*a*t); I = sqrt(b)*integral(f,0,Inf)

5 days ago | 0

Solve system of equations without Symbolic Math Toolbox for Compiler

Did you try "solve" on your system of equations with the numerical vlist coeffcients replaced also by symbolic variables ? syms...

5 days ago | 1

| accepted

Error fit: Complex value computed by model function

Fit under the constraint b^2-a^2 >= 0

5 days ago | 0

pde mass transfer with reaction

clc;clear all;close all; L = 1; x = linspace(0,L,50); t = linspace(0,1,50); % t = [linspace(0,0.05,20), linspace(0.5,5,10)];...

5 days ago | 0

Crank-Nicolson scheme for system of differential equations

y0(1,1) = 0.2; y0(2,1) = 0; x = 0:0.125:100; h = x(2)-x(1); y = zeros(2,numel(x)); y(:,1) = y0; for i4 =...

5 days ago | 0

| accepted

Need help plotting this

n = 36; T=20*10^-3; t=0:T/1000:T; w=2*pi*1/T; %f=i(t) f=@(t)0+100.*(t>1*T/12).*(t<5*T/12)-100.*(t>7*T/12).*(t<11*T/12); a0...

6 days ago | 1

| accepted

Help in looping and iterating constant values in differential equation which is being solved my numerical method.

Define all time-dependent parameters as function handles and apply these handles in your function definition. As an example for...

6 days ago | 0

Using Implicit Euler Method with Newton-Raphson method

f = @(t,y) -20*t*y^2; f_y = @(t,y) -40*t*y; t0 = 0; T = 1; y0 = 1; h = 0.01; tol = 1e-8; N = 100; y = imp_euler(f,f_y,t0...

6 days ago | 0

| accepted

How do I create a loop for this calculation in MATLAB?

%constants c = 1.1e-10; m = 3.3; p500 = [ -0.3142 2.3172 -7.0249 11.2867]; p2000 = [ -1.2577 9.2728 -28.1056 45...

6 days ago | 0

| accepted