Đề tài bài tập lớn môn giải tích 2 học kì 2 năm học 2016 2017 matlab Đề 3 và 5

y.m=doublem;n=doublen;l=doublel; %chuyen sang dang double.

ĐỀ TÀI BÀI TẬP LỚN MÔN GIẢI

TÍCH 2 HỌC KÌ 2 NĂM HỌC 2016 - 2017

ĐÁP ÁN

ĐỀ 3 Câu 1:

>> syms x y

[x,y1]=meshgrid(-1:.1:1,0:.1:3);z=1-x.^2;

for i=1:length(x)

for j=1:length(z)

if x(i,j).^2+z(i,j).^1 > 1

x(i,j)= NaN;

y0 =

0

>> [x1 y1

lamda]=solve(diff(L,x),diff(L,y),phi,x,y,lamda)x1 =

hold on;

plot3(5,0,25,'marker','*');

text(5,0,25,'CD(5,0,25)');hold on;

plot3(-5,0,25,'marker','*');text(-5,0,25,'CD(-5,0,25)');hold on;

plot3(0,5,-25,'marker','*');text(0,5,-25,'CT(0,5,-25)');hold on;

plot3(0,-5,-25,'marker','*');text(0,-5,-25,'CT(0,-5,-25)');hold on;

f =

(- x^2 + 2*x)^(1/2) - (x*(2*x - 2))/(2*(2*x - x^2)^(1/2))

Câu 5

>> syms n

>> symsum((3*n-1)/(n*4^n),1,inf)ans =

log(3/4) + 1

ĐỀ 5

Câu 1

syms x y z

f=80\(1+x^2+2*y^2+3*z^2);nabla=gradient(f,[x,y,z])

Câu 3

syms x y lamda real

f=input('nhap ham f(x,y)= ');f=sym(f);disp(' nhap toa do I(x0,y0) va R :')

<=> y.

m=double(m);n=double(n);l=double(l);

%chuyen sang dang double

hold on

t=linspace(0,2*pi,50);

x=r*cos(t)+x0;

y=r*sin(t)+y0;

a=max(A(:,1));

b=min(A(:,1));

if a==b

if a>subs(f,[x y],[x0+0.12 y0+0.12])

disp(['GTLN f(' num2str(A(1,2)) ',' num2str(A(1,3)) ')= ' num2str(A(1,1))]) else

disp(['GTNN f(' num2str(A(1,2)) ',' num2str(A(1,3)) ')= ' num2str(A(1,1))]) end

text(A(i,2),A(i,3),a+1,['(' num2str(A(i,2)) ','num2str(A(i,3)) ',' num2str(a) ')' ])

elseif A(i,1)==b

GTNN=[GTNN 'f(' num2str(A(i,2)) ','

num2str(A(i,3)) ')= '];

text(A(i,2),A(i,3),b-1,['(' num2str(A(i,2)) ',' num2str(A(i,3)) ',' num2str(b) ')' ])

[x, y, ~] = ellipsoid(0,0,0,2,2,2,100); % tao luoiz=sqrt(4-x.^2-y.^2);

z(imag(z)>0)=0;

set(surf(x,y,z),'facecolor','b','edgecolor','non','facealpha',.3)

hold on

% ve vector phap tai M0

x=x0;y=y0;z=z0;

n=eval(n);% the toa do diem vao vector

%bat dau ve vector phap