Community Profile # Davide Masiello

### Paul Scherrer Institut (PSI)

Last seen: Today Active since 2017

Programming Languages:
MATLAB
Spoken Languages:
English, Italian
Professional Interests:
Chemical Kinetics, Fluid Dynamics, Chemical Engineering, Heat and Mass Transfer

#### Statistics

All
• • • • • • • #### Content Feed

View by

Help with Graphing Difference Equations
n = 10; x = zeros(1,10); x(1) = 2; x(2) = 5; for idx = 3:n x(idx) = 5*x(idx-1)-6*x(idx-2); end semilogy(1:n,x)

1 day ago | 0

| accepted

Running for loop on an equation
I suspect you want to do something like this S = [ 1/2, 1/3, 1/2, 1/3, 1/3, 1/2;... 1/2, 0, 0, 0, 0, 0;... ...

1 day ago | 1

perimeter , area of circle
You just need to change the arguments in fprintf, since you must have accidentally swapped them. r = 1; P = pi * 2 * r; S = p...

2 days ago | 0

| accepted

How to apply stop time condition with multiple time intervals?
n=100; time = 0; for i = 1:n if sum(time) >= 35 fprintf('You run out of time') break end ...

3 days ago | 1

| accepted

How do I find a maximum within a certain range
Take this example dataset x = linspace(0,10,100); y = randi(10,1,100).*sin(x); plot(x,y) Now assume you want to find the max...

3 days ago | 0

Method of lines to solve unsteady state PFR with multiple reactions
I deleted many commented lines for clarity. clear, clc k0 = [0.1023 0.0011 0.0129 0.0000001 0.0008 0.000001 0.163 0.5 0.01 0...

4 days ago | 0

Help with time dependent ODE with a piece-wise component
I think that the problem is that you already defined the ODE system as an anonymous function, therefore there's no need to call ...

4 days ago | 0

| accepted

Solve 2nd order ODE using Euler Method
Hi Matt - a second order ODE can be decomposed into two first order ODEs. The secret is to set 2 variables y as The you have...

5 days ago | 0

convert char to double
var = '0'; out = str2double(var) whos out No need to use regular expressions at all, at least in Matlab.

5 days ago | 0

How to plot implicit ODE?
You need to solve it first. The ODE can be easily made explicit Then you can use any Matlab ODE solver (see this documentati...

5 days ago | 0

| accepted

How to plot this feather using a mathematical equation
There: clear,clc [x,y] = meshgrid(linspace(-0.4,0.9,1000),linspace(-1.5,1.1,1000)); xy = x-y/8-(1/8)*(y+1).^2; A = 1+1/...

6 days ago | 1

| accepted

How to solve systems of non linear equation of dimensions 100 using ode45 matalb.
If I were you, I would proceed substantially differently. Since you need to solve 100 equations, it is unthinkable to cde them ...

7 days ago | 2

| accepted

I would like to plot row 1 and row 3 of B as (x vs y plot graph).
I suspect you want something like this syms x y z for i=1:10 eqn1 = 2.*x + 2.*i.*y + z == 2; eqn2 = -x + y - i.*z ==...

7 days ago | 0

Force a starting point on exponential graph
See if this can help counts = [2423970,2171372,2065862,1830553,1100899,1037972,914015,752138,684123,606126]; normalized_counts...

8 days ago | 0

Sum groups of columns
Example M = rand(10,4600); n = 50; for idx = 1:size(M,2)/n S(idx) = sum(M(:,n*(idx-1)+1:n*idx),'all'); end size(S) ...

8 days ago | 0

I'm new in using Matlab and I'm very confused of how I would solve a derivative problem. Here is the equations and problem
I think what you have in your snapshot is already ok. You just need to save the functions in y1, y2 and y3. There are also some ...

8 days ago | 0

How to terminate the MATLAB code?
n = 100000000; a = 3.8; x(1) = 0.5; tic Time = 0; for i=1:n if Time > 5 fprintf('You run out of time') ...

9 days ago | 1

| accepted

Creating a matrix from one column
A = rand(1,1000)' M = repmat(A,[1 length(A)]) writematrix(M,'yourfile.csv')

9 days ago | 0

| accepted

How I can add an exponential fit to my cdfplot?
I guess you could do something like this load z x = linspace(0,1,length(z)); fitfunc = fittype(@(A,B,C,x) A*exp(B*x.^C)); ...

9 days ago | 0

| accepted

Use of 'ArrayValued' in Matlab numerical integration
Let's take a look at the error message f = @(x) 5; integral(f,0,2) The important line here is "Output of the function must be...

9 days ago | 0

| accepted

How to display only few values in a plot rather than for whole points in a big array?
n = 10; plot(x,T(:,1:n:end)); By increasing the value of n you decrease the number of lines shown in your plot.

9 days ago | 1

how to calculate biokinetic parameters?
You can use MatLab's fit function. Let me give you an example. Let us assume you have the following set of 50 experimental valu...

9 days ago | 0

| accepted

Expanding Array as following;
Try this A = [1 0 0 1]; B = repelem(A,3) For more info see the documentation for repelem.

10 days ago | 0

| accepted

Determine the magic sum from a magic square
I think you are inputting the magic square itself to your function, rather than its order. I believe the relation between a mag...

10 days ago | 0

| accepted

remove indent in command window
I do not know you can change the way the editor looks, but you can print the result so that it looks the way you want. b = -3; ...

10 days ago | 0

SIR model with recovered individuals may lose their immunity and become reinfected with the disease. But came with a failure about integration tolerances
I guess the problem is in this three lines S = beta*S*I + delta*R; I = beta*S*I - gamma*I; R = gamma*I - delta*I; where you ...

12 days ago | 1

Find Index of Max Value in One Matrix and find value of that index in another
You can do it this way force = rand(1,10); length = rand(1,10); [max_force,idx] = max(force); max_length = length(id...

15 days ago | 0

| accepted

max value in each row with its index
inputmatrix= [8.0000 0 7.3398 0 8.0000;... 1.6635 0.7103 3.2000 3.2000 3.2...

4 months ago | 0