Error using / Matrix dimensions must agree.

I have include my coding below which is showing operator error. I am trying to do 3D plotting of displacement.

P=[0:300]*10^3;
h=[5:20]*10^-6;
E=130*10^9;
a=150*10^-6;
for Z=0:100
    Z(:,Z)= (0.216*P.*a.^4)/(E.*h.^3);
end
surf(P,h,Z)

Answer 1

[David Hill]

[P,h]=meshgrid(linspace(0,300,300)*10^3,linspace(5,20,15)*10^-6);
E=130*10^9;
a=150*10^-6;
Z= (0.216*P*a^4)./(E*h.^3);
surf(P,h,Z);

[Kaushik Domadiya]

Thanks for answer!! But it is showing black surface.. Is there any idea how to make colour graph..