.

ltd13002 · Dec 14, 2017 · 3fbaa55 · 3fbaa55
1 parent be8747c
commit 3fbaa55
Show file tree

Hide file tree

Showing 2 changed files with 11 additions and 14 deletions.
diff --git a/Part G/part_g.asv b/Part G/part_g.asv
@@ -9,12 +9,9 @@ for i = 1:length(P)
 end
 clf
 setDefaults
-plot(wmax,P,'o')
-%P_out = interp1(wmax,P,linspace(.1,.3),'pchip');
-hold on 
-%plot(linspace(.1,.3),P_out)
-x = wmax;
-y = P;
-Z = x.^3;
-a = Z/
-hold off
+x = wmax';
+y = P';
+Z = [x.^0, x.^3];
+a = Z\y;
+x_fcn = linspace(min(x),max(x));
+plot(x,y,'o',x_fcn,a(1)+a(3)*x_fcn.^2)
diff --git a/Part G/part_g.m b/Part G/part_g.m
@@ -9,9 +9,9 @@ for i = 1:length(P)
 end
 clf
 setDefaults
-x = wmax;
+x = wmax';
 y = P';
-Z = [x.^3]';
-a = Z\y;
-x_fcn = linspace(min(x),max(x));
-plot(x,y,'o',x_fcn,a.*x_fcn)
+Z=x.^3;
+a=Z\y;
+x_fcn=linspace(min(x),max(x));
+plot(x,y,'o',x_fcn,a*x_fcn.^3)