Calc. Exp 1a and 1b

EXPERIMENT 1-A
1. Using MATLAB find the tangent to the curves y=√ x at x=4 at show graphically.
Solution:
The tangent to the line of the curve y=f(x) at the point P (a, f(a)) is the line
f ( a+h ) −f ( a)
through P with the slope m=f’(a)=lim
h→ 0 h
The following code illustrate the plotting of the tangent to the curve y= √ x at the
point x=4.
clear all
clc
syms x y
f = input (‘Enter the given function in variable x:’);
x0 = input (‘Enter the x-coordinate of the point:’);
y0 = subs (f, x, x0);
fx = diff (f, x);
m = subs (fx, x, x0);
tangent = y0+m*(x-x0);
t_line = y-tangent;
plotrange = [x0-3, x0+3];
ezplot (f, plotrange);
hold on;
ezplot (tangent, plotrange)
title (‘Tangent line plot’)
t = sprintf (‘The tangent to the curve y=%s at (%d, %d) is y=%s’, f, x0, y0, tangent’);
disp (t)
Input : Enter the given function in variable x: sqrt(x)

Enter the x-coordinate of the plot: 4
Output: The tangent to the curve y=x^ (1/2) at (4,2) is y=x/4 + 1
2. Using MATLAB find the tangent to the curves y=-sin(x/2) at the origin and
show graphically.
Solution:
The tangent to the line of the curve y=f(x) at the point P (a, f(a)) is the line
f ( a+h ) −f ( a)
through P with the slope m=f’(a)=lim
h→ 0 h
The following code illustrate the plotting of the tangent to the curve y=-sin(x/2)
clc
syms x y
f = input (‘Enter the given function in variable x:’);
x0 = input (‘Enter the x-coordinate of the point:’);
y0 = subs (f, x, x0);
fx = diff (f, x);
m = subs (fx, x, x0);
tangent = y0+m*(x-x0);
t_line = y-tangent;
plotrange = [x0-3, x0+3];
ezplot (f, plotrange);
hold on;
ezplot (tangent, plotrange)
title (‘Tangent line plot’)
t = sprintf (‘The tangent to the curve y=%s at (%d, %d) is y=%s’, f, x0, y0, tangent’);
disp (t)
Input : Enter the given function in variable x: -sin(x/2)

Enter the x-coordinate of the point : 0
Output : The tangent to the curve y= -sin(x/2) at (0,0) is y=-x/2
3 4
3. Verify Rolle’s theorem for the function (x +2) (x −3) in the interval [2,3]. Plot
the curve along with the secant joining the end points and the tangents at points
which satisfy Rolle’s theorem.
Solution:
Rolle’s theorem: Suppose that the function y  f (x) is continuous at every
point of the closed interval [a, b] and differentiable at every point in (a, b) and if
f (a)  f (b), then there is at least one number c in (a, b) at which f (c)  0.
The below code illustrates the verification of Rolle’s theorem for the function
3 4
f(x)=(x +2) (x −3) on the interval [-2, 3].
clear all;
clc;
syms x c;
f=input ('Enter the function: ');
I=input ('Enter the interval [a, b]: ');
a=I (1); b=I (2);
fa=subs (f, x, a); fb=subs (f, x, b);
df=diff (f, x); dfc=subs (df, x, c);
if fa==fb
c=solve(dfc);
count=0;
for i=1: length(c)
if c(i)>a && c(i)
count=count+1;
r(count)=c(i);
end
end
values=sprintf ('The values of c between %d and %d which satisfy Rolle’s theorem are x=', a, b);
disp(values)
disp(r)
else
disp('f(a) and f(b) are not equal, function does not satisfy conditions for Rolle’s theorem');
end
tval=subs (f, x, r);
xval=linspace (a, b, 100);
yval=subs (f, x, xval);
plot (xval, yval);
[p, q] =size(xval);
for i=1: length(tval)
hold on;
plot (xval, tval(i)*ones (p, q),'r');
hold on;
plot(r(i), tval(i),'ok');
end
hold off;
legend ('Function', 'Tangent');
title ('Demonstration of Rolles theorem');
Input : Enter the function: (x+2) ^3*(x-3) ^4

