//example 3.8
//page 143
clc; funcprot(0);
// Initialization of Variable
k1=-9//dv/dx;
k2=16//du/dy
R=k1-k2;
if R~=0 then
disp("flow is rotational");
end