Community Profile

photo

Alan Stevens


Last seen: Today Active since 2020

Statistics

  • 12 Month Streak
  • Guiding Light
  • Knowledgeable Level 4
  • First Answer

View badges

Content Feed

View by

Answered
Runge-Kutta integration of the Taylor-Maccoll eq
Shouldn't k1 = h*F1(t, y1(i), y2(i)); m1 = h*F1(t, y1(i), y2(i)); k2 = h*F2(t+0.5*h, y1(i)+0.5*h, y2(i...

5 Tage ago | 0

| accepted

Answered
How can I shade the area between a curve and a vertical line?
Try help fill or help patch Although they refer to polygons, you could easily represent your curve as a multi-segment polyg...

13 Tage ago | 0

Answered
Finding best parametric function estimation for ODE of first order
First plot the points, ui vs xi to see what sort of curve it might be. If it looks like it could be a polynomial (as in the cas...

13 Tage ago | 1

Answered
Impose a numerical limit on the value of a sum of two vectors
One possibility is min(max(a1),a1+a2)

15 Tage ago | 1

| accepted

Answered
How read Matrix Market data into MATLAB?
Do you mean like this: %%MatrixMarket matrix coordinate real symmetric % System: M x_dotdot + E x_dot + K x = B u % y...

16 Tage ago | 0

Answered
rounding the imaginary part of the number
Z = 0.0000e+00 + 7.8765e-04*1i; format shortG round(Z,4)

16 Tage ago | 1

| accepted

Answered
How can I model second order ODE with matrices and external forcing?
Like this perhaps (I've made up arbitrary data; you will obviously have to replace it with your own) x0 = [1; -1]; v0 = [0; 0]...

16 Tage ago | 0

Answered
Randomisation of ramberg osgood equation
You don't need a loop. Yu can just do deltaepsilon = r/E + 2*(r/(2*Kp)).^(1/np); where Kp = K' and np = n'. Note that you nee...

19 Tage ago | 0

Answered
How to solve a ode within a for loop?
Like this A = 0:0.05:0.35; tspan = 0: 0.1: 1; p32 = zeros(length(tspan),length(A)); r = 0.7; r0 = 2; p0 = [r; r0]; for k ...

21 Tage ago | 0

| accepted

Answered
How can I solve a system of equations with exponential terms in Matlab?
There is almost certainly no analytical solution here! For a numerical solution (given the values of the other parameters), try...

22 Tage ago | 0

Answered
To find integral of a function in loop
Try replacing q = integral(f,6e-7, x); with q(i) = integral(f,6e-7, x(i));

25 Tage ago | 0

| accepted

Answered
Numerical integration of the missile dynamic model
I don't know the actual problem, but if you look at your values of Y against time you will see that some of them are clearly div...

25 Tage ago | 0

| accepted

Answered
Monte Carlo simulation for probability of a multi part machine working
Here's one possibility: % each A part works 98% of the time. Each B works 91% of the time. % The machine has 12A and 7B parts....

26 Tage ago | 0

| accepted

Answered
Need help on using trapz on a definite integral
Replace Final = Nach + Vor; by Final(i) = Nach + Vor;

27 Tage ago | 0

Answered
Multiple parameter Non Linear Curve Fit- Error using ^
You probably need to use element by element multiplication and division etc in your definition of fun. I.e. use a.*b and a...

27 Tage ago | 0

| accepted

Answered
How to separate plot titles on single line? [solved]
Define the titles with an extra space at the end of each word.

27 Tage ago | 0

Answered
Plotting to test probability convergence
See my reply to your question at https://uk.mathworks.com/matlabcentral/answers/871048-how-to-do-integration-of-this-equation-an...

28 Tage ago | 0

| accepted

Answered
How to do integration of this equation and plot graph between T and y
Like this perhaps: u=0.020; x=0.025; z=0; t=5; eta= 0.6; P=200; C=500; ro= 1000; l=0.008; Alpha=17e-6; r = 0.01; % Gu...

28 Tage ago | 0

| accepted

Answered
Making relative histogram for data, and a convergence diagram to see if enough trials have been run for monte carlo simulation.
Something like the following? N = 1000; Throwsum = zeros(N,1); Relativefreq = zeros(N,1); count = 0; for i=1:1000 throwD...

28 Tage ago | 0

| accepted

Answered
Runge Kutta and Binary Search to analyze system of ODEs
At the moment you have dj/dt = -K*d0^2; which means rate of reduction of radius is constant. Since d0^2 is a constant this m...

29 Tage ago | 0

| accepted

Answered
Gaussian Curve Plotting Issues
If you have 10000 individuals with a probability of 0.013 of being in group A, you would expect your curve to peak close to 130 ...

30 Tage ago | 0

Answered
Solve equation with 2 unknowns using solve
There are an infinite number of solutions! Restrict the range. One possibility is as follows bfn = @(a) acos(1.6-cos(a)); fn...

etwa ein Monat ago | 1

Answered
Is this formula correct for simulate 3D Water rocket
Probably vy = v0*cos(alpha); should be more like vy = v0*cos(alpha)*cos(beta); where beta is an angle that will point the r...

etwa ein Monat ago | 0

Answered
Need to Produce a different size matrix on each iteration
Try the following wqe=[1;2;3;4;5;6;7;8;9;10;11;12;13;14;15]; kli=[51;52;53;54;55;56;57;58;59;60;61;62;63;64;65]; hsd=size(kli...

etwa ein Monat ago | 0

| accepted

Answered
Getting error using syms to implement power formula
Like this? % Arbitrary values - replace with your own fs = 1; x_h = @(n) 1/(n+10)^2; %Power formula syms N...

etwa ein Monat ago | 0

| accepted

Answered
Solving a system of 12 ode’s with discrete data using Euler method
Look up interp1 doc interp1

etwa ein Monat ago | 0

| accepted

Answered
How to save all variables in for loop
Try load(i,:)=com(i,1)*sin(com(i,2)*t);

etwa ein Monat ago | 0

Answered
How to handle a variable inside matrix without using syms toolbox?
Are you looking for somethig like this; f = @(w) [1,-0.1409*w^2;0,1]*[1,0;1/286.48,1]*[1,-0.05493*w^2;0,1]... *[1,0;1/1793...

etwa ein Monat ago | 0

| accepted

Answered
Issues in For loop
Have you intialised j (i.e. set j = 1;) before going into the for i = 1:38976 loop?

etwa ein Monat ago | 0

Answered
How can I randomize cube positions in matlab?
Like this? E=1; N=3; Nsamples=10; d=randi(N,1,Nsamples); s=randi(2,1,Nsample...

etwa ein Monat ago | 0

| accepted

Load more