Community Profile

photo

Davide Masiello


Paul Scherrer Institut (PSI)

Last seen: Today |  Active since 2017

Contact

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

Statistics

All
  • 3 Month Streak
  • Knowledgeable Level 4
  • Thankful Level 2
  • Knowledgeable Level 3
  • Solver
  • First Answer
  • Thankful Level 1

View badges

Content Feed

View by

Answered
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)

circa un'ora ago | 0

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

circa 3 ore ago | 1

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

1 giorno ago | 0

| accepted

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

2 giorni ago | 1

| accepted

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

2 giorni ago | 0

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

3 giorni ago | 0

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

3 giorni ago | 0

| accepted

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

4 giorni ago | 0

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

4 giorni ago | 0

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

4 giorni ago | 0

| accepted

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

5 giorni ago | 1

| accepted

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

6 giorni ago | 2

| accepted

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

6 giorni ago | 0

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

7 giorni ago | 0

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

7 giorni ago | 0

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

7 giorni ago | 0

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

8 giorni ago | 1

| accepted

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

8 giorni ago | 0

| accepted

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

8 giorni ago | 0

| accepted

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

8 giorni ago | 0

| accepted

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

8 giorni ago | 1

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

8 giorni ago | 0

| accepted

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

9 giorni ago | 0

| accepted

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

9 giorni ago | 0

| accepted

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

9 giorni ago | 0

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

11 giorni ago | 1

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

14 giorni ago | 0

| accepted

Answered
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 mesi ago | 0

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

4 mesi ago | 0

Answered
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 mesi ago | 0

Load more