# Alan Stevens

Last seen: Today Active since 2020

Followers: 1   Following: 0

#### Feeds

View by

Implementation of the analytical expression for the magnetic field of a circular current loop and interpretation/representation of the results
When I run it I get different signs as well as slightly different magnitudes: %% Analytical calculation of the magnetic field o...

1 day ago | 0

How to remove gaps in timeseries plot to have a continuous graph
Include a NaN in the data list between each data set.

4 days ago | 0

what the reason behind roundoff error inside while loop
Compare the following dt = 0.1*16^-1; t = 0; while t<1 t = t+dt; end disp(t) t = 0; c = 1; while t<1 t = c...

12 days ago | 0

Something doesn't work for me in fit
Here's a fit using fminsearch (I just used a quick but coarse approach- the code can be made much cleaner!). v=[10 30 5...

17 days ago | 0

Plotting while skipping the middle point in a vector
Something like this? tic v1 = [1 2 3 4 5 8 9 10 20 21 22 30 31 32]; dv = diff(v1); ix = find(abs(dv)>1); v = []; s2 = 0; f...

1 month ago | 1

| accepted

help with nested loop and plot
Here's another way, still using nested loops: n= 101; x = linspace(-100, 100, n); y = linspace(-100, 100, n); [X, Y] = meshg...

1 month ago | 0

| accepted

How to integrate an implicit acceleration function to a fitting pair of graphs of velocity and displacement
Look up ode45 for the integration (for a fixed set of parameters) and fminsearch for fitting the parameters to known data.

2 months ago | 0

Brush: How can I change the font color of the displayed X and Y range?
Select the Property Inspector, select the text, choose a color. (The following code is just to get the image uploaded) I=imrea...

2 months ago | 0

how to use an integral function inside a for loop? i need to find the value of function that chages as i chages .then using this values i need to find another integral.
Are you looking for something like this? r0 =0.1; u =[1:7]'; Q =zeros(length(u),1); E =zeros(length(Q),1); Q = @(x)0.023*(r...

2 months ago | 0

Simpson integral - problem with writing the formula
Should be x(i) etc. not a(i). Also be careful with brackets. % Test 1 a = 0; b = 1; n = 10; f = @(x)1-x+x.^3; I = simpson(f...

2 months ago | 0

How do I apply iteration method to an integral?
Look up help on function fzero.

3 months ago | 1

| accepted

How can I plot a reflection angle of a projectile in a kinematics model? Overlay projectile motion plots?
Surely you only need to reverse the x-velocity? imperialBallVelocity = 52.20; release_angle = 42; release_height = 64.2; %Ki...

3 months ago | 0

Why am I receiving this error?
You have 4 layers layer_angles = [45, -45, -45, 45]; but data for only three mid_plane_strains = [0.05, 0.00, 0.00]; % Mid-p...

4 months ago | 0

Function g Returns Only 2 Values Instead of Array - Vectorized Operation Issue in MATLAB
You have g as a function of x2, but x2 doesn't appear on the right-hand side!

4 months ago | 0

| accepted

Solve for x with given equation
Look up the help on function fzero. And don't forget that Matlab doesn't recognise implicit multiplication. x(...) should...

4 months ago | 0

power method with rayleigh coeff
Make sure you call the function with the arguments in the same order as those defined in the function! A=[4 1 -1 0;1 3 -1 0;-1 ...

4 months ago | 0

| accepted

Fill area between two vertical curves
Like this: value = [186.7646 198.4115 191.1406 180.8430 175.7136 ... 167.7459 151.2865 144.5964 139.8148 139.4305 188.2865 204...

4 months ago | 2

| accepted

Simplificar el resultado como valor numerico
Why not simply: p1=[2;3] ; x0_1=@(theta) [cos(theta); sin(theta)]; % =[x1*x0; x1*y0] y0_1=@(theta)[-sin(theta); cos(theta)];...

5 months ago | 1

Optimizing code for HAC Weight matrix generation - How to be faster?
Try replacing the loops with something like this: i = (1:size(Weight_correlation,1))'; j = 1:size(Weight_correlation,2); c = ...

5 months ago | 0

How to a create a label at the bottom of a vertical line?
Like this? x = linspace(0,6,100); y = exp(x); plot(x,y) xline(4.5,'-',{'Acceptable','Limit'}, 'LabelVerticalAlignment', 'bot...

5 months ago | 0

| accepted

Secant method Nan Error
You haven't updated error in your while loop. Try including something like error = abs(X0-X1);

5 months ago | 0

| accepted

Function error in ODEs RK4 method, solving 6 unknowns
It works for me (I've commented out the fprint statements in the code below - but it works with them in!): global CDin psi thet...

5 months ago | 0

| accepted

I keep getting the error using '/'.
Try mu_ratio=(1-(vol_frac(c,1)/0.605).*((particle_ratio(c,1)).^1.2)).^(-1.5125); Notice the dots! .* .^ Difficult to be...

5 months ago | 0

Plateau followed by one phase decay
Like this? x = 0:0.5:20; % time in seconds Y0 = -0.6; % signal baseline value Plateau = -1; % singnal plateu after trigger/st...

5 months ago | 1

Radiation view factor using Monte Carlo
I think you probably need cos(theta), rather than theta, to be linearly spaced.

6 months ago | 0

How to fix a syntax in the oscillator's plot code?
Like this (both x_1 and x_2 are contained in y) T0 = 10; omega0 = 1; M = 1; F = 1; x_1 = 0; x_2 = 0; x_0 = [x_1 x_2]; ts...

6 months ago | 0

| accepted

ODE45has long runtime and graph will not plot
Here's my attempt to make sense of your equations tauspan = 0:100:10000; % p0 = [0.5 0.5 0 0]; % dp1dt = -rf + rr % dp2...

6 months ago | 0

Two graphs in the same plot, why?
Something like this data = xlsread('sunHoursInTrollhattan.xlsx'); mins = data(:,1)*60 + data(:,2); plot(1:12, mins,'o'),grid ...

6 months ago | 0

| accepted