Enter the interval [a, b]: [-2, 3]
Output : The values of c between -2 and 3 which satisfy Rolles theorem are
x=1/7
4. Verify Lagrange’s mean value theorem for the function f(x)=x+e 3 x in the interval
[0,1]. Plot the curve along with the secant joining the end points and the tangents
at points which satisfy Lagrange’s mean value theorem.
Solution:
Suppose that the function y  f (x) is continuous at every point of the closed
interval [a, b] and differentiable at every point in (a, b), then there is at least
one number c in (a, b) so that
f ( b ) −f ( a)
f’(c)=
b−a
The below code illustrates the verification of Lagrange’s theorem for the function
f(x)=x+e 3 x on the interval [0,1]
clear all;
clc;
syms x c;
f=input ('Enter the function: ');
I=input ('Enter the interval [a, b]: ');
a=I (1); b=I (2);
fa=subs (f, x, a); fb=subs (f, x, b);
df=diff (f, x);
dfc=subs (df, x, c);
LM=dfc-(fb-fa)/(b-a);
c=solve (LM);
count=0;
for i=1: length(c)
if c(i)>a && c(i)
count=count+1;
tx (count)=c(i);
end
end
fprintf ('The values of c between %d and %d which satisfy LMVT are x=', a, b);
disp(double(tx))
xval=linspace (a, b,100);
yval=subs (f, x, xval);
m=subs (df, tx); % Slopes of tangents at the points between a and b.
ty=subs (f, x, tx);
plot (xval, yval);
hold on;
secant_slope=(fb-fa)/(b-a);
secant_line=fa+secant_slope*(x-a);
sx_val=xval;
sy_val=subs (secant_line, x, sx_val);
plot (sx_val, sy_val);
for i=1: length(tx)
txval=linspace(tx(i)-1, tx(i)+1,20);
t_line=ty(i)+m(i)*(x-tx(i));
tyval=subs (t_line, x, txval);
plot (txval, tyval,'k');
hold on
plot(tx(i), ty(i),'ok');
end
hold off;
grid on
legend ('Function', 'Secant', 'Tangents');
title ('Demonstration of LMVT');
Input : Enter the function: x+exp(3*x)

Enter the interval [a, b]: [0,1]
Output : The values of c between 0 and 1 which satisfy LMVT are x=0.6168
EXPERIMENT 1-B
1. Find the local and global maxima and minima for the function x 3-12x-5, x ∈
(-4,4).
Solution:
The following MATLAB code illustrates the evaluation and visualization of
global and local extrema of a function f (x) in an interval (a, b).
clear all
clc
syms x
f = input ('Enter the function f(x):');
I = input ('Enter the interval: ');
a=I (1); b=I (2);
df = diff (f, x);
ddf = diff (df, x); f = inline(vectorize(f));
df = inline(vectorize(df));
ddf = inline(vectorize(ddf));
range = linspace (a, b,100);
plot (range, f(range),'-b','LineWidth',2);
legstr = {'Function Plot'}; % Legend String
hold on;
guesses = linspace (a, b,5);
root = zeros(size(guesses));
for i=1: numel(guesses)
root(i) = fzero (df, guesses(i));
end
root = root (a <= root & root <=b);
root = unique(round(root,4));
plot (root, f(root),'ro','MarkerSize',10);
legstr = [legstr, {'Critical Points'}];
disp (['Critical Points of f(x) are: ', num2str(root)])
maxp = root(ddf(root) < 0);
if(numel(maxp) ~= 0)
disp (['Local maximum of f(x) occurs at: ', num2str(maxp)])
end
minp = root(ddf(root) > 0);
if(numel(minp) ~= 0)
disp (['Local minimum of f(x) occurs at: ', num2str(minp)])
end
fval = f(root);
gmax = root (fval == max(fval));
disp (['Global maximum of f(x) occurs at: ', num2str(gmax),' and its value is:', num2str(max(fval))]) plot (gmax,
f(gmax),'m+','MarkerSize',10);
legstr = [legstr, {'Global Maximum'}];
end
gmin = root (fval == min(fval));
disp (['Global minimum of f(x) occurs at: ', num2str(gmin),' and its value is: ', num2str(min(fval))]) plot (gmin,
f(gmin),'m*','MarkerSize',10);
legstr = [legstr, {'Global Minimum'}];
end
legend (legstr, 'Location', 'Best')
Input : Enter the function f(x): x^3-12x-5

