surf
[X,Y] = meshgrid(1:0.5:10,1:20);Z =(Y*((3*(X.^2) + Y.^2) )/((X.^2 + Y.^2).^2));surf(X,Y,Z)Error using * Inner matrix dimensions must agree.Error in Mse302 (line 2)Z =(Y*((3*(X.^2) + Y.^2) )/((X.^2 + Y.^2).^2));[x y]=meshgrid(13:120,1:12);igx1=((x.*y.^3)/12)+(x.*y).*((((x.*y).*(150-(y./2))+113664+6000)./((x.*y)+1536+1200))-(150-(y./2))).^2;igx2=2097152+1536.*((((x.*y).*(150-(y./2))+113664+6000)./((x.*y)+1536+1200))-74).^2;igx3=10000+1200.*((((x.*y).*(150-(y./2))+113664+6000)./((x.*y)+1536+1200))-5).^2;z=igx1+igx2+igx3;surf(x,y,z)
Y0 = [2 5]; % intial conditions
xspan = [0 5]; % x-boundary values
% ode solver
[x, Y] = ode45(@rates, xspan, Y0);% note that Y(:,1) are the y values and Y(:,2) are the y' values
% Plot results
plot(x, Y) function dYdx = rates(x,Y) y = Y(1); v = Y(2); dydx = v; dvdx = -v + 2*y +x.^2 + x; dYdx = [dydx; dvdx];end
Best Answer