# How to code equation

15 Ansichten (letzte 30 Tage)
Cesar Cardenas am 24 Jan. 2023
Kommentiert: Cesar Cardenas am 3 Feb. 2023
Hello, just wondering if this equation could be coded like this:
Any feedback will be greatly appreciated. Thanks
m = 10e-3
a = 0.05
c = 0.1
x(t) = m/a e^-(c/a)t
##### 0 Kommentare-2 ältere Kommentare anzeigen-2 ältere Kommentare ausblenden

Melden Sie sich an, um zu kommentieren.

### Akzeptierte Antwort

Sulaymon Eshkabilov am 24 Jan. 2023
A couple of small errors have been corrected:
m = 10e-3;
a = 0.05;
c = 0.1;
x = @(t)(m/a)*exp(-(c/a)*t);
t = linspace(0,1.3);
plot(t, x(t))
xlabel('t')
ylabel('x(t)')
grid on
##### 1 Kommentar-1 ältere Kommentare anzeigen-1 ältere Kommentare ausblenden
Cesar Cardenas am 27 Jan. 2023
Thanks much for your help. I added another equation and would like to know how to plot both on the same graph? not sure about it. Thanks
m = 10e-3;
a = 0.05;
c = 0.1;
x = @(t)(m/a)*exp(-(c/a)*t);
x = @(t)(m/a)*1-exp(-(a/m)*t);
t = linspace(0,1.5);
plot(t, x(t))
xlabel('t')
ylabel('x(t)')
grid on

Melden Sie sich an, um zu kommentieren.

### Weitere Antworten (2)

Torsten am 24 Jan. 2023
Verschoben: Torsten am 24 Jan. 2023
m = 10e-3;
a = 0.05;
c = 0.1;
x = @(t) m/a*exp(-c/a*t);
t = 0:0.01:2;
plot(t,x(t))
##### 1 Kommentar-1 ältere Kommentare anzeigen-1 ältere Kommentare ausblenden
Cesar Cardenas am 27 Jan. 2023
Thanks much for your help. I added another equation and would like to know how to plot both on the same graph? not sure about it. Thanks
m = 10e-3;
a = 0.05;
c = 0.1;
x = @(t)(m/a)*exp(-(c/a)*t);
x = @(t)(m/a)*1-exp(-(a/m)*t);
t = linspace(0,1.5);
plot(t, x(t))
xlabel('t')
ylabel('x(t)')
grid on

Melden Sie sich an, um zu kommentieren.

Sulaymon Eshkabilov am 27 Jan. 2023
Put both functions in plot() command or use hold on:
m = 10e-3;
a = 0.05;
c = 0.1;
x1 = @(t)(m/a)*exp(-(c/a)*t);
x2 = @(t)(m/a)*1-exp(-(a/m)*t);
t = linspace(0,1.5);
plot(t, x1(t), 'r-', t, x2(t), 'b-.', 'linewidth', 2)
xlabel('\$t\$', 'Interpreter', 'latex')
ylabel('\$x(t)\$', 'Interpreter', 'latex')
legend('\$x_1(t)\$', '\$x_2(t)\$','Interpreter', 'latex', 'location', 'SE')
grid on
%% Alternative solution:
m = 10e-3;
a = 0.05;
c = 0.1;
x1 = @(t)(m/a)*exp(-(c/a)*t);
x2 = @(t)(m/a)*1-exp(-(a/m)*t);
figure
fplot(x1, [0, 1.5], 'r-', 'linewidth', 2)
hold on
fplot(x2, [0, 1.5], 'b-.', 'linewidth', 2)
xlabel('\$t\$', 'Interpreter', 'latex')
ylabel('\$x(t)\$', 'Interpreter', 'latex')
legend('\$x_1(t)\$', '\$x_2(t)\$','Interpreter', 'latex', 'location', 'SE')
grid on
##### 1 Kommentar-1 ältere Kommentare anzeigen-1 ältere Kommentare ausblenden
Cesar Cardenas am 3 Feb. 2023
Hello, thanks much for your response, really useful. Now, I'm trying to code these equations. This is my attempt. I would like to plot Cf vs x/L, but x/L as I have it shows me an error..the same with delta/x, not sure why? any help will be greatly appreciated. Thanks much.
Re = 0:500:1000;
x/L = 0:0.2:1;
Cf = 0.664*(Re)^-0.5;
delta/x = 5.0*(Re)^-0.5;

Melden Sie sich an, um zu kommentieren.

### Kategorien

Mehr zu Specifying Target for Graphics Output finden Sie in Help Center und File Exchange

### Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by