P=[0,0,2]; u=[1,0,-1]; Q=[0,0,1]; v=[0,1,1]; n=cross(u,v); a=n(1);b=n(2);c=n(3); d1=P*n'; d2=Q*n'; %tasot muotoa ax+by+cz=dj t=linspace(-3,3,100); [x,y]=meshgrid(t); z1=(d1-a*x-b*y)/c; z2=(d1-a*x-b*y)/c; surf(x,y,z1), hold on surf(x,y,z2) shading interp