fixing prob4

README.md

@@ -249,13 +249,9 @@ F = @(dx) 4*epsilon*6*(sigma^6./(xmin+dx).^7-2*sigma^12./(xmin+dx).^13)
 plot(dx_full,F(dx_full),dx,F_applied)
 
 [K,SSE_min] = fminsearch(@(K) sse_of_parabola(K,dx,F_applied),[1,1]);
-%fprintf('\n Nonlinear spring constants = K1=%1.2f nN/nm and K2=%1.2f nN/nm^2\n',K)
-%fprintf('The mininum sum of squares error = %1.2e \n',SSE_min)
 
 plot(dx,F_applied,'o',dx,K(1)*dx+1/2*K(2)*dx.^2)
 ```
 
 ![Plot showing LJ vs Bond Length](./figures/figure02.png)
 ![Plot showing Force vs dx](./figures/figure03.png)
-####Nonlinear spring constants = K1=%1.2f nN/nm and K2=%1.2f nN/nm^2
-####The mininum sum of squares error = %1.2e