Ðề: tính tích phân kép bằng Matlab
nhập hàm số z=z(x,y).tính diện tích phần mặt cong z=z(x,y) nằm trong hình trụ
có đáy là hình tròn tâm I(a,b) bán kính R bất kì nhập tứ bàn phím. vẽ hình minh họa.
sau khi mày mò tất cả tài liệu có được, em đã viết ra đoạn code này . mong các sư huynh chỉ cho em lỗi sai hay thiếu sót chưa bao quát tất cả trường hợp. cám ơn các sư huynh :
function [kq]=dientichmatcong(I)
syms x y r phi
g=input('nhap ham so z=');
z=sym(g);
disp('nhap tam I(a,b) va ban kinh day hinh tru :');
a=input('nhap a');
b=input('nhap b');
R=input('nhap R');
f=(x-a)^2 + (y-b)^2
f= R^2
z'x=diff(z,x);
z'y=diff(z,y);
x=a + r*cos(phi) ;
y=b + r*sin(phi);
u=subs(z'x,x);
v=subs(z'y,y);
s=sqrt(1 + u^2 + v^2);
I=int( int(s*r,r,0,R) phi,0,2*phi );
disp('dien tich mat cong la :');
disp(I);
mong được các huynh chỉ bảo:22: