Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                
Scribd
Upload
31 views

Matlab Assignment 3

nco ano paoubqjep

Uploaded by

gineethrs696
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
31 views

Matlab Assignment 3

nco ano paoubqjep

Uploaded by

gineethrs696
Copyright
© © All Rights Reserved
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 10

Calculus Assignment – 3

Name: - Gineeth.R

Registration Number: - 24BAI1017


%To evaluate the definite integrals , Riemann sums and compares it.

MATLAB CODE :

clc

clear all

syms x

f = input('Enter the function f(x):');

a = input('Enter the lower limit of x:');

b = input('Enter the upper limit of x:');

n = input('Number of intervals:');

value = 0;

dx = (b-a)/n;

for k=1:n

c = a+k*dx;

d = subs(f,x,c);

value = value + d;

end

value = dx*value

figure(1)
ezplot(f,[a,b])

z = int(f,a,b)

figure(2)

rsums(f,a,b)

Q1) F(x) = (x^3 cos(x/2)+1/2)( √(4-x^2)) from x=-2 to x=2?

Ans)

In the Command Window:

Enter the function f(x):

(x^3*cos(x/2)+1/2).*(sqrt(4-x^2))

Enter the lower limit of x:

-2

Enter the upper limit of x:

Number of intervals:

10

value =

(2*96^(1/2)*((8*cos(1/5))/125 + 1/2))/25 - (2*96^(1/2)*((8*cos(1/5))/125 -


1/2))/25 - (2*84^(1/2)*((64*cos(2/5))/125 - 1/2))/25 +
(2*84^(1/2)*((64*cos(2/5))/125 + 1/2))/25 + 38/25

z=

int((4 - x^2)^(1/2)*(x^3*cos(x/2) + 1/2), x, -2, 2)

In the Figure Window :


%To find the area of the regions enclosed by curves and visualize it. (Note:The
following code will work if the input functions are polynomials (not
trigonometric)).

MATLAB CODE :

clc

clear all

syms x y real

y1 = input('Enter the first curve:');

y2 = input('Enter the second curve:');

fg = figure;

ax = axes;

t = solve(y1-y2);

k = double(t)

n = length(k)

m1 = min(k)

m2 = max(k)

ez1 = ezplot(y1,[m1-1,m2+1]);

hold on

TA = 0

ez2 = ezplot(y2,[m1-1,m2+1]);

if n>2
for i=1:n-1

A = int(y1-y2,t(i),t(i+1))

TA = TA + abs(A)

x1 = linspace(k(i),k(i+1));

yy1 = subs(y1,x,x1);

yy2 = subs(y2,x,x1);

x1 = [x1,fliplr(x1)];

yy = [yy1,fliplr(yy2)];

fill(x1,yy,'g')

grid on

end

else

A = int(y1-y2,t(1),t(2))

TA = abs(A)

x1 = linspace(k(1),k(2));

yy1 = subs(y1,x,x1);

yy2 = subs(y2,x,x1);

x1 = [x1,fliplr(x1)];

yy = [yy1,fliplr(yy2)];

fill(x1,yy,'g')

end
Q2) Find the area of the regions enclosed by the curves y = -x^2+4x , y =
x^2-6x+5.

Ans)

In the Command Window :

Enter the first curve:

-x^2+4*x

Enter the second curve:

x^2-6*x+5

k=

0.5635

4.4365

n=

m1 =

0.5635

m2 =

4.4365

TA =

A=

5*15^(1/2)
TA =

5*15^(1/2)

In the Figure Window :-

Q3) Find the area of the regions enclosed by the curves y = 2x^5+3x^4-
5x^3+2x^2-x=1 , y = x^3-2x^2-5x+6.

Ans)

In the Command Window :

Enter the first curve:

2*x^5+3*x^4-5*x^3+2*x^2-x+1

Enter the second curve:

x^3-2*x^2-5*x+6

k=

0.8443

-0.9099
-2.7258

n=

m1 =

-2.7258

m2 =

0.8443

TA =

A=

2*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^2 - (4*root(z^5 +


(3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 1)^3)/3 - 2*root(z^5 + (3*z^4)/2 - 3*z^3 +
2*z^2 + 2*z - 5/2, z, 1)^2 + (3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z,
1)^4)/2 + (4*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^3)/3 -
(3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 1)^5)/5 - (3*root(z^5 +
(3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^4)/2 - root(z^5 + (3*z^4)/2 - 3*z^3 +
2*z^2 + 2*z - 5/2, z, 1)^6/3 + (3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2,
z, 2)^5)/5 + root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^6/3 +
5*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 1) - 5*root(z^5 + (3*z^4)/2 -
3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)

TA =

2*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^2 - (4*root(z^5 +


(3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 1)^3)/3 - 2*root(z^5 + (3*z^4)/2 - 3*z^3 +
2*z^2 + 2*z - 5/2, z, 1)^2 + (3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z,
1)^4)/2 + (4*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^3)/3 -
(3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 1)^5)/5 - (3*root(z^5 +
(3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^4)/2 - root(z^5 + (3*z^4)/2 - 3*z^3 +
2*z^2 + 2*z - 5/2, z, 1)^6/3 + (3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2,
z, 2)^5)/5 + root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^6/3 +
5*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 1) - 5*root(z^5 + (3*z^4)/2 -
3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)

A=

(3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^4)/2 - (4*root(z^5 +


(3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^3)/3 - 2*root(z^5 + (3*z^4)/2 - 3*z^3 +
2*z^2 + 2*z - 5/2, z, 2)^2 - (3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z,
2)^5)/5 + 2*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 5)^2 - root(z^5 +
(3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^6/3 + (4*root(z^5 + (3*z^4)/2 - 3*z^3 +
2*z^2 + 2*z - 5/2, z, 5)^3)/3 - (3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2,
z, 5)^4)/2 + (3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 5)^5)/5 +
root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 5)^6/3 + 5*root(z^5 +
(3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2) - 5*root(z^5 + (3*z^4)/2 - 3*z^3 +
2*z^2 + 2*z - 5/2, z, 5)

TA =

4*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^2 - (4*root(z^5 +


(3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 1)^3)/3 - 2*root(z^5 + (3*z^4)/2 - 3*z^3 +
2*z^2 + 2*z - 5/2, z, 1)^2 + (3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z,
1)^4)/2 + (8*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^3)/3 -
(3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 1)^5)/5 - 3*root(z^5 +
(3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^4 - root(z^5 + (3*z^4)/2 - 3*z^3 +
2*z^2 + 2*z - 5/2, z, 1)^6/3 + (6*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2,
z, 2)^5)/5 - 2*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 5)^2 +
(2*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2)^6)/3 - (4*root(z^5 +
(3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 5)^3)/3 + (3*root(z^5 + (3*z^4)/2 - 3*z^3
+ 2*z^2 + 2*z - 5/2, z, 5)^4)/2 - (3*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z -
5/2, z, 5)^5)/5 - root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 5)^6/3 +
5*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z - 5/2, z, 1) - 10*root(z^5 + (3*z^4)/2
- 3*z^3 + 2*z^2 + 2*z - 5/2, z, 2) + 5*root(z^5 + (3*z^4)/2 - 3*z^3 + 2*z^2 + 2*z -
5/2, z, 5)

In the Figure Window :

--- THE END ---

You might also like