MATLAB: Plotting the following equation L=((N^2)muoAg1)/(x+((gR)/(2h)))

plot

I am new to matlab and need help plotting the following equation L=((N^2)*muo*Ag1)/(x+((g*R)/(2*h))) against x where x varies from xmin to xmax in steps of 0.0001m.

My parameters are as follows: R=0.065; h=0.0319; g=0.00055; xmin=0.0011; xmax=0.0187; N=2470; muo=4*pi*1e-7; Ag1=pi*R^2;

I have tried the following code with no luck: x=[xmin:0.0011:xmax] plot=(x*1000,L);xlabel=('position x (mm)'); x is multiplied by 1000 here to get value of x in mm rather than in m which are used in the equations.

Any help appreciated.

Regards

Mark

Best Answer

Your code works except for your not vectorising the division (replace â€˜/â€™ with â€˜./â€™) in the â€˜Lâ€™ calculation:

L=((N^2)*muo*Ag1)./(x+((g*R)/(2*h)));

Related Solutions

MATLAB: Integrals with syms over the variable x

If you have the symbolic toolbox this can be done:

syms E A0 L P c1 c2 c3 P x
e1=E*A0*x * (c1*(2 - (x/L)) + c2*(4*x - 2* (x^2) *(1/L)) + c3*(6*(x^2) - 3*(x^3)*(1/L)) - P);
F1=int(e1, L ,0) 
% Returns: 
% F1 =
% 
% -(A0*E*L^2*(27*c3*L^2 + 25*c2*L - 15*P + 20*c1))/30

The symbolic toolbox should also be able to solve the 3 eqs for c1, c2, and c3.

But considering that you have 3 polynomials in x, you should be able to calculate the integrals by hand easily.

HTH

MATLAB: How to plot a box on an existing plot

xmin=0;
xmax=1;
ymin=0;
ymax=1;
plot([xmin xmax xmax xmin xmin],[ymin ymin ymax ymax ymin])

rectangle('position',[xmin ymin xmin xmax])

Best Answer

Related Solutions

MATLAB: Integrals with syms over the variable x

MATLAB: How to plot a box on an existing plot

Related Question