n=[1,2,3];d=2;
a=n(1);b=n(2);c=n(3);
t=linspace(-3,3,100);
[x,y]=meshgrid(t);

z1=(d-a*x-b*y)/c;%z ratkaistu yhtälöstä ax+by+cz=d

%figure('units','normalized','outerposition',[0 0 1 1])
%hold on
surf(x,y,z1),shading interp, hold on

piste1=[0,0,d/c];
piste2=piste1+n;
normaali=[piste1',piste2'];

plot3(piste1(1),piste1(2),piste1(3),'o')
plot3(piste2(1),piste2(2),piste2(3),'o')
plot3(normaali(1,:),normaali(2,:),normaali(3,:))
axis tight
axis equal