Enter the interval: [-4,4]
Output : Critical points of f(x) are: -2,2
Local maximum of f(x) occurs at: -2
Local minimum of f(x) occurs at: 2
Global maximum of f(x) occurs at: -2 and its value is: 11
Global minimum of f(x) occurs at: 2 and its value is: -21
x
2. Find the global extrema of the function f(x)= x 2 e sinx - 3 on the interval [0,5].
x +1
Solution:
The following MATLAB code illustrates the evaluation and visualization of
global extrema of a function f (x) in an interval (a, b).
clear all
clc
syms x
f = input ('Enter the function f(x):');
I = input ('Enter the interval: ');
a=I (1); b=I (2);
df = diff (f, x);
ddf = diff (df, x); f = inline(vectorize(f));
df = inline(vectorize(df));
ddf = inline(vectorize(ddf));
range = linspace (a, b,100);
plot (range, f(range),'-b','LineWidth',2);
legstr = {'Function Plot'}; % Legend String
hold on;
guesses = linspace (a, b,5);
root = zeros(size(guesses));
for i=1: numel(guesses)
root(i) = fzero (df, guesses(i));
end
root = root (a <= root & root <=b);
root = unique(round(root,4));
plot (root, f(root),'ro','MarkerSize',10);
legstr = [legstr, {'Critical Points'}];
disp (['Critical Points of f(x) are: ', num2str(root)])
maxp = root(ddf(root) < 0);
disp (['Local maximum of f(x) occurs at: ', num2str(maxp)])
end
minp = root(ddf(root) > 0);
disp (['Local minimum of f(x) occurs at: ', num2str(minp)])
end
fval = f(root);
gmax = root (fval == max(fval));
disp (['Global maximum of f(x) occurs at: ', num2str(gmax),' and its value is:', num2str(max(fval))]) plot (gmax,
f(gmax),'m+','MarkerSize',10);
legstr = [legstr, {'Global Maximum'}];
end
gmin = root (fval == min(fval));
disp (['Global minimum of f(x) occurs at: ', num2str(gmin),' and its value is: ', num2str(min(fval))]) plot (gmin,
f(gmin),'m*','MarkerSize',10);
legstr = [legstr, {'Global Minimum'}];
end
legend (legstr, 'Location', 'Best')
Input : Enter the function: x^2*exp(sin(x))-x/x^3+1

Enter the interval: [0,5]
Output : Critical points of f(x) are: 0
Global maximum of f(x) occurs at: 0 and its value is: -∞

Calc. Exp 1a and 1b

Uploaded by

Copyright:

Available Formats

Calc. Exp 1a and 1b

Uploaded by

Document Information

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Calc. Exp 1a and 1b

Uploaded by

Copyright:

Available Formats

EXPERIMENT 1-A

Input : Enter the given function in variable x: sqrt(x)

Input : Enter the given function in variable x: -sin(x/2)

Input : Enter the function: (x+2) ^3*(x-3) ^4

Input : Enter the function: x+exp(3*x)

Input : Enter the function f(x): x^3-12x-5

Input : Enter the function: x^2*exp(sin(x))-x/x^3+1

You might also like