1 plot( [ solutionData(:).t ],[ solutionData(:).n3 ] );
3 plot( [ solutionData(:).t ],[ solutionDataCopy(:).n3 ] );
4 plot( [ solutionData(:).t ],
abs( [solutionDataCopy(:).n3] - [solutionData(:).n3] ) );
8 plot( [ solutionData(:).t ],[ solutionData(:).n4 ] );
10 plot( [ solutionData(:).t ],[ solutionDataCopy(:).n4 ] );
11 plot( [ solutionData(:).t ],
abs( [solutionDataCopy(:).n4] - [solutionData(:).n4] ) );
plot([solutionData(:).t],[solutionData(:).n3])
end if abs(real(dotprod))>rstoerst rstoerst