Community Profile

photo

Are Mjaavatten


Last seen: Today Active since 2015

Retired physicist with background in the petroleum and process industry. My main interests are mathematical modelling, numerical methods, and thermodynamics.

Statistics

All
  • MATLAB Central Treasure Hunt Finisher
  • Thankful Level 1
  • Knowledgeable Level 3
  • Leader
  • Project Euler I
  • Personal Best Downloads Level 4
  • 5-Star Galaxy Level 5
  • Computational Geometry I Master
  • Indexing I Master
  • Sequences And Series I Master
  • Computational Geometry IV Master
  • Community Group Solver

View badges

Content Feed

Answered
Error in setdiff function
One workaround is to replace all NaNs with some spceific value that is not present in your data, say -9999: >> S1 = [1,2,3,NaN,...

3 månader ago | 0

Answered
How can I solve an ode where one of the variables is a matrix?
A(t) should be a function, not a vector. Below I define A as an anonymous function. I also use an anonymous function for your ...

3 månader ago | 0

| accepted

Answered
4D data /Index exceed array bounds
I think it may be better to let M_Val be a cell array of 5 ordinary arrays, where M_Val{k} is of dimension 134 by 2 by nsyll(k)....

3 månader ago | 0

Answered
4D data /Index exceed array bounds
I think you may have inadvertently created a cell array containing your whole data matrix in element 1. Something like this: >...

3 månader ago | 0

Answered
How to find index range of range of values in an array
i = find(x>=69 & x<=72);;

6 månader ago | 1

Answered
y coordinates of center of buoyancy
Is this what you want? x = [-8,-8,-6,0,6,8,8,-8]; y = [10,5,1,0,1,5,10,10]; ship=polyshape(x,y); figure(1) clf hold on p...

6 månader ago | 0

Answered
Is it possible to convert a cell array into individual variables on workspace?
Dynamically creating variables in Matlab is not encouraged. See Stephen Cobeldick's tutorial entry. A better idea is to create ...

7 månader ago | 0

Answered
How to plot a trinomial graph
fun = @(x) x^3; fplot(fun,[-3,3]) Alternatively: x = linspace(-3,3); plot(x,x.^3); % Note the use of .^, for element by ele...

7 månader ago | 0

Answered
Fitting a 3D gaussian to a 3D matrix
meshgrid arranges the xx and yy output matrices to fit the convention where the surface plots have x on the horizontal axis and ...

8 månader ago | 0

Answered
For loop assistance- incorrect data
You give DATA.s(3,1) the desired value in line 8 of exec2w. Then you overwrite it with the wrong value in line 22. I found deb...

10 månader ago | 0

Answered
How can resolve r0<=r0_max?
I interpret your question to mean that you want copertura{mm,nn}(k,1) to take the value 1 if r0{mm,nn}(k,1) <= r0_max(mm,nn), an...

10 månader ago | 0

Answered
Increasing precision of second order boundary value problem
Your asymptotic expansion solution seems to have zero slope at y = 0. This implies that the second derivative is zero and the s...

ungefär ett år ago | 0

Answered
Error using + Matrix dimensions must agree.
In line 43, x(:,i) is a column vector, while k1' is a row vector. Summing a column vector and a row vector is not allowed in ol...

ungefär ett år ago | 0

Answered
Golden Section Search algorithm, golden number why is 0.381
Let and be the end points of the current range. In golden ratio search we choose a new point such that the ratio of the two ...

mer än ett år ago | 1

Answered
Problem relating to mathematical expression
Ameer's answer is correct, of course. Here is another approach, without the elegant Matlab functions: N = 5; % Max window size...

mer än ett år ago | 0

Answered
Help text does not work for user-defined functions
Problem solved - finally! My Matlab 2019b and 2020a installations had several other odd behaviours in addition to the help prob...

mer än ett år ago | 3

| accepted

Answered
Solution for system of two non linear equation
You can rewrite you equations as (Nd-x-y)-a*(x*(x+y)) = 0 (Nd-x-y)*(No-y)-b*(y*(x-y)) = 0 where x = Nd, y = Nonegative, a = ...

mer än ett år ago | 1

| accepted

Answered
Problem of erf when using fmincon
Your problem has infinitely many solutions. Whenver x(1)-x(3) = 1 and x(2)-x(4) = 1, fun = 0, which is the minimum value possib...

mer än ett år ago | 0

Answered
how to Nonlinear equations?(two variables)
Although R1 and R2 are nonlinear finctions of t, your equations are linear in x. So you write the equations on the form A(t)*...

mer än ett år ago | 2

| accepted

Answered
How to count number of times a number appears in column B based on column A
letters = 'ABAACD'; numbers = [1,2,2,1,1,3]; uniqueletters = unique(letters); uniquenumbers = unique(numbers); n = length(un...

mer än ett år ago | 0

| accepted

Answered
Van der Waals equation - Newton's method
In your code, newton_ndim does not know the value of T and so cannot call fun(X,T). You can probably fix this by defining an in...

mer än ett år ago | 0

| accepted

Answered
How do I solve system of equations?
If I understand you correctly, your system consist of 128 eqations. a and b are arrays of 64 unknowns each. You must then solv...

mer än ett år ago | 0

| accepted

Answered
Join repeated continuous elements of a vector
in = [1 1 1 2 2 3 6 6 6 1 1 1 4 4 4 2 2 2 6 4 4]; i = 1; count = 1; out = []; a = in(1); while i < length(in) i = i+1; ...

mer än ett år ago | 0

Answered
Fsolve don't work good with trigonometric
Your function has several different zeros. Which one fsolve finds will depend on the starting point. I ran: f = 100; while no...

mer än ett år ago | 0

Question


Help text does not work for user-defined functions
I just installed R2019b and the help function does not work. Example: I created the addme function from the Add Help for Your P...

mer än ett år ago | 2 answers | 2

2

answers

Question


Is there a way to avoid repeating time-consuming calculations from odefun in the events function?
The right-hand side in my ODE is quite time-consuming. In the course of calculating dy/dt I calculate a variable u that I also u...

nästan 2 år ago | 2 answers | 0

2

answers

Answered
First order ODE parameter fit to dataset
The challenge is how to use your table of temperatures in your updateStates function. At a given time t, you must interpolate y...

ungefär 2 år ago | 0

| accepted

Answered
How to calculate a curve for different variables?
Here is one way to do it. I have reduced the number of elements in your Time and V vectors in order to make a nicer surface plo...

ungefär 2 år ago | 0

| accepted

Answered
datetime default stetings 12/24
Use HH instead of hh for the hours: >> datetime.setDefaultFormats('default','hh:mm:ss yyyy-MM-dd') >> t = datetime(2019,8,10,...

ungefär 2 år ago | 4

| accepted

Answered
How to read a file containing non numerical characters and commas?
The 22 lines of text are easily dealt with using the 'HeaderLines' option. Next you must convert your strings to doubles. Below...

ungefär 2 år ago | 0

| accepted

Load more