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

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

17 hours 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 ...

21 hours 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...

1 day 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.

2 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...

2 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/...

3 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 ...

4 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 ==...

4 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...

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

5 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 ...

5 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') ...

6 days ago | 1

| accepted

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

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

6 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...

6 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.

6 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...

6 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.

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

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

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

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

11 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

Save for loop outputs into structure
Does allBIN(n).corr work?

4 months ago | 0

How to replace non consecutive value on a vector?
You could write a function that scans your array in search of the pattern you specify and replaces it with another. clear clc ...

4 months ago | 0

How to solve coupled partial differential equations with method of lines?
Hi @Ari Dillep Sai, take a look at the code below. The breakthrough now is found to be sometimes after 4000 s (or ~67 mins), so...

4 months ago | 0

| accepted

How to import multiple text files from multiple folders and take maximum from each text file
Hey @Nazanin Farsi, based on your replies, maybe the following code will work. n1 = 1:22; n2 = 0.1:0.1:1.3; max_value = zer...

4 months ago | 0

How to import multiple text files from multiple folders and take maximum from each text file
clear,clc n = 0.1:0.1:1.3; for k = 1:length(n) filename = [num2str(n(k)),'pga/AllMaxDrift.out']; data = read...

4 months ago | 0